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SYNOPSIS 


R.S, Karmarkar 
Ph.D. Thesis 

Department of Civil Engineering 
Indian Institute of Technology.Kanpur, India 
(July 1982) 

SETTLEMENT AmiYSIS OF ANMJLAR FOOTINGS 


Footings in the form of an annular ring, with or 
without piles, are highly economical iwid well suited for 
tall isolated structures such as overhead water tanks, 
towers, chimney stacks .©%c* Such structures are usually 
centrally loaded. However forces such as those due to 
wind load may cause eccentric loading on the footing. 
Solutions available to-date for analysis of settlement 
of such footings are limited to those based on Elastic 
theory alone. Solutions for consolidation settlement 
of ring footings are not available and need to be developed. 

Except under a fully flexible footing, the contact 
pressures near the outer and the inner edges of a ring 
footing are very high even under small loads and are 
always likely to be beyond the yield capacity of the soil. 
Therefore, even under working load conditions, the soil 
in contact below the footing will be in plastic state over 
some zone. The case of centrally loaded rigid ring footing 



xix 


resting on a semi -inf inite mass of soil is taken up 
first. Using the method of numerical integration of 
available solutions from elasticity theory and the 
concept of limiting yield stress, an elasto-plastic 
analysis has been carried out to obtain an entire load- 
deformation response to failure. The available results 
from some experimental and Finite Element investigations 
have been compared with those o:*' the present v/ork and 
good agreement has boon found to exist. The validity of 
the proposed olasto-plastic analysis is thus established. 
The deformation is found to be nonlinear throughout and 
tends to infinity under high loads. 

The method has been extended to study the 
effect of embedment of a centrally loaded rigid ring 
footing on settlement response. Correction factors are 
given f02: different values of ratios of depth of embedment 
to outer radius, and inner radius to outer radius and 
for different values of Poisson's ratio. 

The settlement response of a centrally lofki.od 
rigid ring footing; resting at the surface of a compressible 
soil layer of finido depth is also studied. 



XX 


The case of axisymniotrically loaded ring footing 
of finite rigidity, defined in terms of dimensions of 
the footing and the relative elastic properties of the 
soil and the footing, has been investigated. 'The effect 
of finite rigidity on settlement and bending moments in 
the footing is studied under elasto-plastic soil conditions. 
The footinjr has been idealised as a plate and an integrated 
f ounclation-soil interaction analysis has been carried out. 

The plastic yielding of the soil is ■:^ound to have significant 
effect on the bonding momenta in the footing. A study of 
shear stresses in the footing has also been made. The 
e’^fects of embedment of the footing and tho finite layer 
of the soil are also studied. The positionimr of the 
applied load on the footing is another significant p"’rametcr 
studied in this analysis. 

To restrict the settlement, it is advantageous to 
provide a few short piles below va ring footing. In the 
present work, the axis.-'/mmetric case of centrally loaded 
rigid rih'j footing with one or more ooncenti'ic rows of 
incompressible piles has been investigated under clasto- 
plastic conditions of the soil. The technique of numerical 
integration of available elastic, solutio.no is again used 
coupled with plastic yield of soil in contact v\^ith the 
piles and the footing. For tho footing and tho piles. 



the load sharing and the extent of plastic yielding of 
soil in contact, depends on a n-umber of parameters including 
the imposed total load. Therefore, an integrated elasto- 
plastic analysis of the entire pile-raft system has been 
carried out. 

Eccentric loading on a ring footing, whether 
permanent or occassional, increases the possibility and 
extent of plastic yielding of soil especially under the 
edges. The present work includes elasto-plastic analysis 
of eccentrically loaded rigid ring footings either placed 
on or in semi -infinite or on finite soil. Under all loads, 
the settlement as well as rotation of the footing is found to 
vary nonlinearly with load. Under higher eccentricity 
and/or load, the settlement/rotation values are found to 
tend to Infinity, Just vas under all loads, it is 
necessary to consider the plastic yielding of the soil, 
it is found necessary to consider separation which occurs 
under some conditions, and which accentuates rotation. 

This aspect is also studied coupled with elasto-plastic 
conditions of soil. 

In addition to the final settlements, it is 
necessary to study the time rate of settlement of ring 
footings. The extreme cases of fully flexible and fully 
rigid ring footings undergoing axisymmetric consolidation 



have been investigated. Except for the Mandel-Oryer effect 
in the initial stages of consolidation, the Terzaghi-Rendulic 
theory predicts the overall consolidation response in a 
computationally economical manner with sufficient accuracy. 
Therefore, it has been adopted in the present study. The 
effects of various drainage boundary conditions, the radii 
ratio, the depth ratio, construction time and anisotropy 
with respect to permeability have been studied. In case of 
rigid footing, although the main aim was to study the 
consolidation settlements, a study of changes in contact 
pressures with time has also been included. 

During the course of the investigations on annular 
footings in this work, it has been observed that from the 
point of view of settlement at rail load levels, annular or 
ring footings are advantageous over circular footings and 
that it is essential to consider the plastic yielding of 
the soil at all load levels, and separation of footing 
from the soil wherever it occurs. 



CHAPTER 1 


INTROIDUGTIOH 


1 . 1 GE]\TERAI> 

A structure is only as good as its foundation. A 
proper choice of type of foundation suitable for the structure 
as well as for the foundation soil and a proper design of 
foundation are of utmost importance. A good design embodies 
a iudicious combination of requirements of safety and economy 
for ?/hich, an accurate geotechnical prediction of performance 
of the structure under all loads and at all times is 
essentially required. 

In general, the aspects of design of a foundation ares 

1) The stability of the structure with reference to 

its ultimate load hearing capacity( strength problem). 

2) The settlement (Deformation behaviour) 

(a) Maximum settlement 

(b) Differential settlement. 

3) Poundation^-soil interaction. 

The research work done in regard to the first aspect 
mentioned above has now reached a limiting stage where the 
methods of predicting the ultimate capacity are already well 
established. The second aspect in regard to deformation 



behaviour is however still a challenging research activity- 
A need for prediction of the entire load settlement response 
to failure has briskly activated the research in this field 
since late sixties of this century. (Davis, 1969, Roscoe, 1970) . 
The third aspect, although of a relatively recent origin, has 
brought in radical changes in the outlook on foundation design 
recognising the fact that the geotechnical and structural 
aspects of foundation design should not be uncoupled. 

The research work in regard to annular footings 
foundations is however, surprinsingly though, quite scanty. 

1 . 2 ANWJbA^ OR RIIG FOOTINGS AS POUNDAl’IONS 

For any isolated tall and heavy civil engineering as 
structure, a circular shape in plan is ideally suited 
structurally, functionally and aesthetically. Presently, 
all types of storage tanks, chimney stacks, silos and towers 
for several purposes such as aviation control, observatory, 
transmission , television etc. are the structures for vj’hich 
circular or annular footings are being or could be advantage- 
ously used. Annular footings are also suited where the columns 
or 'legs’ of a structure lie on a circle in plan. 

1.3 BIiASTO -PLASTIC AITALYSIS 

The importance of existance of very high stresses in 
soil near the edges of rizid footings has been recognized since 
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late thirties (Borowicka, 1936, 1939). Analytical and 
numerical solutions based on the theory of elasticity have 
shown that except for very flexible footings, the contact 
pressures near the edges of a footing are extremely high 
(theoretically, of infinite magnitude) . (Brown, 1968). 
Therefore, even under small loads, a portion of the soil 
is bound to yield plastically. This fact has been recognised 
and some values of limit inij edge pressure are suggested by 
various researchers , for example , Schultsse (1961), 
Zeevajert(l972). The exi stance of yieldin^j and the consequent 
limiting of edge contact stress and the dependence of pressure 
distribution on load level has been confirmed experimentally 
also, e.g, Selvadurai and Kempt home (1980), However, 
incorporation of these facts in the entire load -settlement 
analysis of annular footings has not been found in the 
present literature. The technique of using limiting stress 
concept has been used (e.g. Hain, 1975) but it has been 
reported to be not valid when large portions of soil yield 
under the ‘footing. At lower load levels, this approach 
does not give results far from elastic analysis. 

The major portion of work at present is based on 
elastic theory alone and that too pertaining to circular 
footings (and not annular footings), Eventhough elasto- 
plastic analysis is likely to give results close to those of 
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elastic analysis in the lower range of load levels, it is 
necessary to find the extent of agreement quantitatively. 
Further, elasto-plastic analysis becomes necessary and 
important at working loads and in the nearby range of 
loading in which thd actual conditions may exist, The 
load-settlement behaviour, if predicted sufficiently 
accurately, in the range of loading between working load 
and failure, is of utmost importance, to visualise and 
predict the risk involved. All such studies are possible 
only if an elvasto -plastic analysis is employed. It ought 
to be mentioned here that in the present age of competetive 
designs, and low factors of safety, the prediction of the 
entire load settlement response, including situations of 
exi stance of yielding of soil over considerable area of 
the footing, is all the more necessary. In several situations 
the allowable bearing pressure is to be read off against 
allowable settlement for which the entire load-settlement 
curve should be predicted realistically. It is therefore 
appropriate to use a reasonably accurate simplified elasto- 
plastic analysis than relying on rigorous elastic analysis 
by itself, 

1.4 TK15 OHOIGE OF mSEOD OF MJMERICAL AimLYSIS 

For elasto-plastic analysis of the type of problems 
tackled in the present work, the suitable methods are 
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approximate numerical techniques, viz., (a) the Finite 
Element Method (F33M) and (b) the method of numerical 
integration of elastic solutions with the simplified 
concept of limiting the stress to yield stress value of 
soil. The Bounday Element Method (BEM) is not considered 
here because of the complexity involved through iterative 
procedure and the requirement of availability of elastic 
solution involving'!; body forces. 

Assuming? that appropriate stress-strain relation 
and elasto-plast ic model are available and that sufficient 
data of fixed or variable soil parameters, stress paths 
etc, is properly incorporated, a PEM elasto-plast ic solution 
of a three dimensional problem is computationally very 
costly. If a parameteric study or design data is the aim, 

PEM may turn out to be computationally prohibitive. At 
the same time, PEM is a powerful and perhaps the only 
tool to deal with a specific complex problem and should 
be used in such cases when the need arises. 

There are several situations wherein elastic 
solutions are readily available. It is desirable to take 
full advantage of these and the method of numerical inte- 
gration is the one which is based on this. This method 
avoids soil discretization,, which is the main cause for 
using up computer memory and time to a great degree in 
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methods like PEM. Sven in some nonhomogeneous conditions, 
elastic solutions based on some average or equivalent 
values are expected to give fairly good results. 

As far as introducing plastic yield of soil is 
concerned, the method of first finding the stresses from 
elastic theory and then limiting the stresses to yield 
stress value of soil causing redistribution of stresses 
has been already in use. This method involves a computation 
effort which, in comparison to PEM employing even simplified 
elasto-plastic analysis, could be termed as negligible. 

However, as pointed out earlier, the present methods of 
this type do not cover the entire ran^je of load-deformation 
response. 

It may be added here that if PEM is to be used for 
axi symmetric problems involving solution upto collapse load, 
enormous computational cost will be involved since fifteen- 
noded elements would have to be used (Sloan and Randolph, 1982). 

1.5 0RGA.NISATI0N OP THE THESIS 

In the preceeding sections of this chapter, a brief 
introduction has been given regarding the importance of, 
the problem taken up in this thesis and the need to 
employ an appropriate method of analysis and a method of 
solution with a wide parameterio coverage and yielding fairly 
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accurate results at reasonable computational costs. 

In the following paragraphs, a brief description 
of the contents of the succeeding chapters is given. 

Chapter 2 is devoted to a broad review of literature 
available in the field relevant to the present work. 

Chapter 3 deals with the load settlement response 
of a smooth perfectly rigid annular centrally loaded footing. 
(Circular footing is dealt with throughout the thesis as a 
particular case of annulcsr footing). Along with the general 
and particular assumptions, some remarks about the implied 
limitations are given. The analysis presented is based on 
the use, to start with, of available elastic solutions 
for settlements due to a point load. In general, the 
solutions are numerically integrated using appropriate 
discretization of the footing. The following types of 
problems are handled; 

Problem tjrpe A - Footing resting on the surface of a 

serai -infinite soil. 

Problem type B - Footing embedded in semi -infinite soil. 
Problem type C - Footing resting on the surface of a soil 

layer of finite thickness, underlain by 
rough rigid base. 
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For dealing with, elasto -plastic analysis, the 
concept Qf limiting stress has been used. This along with 
partitioning and modification of influence coefficient 
matrix have been employed to obtain the load settlement 
response in the entire range of loading upto failure. 

Chapter 4 deals with the same cases as of Chapter 3 
but taking into account the finite rigidity of the footing. 
All the studies are made in terms of the finite rigidity 
of the footing relative to soil. A broad study of bending 
moments and shear forces (maximum values in particular) 
induced in the footing; has been included in this chapter. 

The footing is idealised as thin plate and the interaction 
an^ysis is carried out by combining the inverted soil 
flexibility matrix with the plate rigidity matrix. The 
configuration of loading on the footing plays an important 
role and its effects have been studied broadly in this 
chapter. In some cases, like a centrally loaded relatively 
flexible footing, separation occurs between some portion 
of the footing and the supporting soil. This has been 
taken into account in the analysis contained in this chapter. 

The effect of supplementing a centrally loaded 
rigid ring footing resting on semi-infinite soil by piles 
has been studied in Chapter 5. The piles are assumed 
incompressible. Along with the footing, the piles are also 
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discretiaed and an integrated elasto-plastic analysis is 
carried out in this chapter. Some results for a single 
floating pile are included for the purpose of comparison 
with available solutions. The effects of length, number 
and positioning of piles are broadly studied and discussed 
in this chapter. 

Chapter 6 deals with eccentrically loaded rigid 
ring footings. The three t3rpes of problems dealt with 
in Chapters 3 and 4 are included for study in this chapter 
also. The problem of seperation or lift-off, which crops 
up due to high eccentricities along with plastic yielding 
of soil has been dealt with in this chapter. The method of 
analysis is similar to that of Chapter 3 or 4, but the 
method of solution has been modified slightly to achieve 
economy in computations. 

Chapter 7 deals with the aspect of consolidation 
settlement under perfectly flexible and perfectly rigid 
annular footings resting on a finite layer of soil. The 
earlier portion of this chapter deals with flexible footings, 
Terzaghi-Rendulic theory has been used, on the lines of 
Madhav and Vitkar (1981) to solve this problem. In the 
latter portion of this chapter, the problem of consolidation 
of rigid footings is handled. The changes in contact pressure 
distribution, which are obtained as a by -product of the 
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solution, axe also studied. The construction of footing 
does take sufficient time and this aspect has been included 
in this chapter assuming linear build-up of load. Several 
tjrpes of soil may have higher permeability in horizontal 
direction than in the vertical direction. The anisotropy 
on account of permeability has also been included in the 
study. 

In all the above mentioned chapters, comparison 
has been included, with available solutions. Also charts 
have been presented in all the chapteis which could be 
directly used for design. 

The concluding chapter (Chapter 8) contains the 
summary of the work covered in this thesis, the conclusions 
drawn on the basis of the study taken up and some suggestions 
for further study. 



CHAPTER 2 


literature review 


2 . 1 GEUBRAL 

Literature in the field of integrated soil -foundation- 
super structure analysis for an accurate prediction of behaTiour 
with the inclusion of time effects also, is scanty, mainly 
because incorporation of details -especially of superstructure- 
amounts to loss of generality. Even in the field of soil- 
foundation interaction, the earlier work proceeded along one 
detail or the another making highly simplifying assumptions 
about the rest of the system. The 'gaps' created in the 
research work of this type-although there is a brisk research 
activity since late sixties to fill them up-still continue 
to exist to a certain extent. A brief discussion about these 
gaps in the field relevent to the present work is presented 
in the following few paragraphs. 

An extract from the lectures by Davis (1969) on the 
analysis of the settlement of foundations is quoted below, 

'The use of elastic theory for the prediction of 
settlement becomes inappropriate once the load to cause 
local yield within the soil mass has been passed. Further 
increase in load produces growth of plastic regions until 
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they are sufficiently extensive to permit continuous 
settlement, at which point the ultimate bearing capacity 
is reached. During the elasto-plastic phase between wholly 
elastic conditions and full plastic failure, the stress 
distribution departs increasingly from the elastic distribution. 
In the plastic regions, the strains at any point are no 
longer determined solely by the stresses at the point, so 
that even if the stress distribution can be approximately 
estimated in an elasto-plastic situation, the settlement 
can no longer be predicted by analytical or numerical 

integration of the stresses The dilemma is illustrated 

in f if^ire I (See Pig. 2.1). 

The exist ance of a similar situation in the case of 
the basic stress-strain models for soil is expressed 
widely by Roscoe (1970) in the tenth Rankine lecture, 

another 'gap ' manifested in spite of the tremendous 
research activity is , a means to obtain at nonprohibitive 
computational costs, a realistic end product in terms of 
stresses and deformation responses of foundations using 
the best available soil and mathematical models. Por 
example, Das and Gangopadhyay (1978) have shown in a particular 
study that even though a bilinear elastic-plastic finite 
element analysis is adopted, it requires different sets of 
assumptions and computer runs to obtain realistic values of 
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stresses and deformations. 

Yet another *gap ' is evident even to-date through the 
detached roles being played by the structural and geotechnical 
aspects of the design of foundations. The use of Vfinkler 
type of model of soil reaction found even to-date i?i literature 
is filso indicative of this. 

The separate facets of the relevant research work, with 
all the gaps, are still important in their own v-axj, For 

Ccjvvji ti. 

example, elastic solutions are still the basis of mo st of - - t j^e 

research activity. Hence the various facets are included, 
atleast for the purpose of citing the examples, in the brief 
review which follows. The review is ,in general , restricted to 
circular or annular footings, 

2,2 AHAIYTICAL WORiC GOiJHElGTED WITH SETTLEIffiiNT OF GIRGUIAR OR 

MIMLAR FOOTINGS IH GOHTAGT WITH HOMOGEHEOUS ISOTROPIC MASS 
OF SOIL-WITH HO OOHSIDERATIOH OF TIDIS EFFECTS 

2.2.1 Linear Elastic Solutions 

Under certain simplif^ring assumptions, a circular ( or 
even annular) loading configuration is relatively easy for 
mathematical treatment, e.g, the work of Volterra (1953) on 
annular footing resting on Winkler type of foundation. 

Leaving aside the analyses based on Winkler type of soil model, 
there are several analytical elastic solutions for fully 



flexible or fully rigid footings. The basis for most of 
these solutions is either 3onssinesq*s (1385) or Miudlin’s 
( 1936 ) solution for point load. In this respect the 
tabulations by Ahivin and Ulery ( 1962 ) are useful. The 
solution for annular flexible footing resting on soil of 
finite thickness is given by Milovic (1973). The solution 
for flexible circular footing embedded in semi -inf inite 
soil based on Mindlin's solution is given by Hishida (1966) 
and that for flexible annular footing on finite layer of 
soil is given by Zinovev (1979). Stress distribution under 
embedded circular footing using a combination of Westergaad- 
Mindlin solutions is studied by Babu Shankar (1977). 

The case of axially loaded rigid circular raft in 
frictionless contact with a homogeneous isotropic elastic 
half -space has been considered by Sneddon (1946), Schiffman 
and Aggarwala (1961), Chan and Ghaung (1979), and the 
cor responding.; case of adhesive contact has been studied by 
Mossakovski (1954), Keer (1967), Spence (1968), Popova (1972) 
and Bgorov et al. (1973). The effect of embedment of rigid 
circular area in semi-infinite soil is investigated by 
Butterfield and Banerjee (1971) and that on the surface 
of a soil of finite thickness is studied by Bgorov and 
Serebrjanyi ( 1963 ) and Poulos (1968), The effect of 
embedment of a rigid circular area in a soil of finite 
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thickness is investigated by Banerjee (1971 ), 

Annular footings are considered by Egorov (1965) 
for smooth contact and by Shibuya et al, (1974, 1975) for 
smooth as well as adhesive contact. Annular footings are 
also studied by Valov (1968) and Dhaliwal and Singh (1977). 
The effect of finite depth is empirically studied by Egorov 
et al. (1979). A modified version of approximation suggested 
by Poulos and Davis (1975) is used by Madhav (1980) to extend 
the Milovic's solution for flexible rings to smooth rigid 
rings resting on a finite layer of soil. 

The flexure of a thick circular raft resting on 
elastic medium is studied by Selvadurai (1980). A method 
of settlement calculation considering statically and 
kinematically admissible solutions for uniformly loaded 
circular area or vertically loaded rigid circular footiiag 
is suggested by Denvar (1981). 

Moment loadin^j on smooth rigid circular footing 
is studied by Borowicka (1943) and Sneddon (1946), for semi- 
infinite soil. Approximate solution for the corresponding 
case for circular footin^j resting on a soil layer of finite 
thickness is given by Egorov and Nichiporovich (1961). The 
case of eccentric loading on smooth rigid annular footiiogs is 
investigated by Kara et al. (1975) and by Egorov et al.(l979, 
1980). Eccentric loading of embedded rigid circular footing 



is considered by Selvadurai (1980) and a study of stress 
and strain state of a subgrade under eccentric loading is 
done by Egorov et al. (1930). Wilson (1965) has suggested 
a novel method of converting as 3 ni]metric problems into 
axisymmetric problems by expressing the loading- as a Courier 
series. 

The above mentioned work, although it pertains to 
fully flexible or fully rigid footings, is good enough to 
predict the overall elastic settlements. However, the 
differential settlements and bending moments and shear forces 
in a footing are highly sensitive to the relative rigidity 
(or flexibility) of the footing and to the loading confi- 
guration. Some significant works connected with circular 
footini; 5 e of finite rigidity are mentioned below, 

A significant work in this field (circular footings 
of finite rigidity) goes back to as early as 1939 . Borowicka 
( 1939 ) gave solutions for uniformly loaded and centrally 
point loaded circular footings. Except for a few works 
e.g. Ishkova (1951), the work in this field appears to have 
taken impetous only from late sixties. Analytical solutions 
are given by Brown (1969) to the problem of uniformly loaded 
circular footing resting on semi-infinite soil or a finite 
layer of soil. Analytical solutions are also advanced by 
Hooper (1976) for parabolic adhesive loading and by 
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Selvadurai (1979) for unifoim loading using variational 
approach. Some numerical solutions in connection with 
circular rafts are those given by Barden (1965), Smith(l970), 
Zbirohowski-Koscia and G-unasekera (1970), Hooper (1974) and 
Ohakravorty and Ghosh (1975). Various types of numerical 
techniques were used such as finite difference (e.g, 
Ghakravorty and Ghosh), finite element (e.g. Smith), 
iterative (e.g. Zbirohowslct-Koscia and Gunasekara) etc. 

Wood (1977) has presented an approximate analysis, 
Ghakravorty and Ghosh have considered as’'/mmetrical loading 
also. 

Lastly, some of the important elastic solutions, 
although not fully covered under the topic of this section 
(such as pertaining to non-homogeneity, anisotropy etc.) but 
pertaining to circular shapes or footings are cited as 
follows. England (1962), Gerrard and Harrison (1970), Besai 
and Reese (1970), Chen and Engle (1972), Gerrard and Wardle 
(1973), Gainrier and Ghristian (1973), Hooper (1975), Boswell 
and Scott (1975), Buragohain and Shah (1977), Varma et al, 
(1977), Wardle (1977), Hanuska et al. (1978), Brown and 
Gibson (1979) and Gaaetas (1981, 1982). Circular footin^^s 
with loads which are transverse, torsional etc, are not 
included in the review. 
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Many important elastic solutions are recorded by 
G-iroud (1972) and Poulos and Davis (1974). A detailed 
discussion about the analytical approaches and methods to 
solve the interaction problem of 'plate on continuum' is 
presented by Selvadurai (1979). 

2.2,2 Uon-l inear it y as a Consequence of Loss of Contact 

There are certain situations in which an elastic 
analysis would indicate the presence of tensile contact 
stresses between the footing and the soil. Since real soil 
is known to be extremely weak in tension, the existance 
of such tensile stresses is unrealistic, and loss of contact 
or separation would occur. This can occur due to moment or 
eccentric loading or even due to a loading near the centre 
of a relatively flexible circular footing. The conventional 
assumption of linear variation of contact pressure can give 
a rough estimate of the cut-off area. Dave (1977) studied 
this aspect for annular footing on the assumption of linear 
variation of contact pressure and concluded that annular 
footings are better than circular ones, 

Gheung and Nag (1968) developed a simple but effective 
procedure to deal with the problem of separation of contact 
surfaces. It is outlined below, 

(i) Obtain a linear elastic solution. 
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(ii) If all the contact pressures are compressiTe the 
problem is solved. Otherwise proceed to next step. 

(iii) Find out the nodes which are associated with tensile 
contact pressures and make the corresponding rows and 
columns in the original flexibility matrix zero. 

,.(iT) ' Invert the new flexibility matrix for the foundation 

, to obtain a new stiffness matrix, 

\ * ' 

(v) ^0,Qmbine this with the plate (footing) stiffness matrix 
by .assuming? compatibility at the nodes that have 
compressive contact stresses and go to step (ii), 

Svec ( 1979 ) also used a method which is basically 
the same as that of Oheung and Hag, except that it replaces 
inversion of the full new foundation flexibility matrix by 
inversion of a small matrix related to the region of loss of 
contact only through a procedure of partitioning.. 

G-ladwell and Iyer (1974) considered the loss of 
contact between a centrally loaded circular plate and the 
elastic half -space. Solutions' using finite difference are 
obtained by Banerjee and Jankov (1975) for asymmetrically 
loaded circular rafts, considering separation. 

In FEM, tension cut-off has been dealt with on the 
basis of 'no-tension' analysis developed by Zienkiewicz et al. 
(1968), e.g, G-ermanis and Yalliappan (1975). 
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2.2,3 Ufon-linearity in Load Settlement Response Under 
Increasing Load 

Bro-wn's (1968) discussion on the effect of local 
bearing failure on behaviour of rigid circular rafts is 
significant of the awareness towards the need to study- 
stages inbetweenthe linear elastic settlements and ultimate 
bearing capacity. (It is now a recognised fact that the 
whole response is nonlinear). Roscoe’s (1970) tenth Rankine 
lecture also highlighted the need to predict load deformation 
characteristics in the entire load range. The research work 
in this aspect can be deemed to have started around this 
period and also it was accelerated due to the employment 
of PBM. 

Desai and Reese (1970) employed PEM along with 
non-linear representation of stress-strain relationship to 
obtain a very realistic non-linear load-settlement response 
of circular bearing plate. The literature related to 
general techniques and analyses incorporating non-linearity 
in PBM can be cited in the works of Oden (1972), Desai 
and Abel (1972), Sienkiewics (1977) and Desai and Christian 
(1977). A brief discussion under classified headings is 
given below in connection with non-linearity in general. 
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(a) Pseudo -Elastic Methods of Non-linearity 

In these raethods the non-linearity is introduced 
through changes in elastic parameters of the soil with no 
consideration of onset of plastic yield. The most 
straight-forward raethods are based on the shape of stress- 
strain curve of the soil. Bilinear is the simplest form. 

An approach in which the second portion of the bilinear 
model (or the last one in multilinear model) represents 
a very small slope is quite popular since it gives elasto- 
plastic like response (Pig. 2.2) , One such work is by 
Dunlop and Duncan (1970) in which the stresses are checked 
against shear strength and if they exceed the strength, 
the corresponding elements are assigned a reduced value of 
modulus. Hyperbolic law (Pig. 2.2) based onKodner's (1963) 
work is one of the most popular mathematical forms which 
follow closely the usual stress-strain curves. The formulation 
for TJ-se in PEM is given by Duncan et al. (1970), 

Ramb erg-0 sgood (1943) presented a three -parameter mathematical 
function. Wong (1971) suggested a parabolic function, Desai 
( 1972 ) developed an approach in which the stress-strain data 
is simulated using cubic and bicubic spline functions (Pig, 2, 2) . 

There are numerous ways in which the elastic 
parameters can be varied, hence a large variety of solutions 
are also seen in the literature, Clough and Woodward (1967) 
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suggested, in the study of embankments, that the deformation 
modulus (=E/2(l+ v))should be changed keeping the bulk 
modulus M|^(=E/2(l+v )(l-.2v )) constant. G-iri javallabhan and 
Reese (1968) selected the value of E depending on the state 
of strains, Huang (1969) assumed that E should be changed 
with the first and second invariants of stress for sand 
and clay respectively. Das and Gangopadhyay (1978) even 
suggested, for the solution of the same problem, two 
different methods of changing elastic parameters, They 
reported that in the bilinear elastic portion, B should be 
changed keeping v constant to obtain stresses correctly, 
whereas both B and v should be changed keeping bulk 
modulus constant to obtain def omations correctly. 

A very simple and interesting approach has been 
adopted recently by Bobe and Pietsch (1981). They have 
suggested that based on a new strength theory by Szabo 
( 1970 / 71 ), a variable modulus V = E(l-qjjj/q^) should be 
used which directly depends on the load level qjn/9f » where 
qjjj is the mean applied pressure on the bearing plate or 
foundation and is the ultimate bearing pressure (which 
depends on geometry of the footing). However this modulus 
is to be used to find out directly the settlement at any 
load level, and the contact pressures etc, are to be worked 
back. The method is suggested to be used for shallow 
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foundations only and expression for for rectangular 
plate or footing is given by them. 

The above method by Bobe and Pietsch yields the 
following results which are qualitatively realistic. 

(1) The load -settlement response is non-linear 
throughout. 

(2) The settlement tends to infinity as load approaches 
ultimate load. 

(3) The contact pressure distribution (worked out) shows 
no peaks of extremely high pressures near the edges 
as predicted by elastic theory* 

They have claimed that their method gives results 
tallying with some experimental results of other investigators* 

In their method, it is found that although non- 
linearity in settlement response is actually a consequence 
of and should depend on the contact pressures changing with 
the load and the subsequent changes in the nature and 
extent of the yield of soil under the footing, they (the 
contact pressures) are not implicitly involved in the 
settlement calculations but are only worked out, 

(b) Rigorous Blast o -Plastic Analysis 

Most of the work in non-linear analyses is found 
in literature to have been based on the spatial separation 
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of elastic and yielded zones and treating them separately. 

In the rigorous analysis, the appropriate flow rules are 
applied to the soil in the regions which satisfy the 
appropriate yield criterian. In the rigorous analysis a 
need arises for dealing with the complexities of volume 
change, work done / extracted, work hardening or softening, 
large strains, tension cut-off etc. and for chosin,g the 
appropriate associated or non-associated flow rule, not only 
to get good predictions but also to get them at resonable 
computational cost. The latter is usually quite difficult 
especially for 3 -T> analysis and it is found fiom literature that 
simplified methods (pseudo-elastic, simplified elastic- 
plastic etc.) are sought for. A comprehensive discussion 
about rigorous elasto -plastic analysis is found in the works 
of Nayak and Zienkiewicz ( 1972 ), Zienkiewicz (1977) and Desai. 
and Christian (1977) . 

(c) Simplified Elasto -Plastic Analysis 

The stress-strain behaviour of soil can be represented 
in a very easy and simplified way as perfectly elastic- 
perfectly plastic (Pig, 2.2). In methods like PBM, it is not 
possible to deal with zero modulus, so either a small value 
is assigned or incremental stress-strain relations are 
introduced at the end of elastic limit (as in one of the 
rigorous elasto-plastic methods). But in the method using 
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numerical integrations, this can be achieved as described 
below in brief. 

Hain. (1975) suggested a method of 'limited contact 
stresses'. In this method the stresses only at the contact 
between the footing and the soil, which are the most signifi- 
cant ones from the points of view of local yielding and 
consequent settlement:, are checked against strength of soil 
and if they exceed the strength value, they are held at that 
value. Further compatibility of displacements is not 
considered in this region which therefore does not violate 
the assumed simplified plasticity rule (large strains 
without increase in stress). Hain's method is outlined 
below in brief. 

Initially a linear elastic analysis is performed. 

If the reaction at any point in the soil mass has exceeded 
the local ultimate bearing pressure it is held at that value 
and the stiffness matrix is modified effectively deleting 
the concerned nodes from the analysis. The excess load on 
the yielded elements is transferred to the non-yielded areas, 
causing redistribution. The iterative procedure is repeated 
until no further points yield during a cycle. 

In this type of analysis, the contact pressure 
distribution is obtained at each stage of loading along with 


settlements 
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A similar analysis is also performed by Poulos and 
Davis (1980) for piles. 

Hain has reported however that an analysis of this 
sort is not expected to predict correctly the deflection 
of the plate when yieldiiog has occurred over a large area 
of the plate. 

Study of centrally loaded annular footings using 
PEM with hyperbolic stress-strain law and Mohr-Coulomb 
yield criterian has been presented by Patankar (1932) for 
some particular sises of footings resting on the surface 
of sand . Por the cases considered it has been concluded 
that lesser total and differential settlements are obtained 
with increasing values of ratio of inner to outer radius of 
the footing for the s;rme intensity of loading. It has also 
been concluded that proximity of rigid layer improves the 
load settlement characteristics, 

(d) i)fo n-1 inear it y Through Stress-path Method 

In this method the stress-strain relations are 
governed by stress-path, For example a so3.ut.lon for circular 
footing on sand using stress -path method is given by 
Varadarajan and Arora (1982), 

(e) 2fon-linear or Higher Order Elastic Theories 

This approach has not been usually adopted as seen 
from the literature, probably because these theories involve 
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more niimber of soil parameters which are difficult and 
cumbersome to obtain, ir-hy ^ e l - a rst-i e f drtriula±±enr-giir^ 
h y Trus d^lX-4.-^-inel:ud:ed ^ - i n i --t h e--woric~of--^rlHgffl 

2.3 AmLITICAL WORK IK GOfflEGTION WITH PILED AMJLAR POOTIKGS 

There is no work found in the literature, directly 
dealing with piled annular footings. Poulos (1968) and 
Poulos and Davis (1980) have done elastic analysis of single 
and twin pile-rigid circular cap units, have extended it to 
pile-raft systems using interaction factors and have also 
suggested an approximate method to account for failure of 
piles. Such an extension from PEM results for a single unit, 
using non-linear interaction factor method is also suggested 
by Farasimhan (1982). 

The Finite element method or the Boundary element 
method for the three dimensional problem has not been 
employed in many cases probably because of very high 
computational cost involved, 

2.4 GOKSOLIDATION SETTLElffilKT OF GIRGULAR FOOTIKGS 
2,4,1 Flexible Gircular Footings 

Most of the available analytical solutions of the 
governing differential equations defining the phenomenon of 
time rate of settlement are for loading on homogeneous 
elastic half space. The solutions by Gibson and Lumb (1953), 
Davis and Poulos (1972) can be cited as those based on simple 
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diffusion theory. Solutions based on rigorous Biot's theory 
are given by De Jong (1957), Mandel (1957., 1961 ), Paria 
(1957), McNamee and G-ibson (i960) atnd Gibson et al, (1970). 
The diffusion theory, ( called as pseudo-consolidation theory) 
is reported to give results of settlement in between those 
from Biot's theory for v = 0 and v = o.5 (Davis and 

Poulos, 1972), 

2.4.2 Rigid Circular Bootings 

Analytical solutions for circular footings on 
semi -infinite soil are given by Agbe2ruge and Deresiewics 
(1975) and Ohiarella and Booker (1975). Agbezuge and 
Deresiewicz used Hankel-laplace transforms to give solutions 
for V =3 0 case. Ohiarella and Booker used integral 

transforms and gave an analytical solution , In both the 
cases it has been concluded that there is a decrease in 
contact pressure at the centre of the rigid footing, with 
time (accompanied by an increase in contact pressure near 
the edge of the footing). Yamaguchi and Murakami (1977) 
gave solution for strip and indicated an analytical method 
for the solution for a circular footing. 

2.4.3 Rafts of Finite Rigidity 

Oonsolidation problems for rafts of finite rigidity 
can be solved using finite difference or finite .element 
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methods, along the lines indicated by Sandhu and Yifilson 
( 1969 ) f Christian and Boehmer (1970), Hwang et al.(l971, 
1972 ), Yokoo et al, (1971), Booker (1973), G-haboussi and 
Wilson (1973), Yalliappan and lee (1975), Asproudas and 
Desai (1975) , Smith and Hobbs (1976), and Sandhu (1977). 

An interesting method, for consolidation under 
footings of buildings, is given by Komornik and Mazurik 
( 1977 ) in which a step-by-step iterative procedure with 
small time steps is adopted to find out and correct the 
compatible settlements. 

2.5 MISCBLLAilEOUS STUDIES 

A very limited literature is found to exist on 
laboratory and /or field tests on circular /annular footings. 
The references of Desai and Reese (1970) and Das and 
Gangopadhyay (1978) can be cited. There are several post- 
failure studies of failure of tanks/towers, e.g, Bozozuk 
( 1972 ) for an annular footing, Bell and Iwakiri (1980) 
for a tank. Failure of circular structures like oil storage 
tanks could be very costly. For example, it has been 
reported by Bell and Iwakiri (1980) that the failure on 
18th December 1974 of a hot oil tank 52.3 d. diameter at 
Mizurhima, Japan amounted to a loss of 150,000,000 dollars, 

A study of effects of period of construction and 
anisotropy with respect to permeability on the consolidation 
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settlement under a fully flexible circular footing is 
conducted by Madhav and Yitkar (1981). 

2.6 GO 'iCILUD I m REIMRKS 

It is observed from the literature that a solution 
for the entire load -settlement response of piled or unpiled 
annular footings with generalized results for direct use 
in the design is not available. The effect of finite 
rigidity of annul;.ar footings is a research area which also 
has remained unattended to* The consolidation under 
annular f ootings-rigid or flexible -taking construction 
time and anisotropy with respect to perraeability into 
account has not been studied. A quantitative investigation 
of the unattended aspects mentioned so far in connection 
with annular footings is essentially required, to obtain a 
better understanding of the behaviour and economical design 
of such footings which are going to be used as foundation 
for a large number of tower like structures to come up in 


future 
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CHAPTER 3 


HLASTO-PLASTIG SHTTLEMEHT aMLYSIS OP 
CENTRALLY LOADED RIGID ANMJLAa FOOTINGS 

3.1 INTRODUCTION 

For isolated structutes such as water tanks, 
cooling and other types of towers, silos, chimney stacks etc, 
footings in the form of an annular ring are well suited. 
Primarily such structures are centrally loaded and in 
several cases the footings can be classified as rigid 
relative to the soil on which they rest. A study of 
centrally loaded rigid ring footings (circular footing 
being a particular case) has been made in this Chapter 
for predicting the settlement response of the footings 
over the entire range of loading up to failure. 

Since plastic yield of soil, as a consequence 
of high stresses introduced under a rigid footing, would 
always exist even under a small load on the footing, the 
elastic analysis alone is not strictly valid. An attempt 
is made herein to introduce a simplified elasto -plastic 
analysis in which the concept of limiting the stresses to 
yield stress of soil has been introduced. Advantage is 
taken of the available elastic solutions covering the 
entire soil mass. For the reasons mentioned in Chapter 1, 
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the method of numerical intregration of the available 
solutions is used aind a technique of partial inversion 
of modified influence coefficient matrix has been adopted. 
The results from the present analysis, in the very low 
range of loading represent elastic behaviour to a 
sufficiently high accuracy. Such results when compared to 
available elastic solutions which are exact in some cases 
and approximate in some other cases, show a good agreement. 
In a few cases, experimental data and results from FEM 
are available and the results of the present analysis when 
compared to these show good agreement. 

Solutions for settlement analysis of rigid annular 
footings which cover an entire range of loading .and take 
elasto-plastic conditions of soil into consideration are 
not available presently. The conventional settlement 
analysis is still based on the application of theory of 
elasticity assuming a linear stress-strain response of soil. 

Relevent elastic solutions compiled by Poulos and 
Bavis ( 1974 ) consists of the works of Borowicka (1939) , 
Schiffman and Aggarwala (1961), Egorov and Serebryani (1963), 
Egorov ( 1965 ), Poulos (1968), Brown (1969) and Butterfield 
and Banerjee (1971). A. judicious extension of Mlovic's 
( 1973 ) solution for flexible annular loadings for predicting 
approximately the elastic settlement under rigid annular 
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footings is done by Mad-hav (1980). Egorov et al. (1979) 
have given correction factors, for depth of layer, to he 
applied to the elastic solution for semi-inf inite soil.- 
These correction factors correlated to field data appear 
to indirectly include the effect of elasto -plastic nature 
of the soil in field. A study of circular footings under 
simplified elasto-plastic conditions is contained in Hain's 
(1975) work. However Hain has reported that the solution 
is not valid when large portion of soil yields under the 
footing. A non-linear load settlement response under all 
loads, a highly non-linear and rapidly increasing settlement 
under higher loads, and a tendency of infinite settlement 
as failure is approached are not obtained through Hain's 
solution. The concept of limiting stress which Hain 
adopted has been used in the present work also, 

Eor simple soil conditions (homogeneous, isotropic) 
and the usual positions of foundation i.e. either on or 
inside a soil mass of considerable depth (treated as semi- 
infinite) or at the surface of a layer of finite depth* elastic 
solutions are readily available either in the form of 
expressions or tables (e.g. Poulos and Davis (1974)). These 
situations are taken up for study in this Chapter, Eor the 
central loading as considered in this Chapter, the footing 
discretization is limited only to the number of annular rings 
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and hence the size of the prohlem in terms of the size of 
the matrix is relatively extremely small, which is further 
reduced as yielding progresses. An overall study of a 
realistic load-settlement response with a significant improve- 
ment over elastic solutions has been achieved at a considerably 
low computational cost, 

3.2 ASSUI\aPTIONS Al® STATEMHllfD OP THE PROBLEM 
3.2.1 As sumpt ions 

(A) The general assumptions applicable throughout tlxLs 
thesis, except for the case of consolidation 
settlement (Chapter 7 ) are: 

(a) (i) The soil is homogeneous, isotropic and 

linearly elastic upto yielding. 

(ii) The footing has smooth contact with soil, 

(b) ■ Por the onset of yield and for the yielded soil, 

the following assumptions are made on the lines 
of Hain (1975), and Poulos and Davis (1980): 

(i) The soil yields plastically when the contact 
stress reaches a yield pressure and the 
stress is held at that value, 

(ii) Although at the points of contact, compatibility 
of foundation and soil displacements at a 
yielded element is no longer possible, dis- 
placements ansrwhere in the soil caused by 
the limiting stress q^. are sill given by 
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elastic theory. 

(B) For the purpose of the study in this Chapter, the 
following assumptions are also applicable: 

(i) The footing is in the form of an annular 

ring (a circle being a particular case) and 
is perfectly rigid. 

(ii) The footing is loaded vertically and centrally 
(giving rise to ?in axisyiDmotric situation) , 

3.2.2 Important Implications of the G-eneral Assumptions 

The general assumptions contained in 5.2.1(A) (a) 
and (b) above imply that the soil basically obeys linear 
elastic stress -strain relations and has a unique value of 
yield stress, for a given soil, irrespective of the geometry 
of the problem. The analysis based on these assumptions, 
although a significant improvement over purely linear 
elastic analysis, would amount to a simplified e last o -plastic 
analysis which is appropriate for clays, especially over- 
consolidated clays. 

The assumption regarding smooth contact between the 
footing and the soil would :amount to a slight overprediction 
of settlement, which would be less than 3 percent for rigid 
circular footings and even lesser for ring footings (Shi buya 
et al., 1975). The tangential forces due to adhesion are 
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more for lower values of Poisson's ratio and relative 
rigidity of the footing (Santos and Quera,198l). Still, 
since they are small and self-balancing, the effect on 
settlement is likely to be less in all the usual cases. 

Except for loading at the outer edge, the assumption 
of smooth contact leads to slight overprediction of bending 
moments in the footing, but is on the safer side. It can 
be said that this assumption is likely to give results, 
on the safer side in most of the cases and within about 
5 percent of those expected for an actual footing in the 
field, 

3,2.3 Statement of the Problem 

A smooth perfectly rigid annular footing of outer 
radius 'a' and inner radius *na’ carries a central vertical. 
load Q. The 'radius ratio' n is given by 
Inner radius of footing 

n = (5,1) 

Outer radius of footing 

n = 0 refers to the particular case of circular 
footing. 

The following types of problems are studied. 

Problem type A : Footing resting on the surface of a 
semi-infinite mass of soil (Pigs, 3.. 1 and 3.2), 
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Problem type B ; Footing is embedded to a depth 'c ' 
in a semi -infinite mass of soil (Pig. 3.2). 

Problem type 0: Footing resting on the surface of 
a soil layer of finite thickness ’h' underlain by a rough 
rigid base (Pig. 3.2). 

3,3 ANALYSIS 

3.3.1 Footing Discretization 

Most of the basic solutions obtained from the theory 
of elasticity, for the analysis of settlement, are for a 
point load. To deal with the variation of contact stress 
over the footing, the footing is discretized into small 
elements (Pig. 3.1) and it is assumed that the stress 
over a small element o is uniform and except when the 
effect of an element as itself is to be found out, it 
amounts to a point load acting at the centre of the element. 
By increasing the number of elements, the inaccuracy involved 
due to this assumption brought down into tolerable limits. 

3.3.2 Elastic Case 

Por expressing the relationship between elemental 
stresses and settlements, the basic form of Boussinesq's 
equation is conveniently adopted. The Boussinesq's equation 
expresses the vertical surface settlement at a point at 
a distance r from a point load Q acting vertically at the 
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surface of a semi -inf inite mass of soil as below(Poulos 
and Davis, 1974). 


Q ( 1- 
A E, r 


(3.2) 


where E and v are the elastic constants of the soil. 


Treating a uniform load q acting on a small area dA 
as point load, Eq, 3.2 can be written as 


q(l- V ) 


z 


E 


dA 

,^Tfr 


(3.3) 


In general, the settlement Pj_j at the centre of an 
element i due to a stress q^ acting on element j (Eig*3.l) 


is expressed on the lines of Eq, 3.3 as 


ID 


^D* 


a 


1-v‘ 

E 




(3.4) 


where is a dimensionless influence coefficient 
representing any available appropriate elastic solution and 
'a' , the outer radius of the footing is introduced for 
the purpose of non-dimensionalising. 

Bq, 3.4 in matrix form, representing settlements of 
all the elements, after non-dimensionalising, reduces to 

{p} = ( 1 -v^) LI] {q} (3.5) 

1 

where {p}= - { p} is the non-dimensional settlement vector, 

1 

tq} = rT {9> non-dimensional contact pressure 

jB . 
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a,nd [l] is the influence coefficient matrix. 


]?or a centrally loaded axi symmetric rigid footing, 
the settlement has to be uniformly equal to, say, p^. 

Eq. 3.5 can therefore be written as 

{p} = Pq {U} = (1- v^) [I] ^5} (3.6) 

where Pq = a *^o unit vector . 

For overall equilibrium of forces, it is necessary that 


Q = Z q^ dA^ (3.7) 

where, Q is the total applied load and dA^ is the elemental 
ll 

area of the j element of the footing over which contact 
stress q. is acting. 


Oombining and solving equations 3.4 through 3.7, the 
settlement of the footing, Pq is obtained as 


'0 


Q (1 

a 




E 


(3.8) 


Where 6 is a dimensionless factor given by 

1 

5 = (3.9) 

{ur[i]“'' {u} 

3.3.3 Elasto -Plastic Case 

Under any load, the contact pressures on some of the 
elements of the rigid footing would reach the value qy, 
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the yield stress for the soil under consideration, and 
it is assumed to be held at that value (Pig. 3.3). 

Let subscripts 'e * and 'y' represent the elastic 
and the yielded or plastic components respectively. 

Eq. 3.6 can now be written in the partitioned form as 

5e 

{ } ( 3 . 10 ) 

o ^ 

Py ^y 

Since compatibility of displacements is satisfied 
for elements in the elastic region only, Eq. 3.10 is 
written as 

2 ^ 

{ P } = { p. } = (1- V ) l_i •] {-_- -> (3.11) 

e ° 

where [l'] is the modified influence coefficient matrix 
corresponding to the coefficients of the elastic elements 
i.e. the first, say kg rows where kg is the number of 
elastic elements, For the yielded elements the stress 
value is known (=qy.). 

Ihe equilibrium equation 5.7 is now written in the 
form of sum of forces on elastic and yielded elements, 

Q == Qg + Qy = 2 qg dAg + S ly (iky (3.12) 

Going through the same procedure as in section 3.3.2 
above, the settlement of the footing once again can be 
expressed as 

III 


Pe 


{- 1“} = (1- v") [I] 
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Q 1 -V 


(3.13) 


^ a E 

where 6 is the dimensionless factor, which nnlike in 
elastic case, is now a function of not only influence 
matrix (Bq, 3.9), but also of the nature and extent of 
yielded region under the footing. 


Yet another convenient form of Eq. 3.13 is 

2 


1 - V 


a a 
^av 


E 


s 


where 


Q 


Q 


Area of footing 7ta^(l-n^) 
and Ig is the settlement influence factor. 


(3.14) 


(3.15) 


3.3.4 Elastic Solutions for Various Problem Types 


The influence coefficient matrix [l] is generated 
by calculating the influence coefficients representing 
the effect of contact pressure on element j on settlement 
of the centre of element i. The values of for various 
types of problems under consideration are obtained as follows. 
The elements i and j are both at the footing level in all 
the cases, 


(i) Problem type A : Footing on the surface of semi-infinite 
mass of soil. 

Under the assumed soil conditions, the Boussinesq's 
solution for a point load is applicable. The load acting on 
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an element 3 of area dA^ can be treated as a point load 
qj X dA^- acting at the centre of the area dA^, and for 
its effect on element i ( i 3 ), is obtained as 

dAj 

I = S_ (3.16) 

where dj_j is the distance between the centres of the 
elements i and 3 . 

For an element bounded by radial lines and concentric 
circles, of radii, r-i and r2 ( r 2 ^ i‘l)» "febe distance of 
centre of area, r , from the centre of the circle is 
obtained as 

ri+r2 (r2-r-^)^ 

2 6 (r2+r.^ ) 


-Sq, 3 . "IS is not applicable for the effect of an 
element on itself ( i = 3 ). The element is taken sufficiently 
small cwad the Bousinesq's solution integrated over a rectangle 
(Poulor and Davis, 1974) is used to get 


b VTl+m )+ m . 

I _ (iji ^ -I- jQ ) for i=3 ( 3 . 17 ) 


7" ( 1 +m ) ■ 


T(l+m )-1 


where b = width of element and mb = length of element. 


(ii) Problem type B : Embedded footing 


In this case, lilndlin’s (Poulos and Davis, 1974) 
solution is used in place of Boussinesq-solution, to get 
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dA4 3-4v 8(l-v^)-(3-4 V ) 

I { + — 

^ Sna (1 - V ) R2 

4(3-4^ )c^- 2 c^ 24 c"*- 

+ — + — } for i 7^ D ( 3 . 18 ) 

Ri Ri 

where R-^ = = distance between centres of elements i and j, 

c = depth of embedment of footing and 

R2 = T(4c^4-Rf). 


For i = j , Eq. 3.18 is not applicable. For this 
case, the Mindlin’s formula integrated over a rectangle 
(Poulos and Davis, 1974) is used to get 
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(3.19) 


where K = — ; — » 3- 4"^ , 

° 8n:(1-v ) 


K. 


5 - 12 V+ 8v , a =2c/l , p = 1/b, 


s = Y'(l+4a^), t = ■\r(1+P^('*+4a^)), 1 = length of longer 


side of rectangle and b = length of shorter side of rectangle, 



45 


(iii) Problem t37pe 0 : Footing on the surface of a soil 
layer of finite depth h. 

In this case the available elastic solution in 
the form of Tables (Paulos and Davis, 1974-) has been used. 
The settlement due to loading on element j is given as 

J 

2 % hE 

from which is obtained as 

Ip 

3 - for i j (3.20) 

^ 2 % ah 

where Ip is the value of the influence coefficient 
obtained from the Table by using lagrangian interpolation. 

For i = 3 and also for elements which are close to 
each other, the Ip values, extrapolated, become inaccurate. 
Therefore the sector curves, following the method given by 
Poulos and Davis (1974), are used, 

3.4 RESULTS AHD DISCUSSION 
3 . 4,1 Discretization 

The annular footing was divided, depending on the 
radius ratio n, into 8 to 16 rings and further into 10 to 20 
sectors. It was found that since the influence coefficients 
are more or less proportional to the inverse of the distance 
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between the elements, it was not necesseiry to subdivide the 
elements farther away from the neighbouring and the next-to- 
t he -neighbouring elements. The effect of subdividing a 
neighbouring element was found to be of the order of about 
2 percent only but it was found necessary to subdivide the 
neighbouring elements since the cummulatiye effect of 
neighbouring elements (which increase in number as distance 
increases) was found to be significant. The neighbouring 
elements and even the next-to-these elements were therefore 
suitably subdivided to obtain the components of the 
influence coefficient matrix. For finding the influence 
of a subelement on itself, formulae for rectangular loads 
or sector curves were used. It was observed that the best 
results are obtained by maintaining consistent order of 
accuracy (and not by increasing the accuracy wherever 
possible) as far as influence coefficients are concerned, 

3.4.2 Computational Economy 

In the present Chapter the elements involved at 
the stage of solution of matrix equations is extremely 
low in number. Further, in the general computer program 
developed in this thesis, a provision is included such that 
yielded elements are successively removed from' the matrix 
equations by substituting the known yield stress values. 
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Due to this and since soil below the footing is not 
discretized, the computation becomes highly economical 
and more so for the simplified e last o -plastic analysis 
used in the present work. For example, the entire set 
of results for annular footings on a semi -infinite mass 
were obtained with OPU time of the ord'er of 2 minutes 
40 seconds only on DEC. 10 system* The time for the other 
two cases was more because of involvement of more parameters 
(depth and Poisson's ratio), 

3.4.3 Footing Resting on the Surface of Semi-Infinite 
Soil Mass 

The complete load settlement response of a circular 
footing ( n= 0 ), as obtained by the present analysis is 
compared (Fig, 3,4) with the results given by Desai (1968), 
The agreement between the results of the present study 
with experimental results and those obtained by FEM 
is good, thus validating the use of the present method. 

Prom Fig, 3.4 it can be seen that all the three 
curves for n = 0 are non-linear throughout. This behaviour 
is attributable to the fact that even under light loading 
the soil near the edge of the footing yields, giving rise to 
the onset of plastic yieldii^ at a very low value of load 
itself. It is also seen that as the load approaches failure 
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load, the settlement increases rapidly with increasing 
load and tends to become infinite at failure load. Thus, 
in the present analysis, the entire load settlement response 
is realistically predicted. 

In the same Figure (Fig, 3.4) , load settlement 
curves obtained from the present analysis for typical 
values of radius ratio n ( n = 0.4, 0.6, 0.8, 0.9) are 
also shown for comparison between circular and annular 
footings. 

For direct comparison of behaviour, under the 
same total load, of a ring footing with that of a circular 
footing using settlement influence factor plotted against 
contact pressure, a correction factor to Ig is introduced 
as below, 

let ’b ' be the outer radius of a ring footing which, 
with a given n value, has the same, area as a circular footing 
of radius 'a* (Fig. 3.5). Then, for equal area, 

% 8^ = rc b^ (1-n^) (3,21 ) 

2 2 2 

Also can be written as Q/ita and Q/nb (1-n ) 
for circle and ring respectively, where Q is the (same) 
total load. Then using Eq, 3. 14, settlement of circular 
footing under load Q is given by 
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and, settlement of ring footing under the same load Q 
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where and Ig^, are for the same average contact pressure, 
so that = lo* 

SC ST fc> 

Denoting the correction factor by Rgg^ and using 
Sq. 3.21 above, 


•^or b 1 

^oc a Y‘(l_n^) 


(3.24-) 


For comparison of the response of a ring footing 
v/ith that of a circular footing under the same total load, 
the Ig values obtained from the charts should be multiplied 

The values of work out to about 1, 1.02, 

1.08, 1.25, 1.67 and 2,3 for n= 0, 0.2, 0.4, 0.6, 0.8 and 
0.9 respectively. Applying Rgg^ to n = 0,8, the settlement 
of equivalent annular footing in shown in Fig. 3.4. It is 
seen that the settlement for the equivalent annular footing 
is also less than that of circular footing. This is true 
for other values of n also. The settlement at all load 
levels for annular footings is less than that for a circular 
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footing and more so for higher values of n. 

Pig. 3.6 shows the values of dimensionless 
settlement influence factor Ig obtained from the present 
analysts, Por the purely elastic case, the dimensionless 
factor 6 is a function of the influence coefficient matrix 
(hq* 3.9), which depends on the geometry of the problem 
only. However, for elasto -plastic case, it is found that, 
since settlement depends also on the nature and extent of 
yielded region under the footing, 6 or Ig is a function of 
radius ratio n as well as of load level indicated by the 
ratio of average applied load to yield stress of soil. 

6 or Ig = f(n, lavAy) (3.25) 

Because of this, curves for Ig in Pig, 3.6 have 
been drawn against q ^/q values for various values of n. 
Pig, 3.6 which is a design chart for settlement, is to be 
used in conjunction with Bq. 3.14 namely, 

1 - 

JCij 

Prom Pig, 3.4 and Pig.. 3.5 it can be seen that the 
settlement curves are non-linear tliroughout , but they are 
less non-linear at lower load levels (lower q^^/oiy values) 
and also over a larger range of load level for higher values 


of n. 
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On studying the settlements in the range of 
ratio ( a form of Factor of Safety) of about 5 to 1.2, it 
can be seen that for n = 0.4, 0.6 and 0.8, the settlements, 
after applying the correction Req’ roughly 80 percent, 

75 percent and 50 percent respectively of settlement for 
circular footing (n=0) for the same total load (and the 
same area of footing). For n = 0.9, this percentage varies 
from about 50 percent to 37 percent as varies from 

5 to 1.2, This clearly shows the advantage of ring footings 
as far as settlement response is concerned. It must be 
mentioned that for narrow rings, the settlement, ixicreases 
almost linearly at a slower pace until qy/^y ratio of 
about 1.1 is reached. Thereafter the settlement is likeljr 
to increase rapidly for q^/q^^ values less than 1.1, If at 
all a situation of loading very near to failure is likely, 
a circular footing will show up excessive settlements well 
before factor of safety of unity is reached, whereas a 
narrow ring footing will collapse suddenly winen the factor 
of safety falls below about 1,1. 


The design chart of Fig. 3.6 can be used in a number 
of ways. For example, for a known value of load and 
dimensions of footing, the ratios q^^^/q^ and n can be worked 
out and the corresponding settlement influence factor Ig 
can be directly read and used in Eq, 3.1 4, to predict settlement, 
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The permissible loading or even the required radius ratio 
can be worked out if other relevent data is available. 

During the course of solution by the present 
analysis, the contact pressures and the progress of yield 
zone also are obtained at various load levels, E‘ig,3.7 
shows typical contact pressure distributions. For the 
sake of comparison they are drawn for various values of 
Q-y/^av* case of a circular footing, the yield starts 
at the outer edge and proceeds inwards under increasing 
loads. In case of an annular footing, the yield starts 
at the outer as well as the inner edges but the progress 
is faster from the outer edge (Fig, 3.7). 

The portion of the load settlement response obtained 
by the present analysis without considering yield represents 
elastic behaviour. The results in this range are compared 
with available solutions from theory of elasticity. The 
settlements and contact pressures obtained by Egorov (1965) 
are compared with the present analysis (Table 3.1 and Fig, 3.8) 
and good agreement is found to exist. 

From Table 3,1, it can be inferred that at a 
sacrifice of 5 percent increase in settlement, a circular 
central portion of radius about two-thirds of the outer 
radius can be simply dispensed with amounting to 45 percent 
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saving in footing material. 

lABLS 3.1 : TAIUSS OF w(n) FOR RIGID RING (ELASTIC CASE) 


n 




.8 .9 


w(n) 

■Egorov 

(1965) 


w(n) 

jPresent 

Work 

(Elastic 

Portion) 


.50 .50 .51 .52 .57 .60 



l^lll^^w(n) 

E a 


.50 .50 .51 .52 .55 .61 


This inference drawn on the basis of elastic theory 
alone could be misleading since yield conditions will always 

exist and they would increase settlement. Thus, under the 

2 

same total load, will increase in proportion to (l-n ) 
and w(n) obtained from e last o -plastic analysis as done in 
the present work will change since it will be a function of 
qg_^/qy,. For example, it has been worked out that, against 
45 percent saving as predicted by elastic theory, only 
16 percent and 11 percent saving is possible for 
ratios of 0.5 and 0.67 respectively. It has to be however 
mentioned that due to plastic yielding of soil, redistribution 
of contact pressure takes place which would reduce the 
maximum bending moment in the footing slab and would lead 
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to a saving in reinforcement if the loading is uniformly 
distributed over the footing, 

5.4.4 Embedded Footing 

The results are obtained for this case exactly 
on the same lines as for the surface footing except that 
the additional parameters 'c ' (depth of embedment) and 
' V * (Poisson's ratio) come into picture. The embedment 
herein only means load transfer from footing to soil at 
depth c. However the load refers to net load at footing 
level and also refers to the yield value at the depth c. 

It has been found from the solutions obtained from 
the present study that the nature of the contact pressure 
distribution beneath an embedded footing is not noticeably 
different (Fig, 3.7) from that beneath a surface footing. 

It has been also found that the nature of the load-settlement 
response is also very similar for the embedded and surface 
footings (Fig, 5.9(a)). Hence it has been possible to 
express the settlement of embedded footing in terms of 
that of the surface footing with only the application of 
correction factors to account for the effects of depth of 
embedment c, the radius ratio n and the Poisson's ratio, v . 


% 


The settlement of embedded footing is expressed as 
1 - 

I _ “ -r. r, n n (3.26) 


^av 


B 


°V 
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where Ig is the settlement influence factor for surface 
footing corresponding to at depth c and G^, 0^ and 0^ 
are correction factors for c, n and v respectively. The 
values of these correction factors are presented in 
Tables 3,2 and 3.3. 

Ig is to be read from the design chart of Pig, 3.6 
given for surface footing only, but against qav/% 
in which should correspond to average net loading at 
the depth c and q^. should be the yield value at depth c 
which would be higher than that for the surface footing, 

TABLE 3.2 : SETTLEIffiKT GORRECTION PAGTOR 0 ^ FOR DEPTH OP 

EMBEDMENT (Eq. 3.26) 


c/a 0 .05 .1 .25 .5 .75 1 


G 


c 


2 3 4 6 8 10 15 or 

more 


1 .99 .98 .95 .91 .87 .83 .71 .65 .62 .58 .56 .55 .53 


TABLE 3.3 : SBTTEEIiElfD OORPv.liaTIOH EAOTORS and Gy EOR RADIUS RATIO (n) AM) 
POISSOF'S RATIO (v ) (Bq. 3.26) 
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The available relations for ultimate bearing capacity 
at depth = 0 and at depth = c can be used to get the 
ratio between these two and finally between q^^ and q at depth 
c. One such relation, for example, based on empirical 



Since settlement of an embedded footing is related 
to that of surface footing (Eq. 3,26), most of the discussion 
given in Section 3.4.3 above holds good for this case also. 

The load settlement response is non-linear throughout, large 
settlements are predicted near total failure and ring footings 
are advantageous over circular ones. 


The effect of embedment is particularly signif leant., 
firstly, as can be seen directly from Tables 3.2 and 3.3. 
Secondly, for higher values of depth of embedment, qy. is 
higher, is lower and hence I 3 , read against qav/'ly 

from Fig, 3.6 is lower, which amounts to lower settlement 
obtained from the expression given in Bq, 3.26. 


The effect of radius ratio as apparent from Table 3.3 
is small, but it is actually quite significant due to its 
inclusion in the settlement influence coefficient itself. 
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The term Ig. 0^ 0^^ 0^ in the expression for settlement 
(Eq. 3.26) is a function of not only c, n and v j but also of 
'lyVly . 


^s * ^n^^ > ly / ^y ^ (3.28) 

O 0 

Thus, the design chart of Eig, 3.6 for surface 
footings can be used in conjunction with Tables 3.2 and 
3.3 and with an appropriate ratio of q /q„ , for 

^ Q ^ 0 

prediction of settlement of embedded rigid ring footings. 


The nature of contact pressure distribution, the 
changes in it at various load levels and the nature and 
extent of yielded zone under the footing are so similar 
to those for surface footing that the difference is not 
noticeable (Eig, 3.7). 

It has been observed that (e.g, Eig. 3.9a) the 
settlement curves obtained by applying corrections to the 
c/a = 0 curves agree with the directly obtained curves 
within about 4 percent at all load levels. 

As discussed in Section 3.4.3 for surface footings, 
a saving is possible by dispensing with some central portion 
of a circular footing and the saving is more for higher values 
of l-^r/^-av* lu case of embedded footing, since q is higher, 

J V 

the ratio qy/q^-y the same load and footing will be higher 
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in comparison to surface footing. Therefore the saving 
in footing material will be more for embedded footing even 
when partial yield conditions exist. 

Pig. 5.9(h) shows the comparison of results from 
the present analysis in the low range of representing 

elastic behaviour, with those of Butterfield and Baner jee 
( 1971 ), for embedded circular footing. It is seen that there 
is a good agreement. 

3 . 4.5 Surface Footing on Finite Layer of Soil 

For this case, the available elastic solution in 
tabular form (Poulos and Davis, 1974) has been used. Unlike 
embedded footing, the load settlement response in this case 
is found to be dissimilar in nature to that of surface footing. 
The settlement is expressed, as 

»- 29 ) 

where is the settlement influence factor, the subscript h 
denoting the case of soil of finite depth h. 

In this case, Ij;^ is a function of depth of layer, h, 
and Poisson's ratio also. 

Ih = f(n, lav/^ly’ ^ ) (3.30) 

Design charts are presented in Fig. 3.10 to readily 
obtain the value of for the required data. 
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From Fig. 3.10, it can be seen that as in the 
earlier two cases, the settlement response is nonlinear 
throughout, there are large settlements near total failure 
and annular footings are advantageous over circular footing. 

The depth of compressible layer has a significant 
effect on settlement. The settlement influence coefficients 
from the present work, as tends to zero, correspond 

to elastic behaviour and agree fairly well with those given 
for circular footing by JPoulos (1968) being 2 to 3 percent 
on higher side and agree within 1 to 2 percent with those 
for annular footings given by Madhav (1980). However the 
effect of depth of layer under all values of (lay/ly is of 
particular interest in the present study. These effects 
can be directly taken into account if the design charts 
in Fig, 3*10 are used. Hot enough data is available for 
comparison other than the results in the elastic range. 
Egorov et al. (1979) have suggested some correction factors, 
to account for finite depth, to be applied in the expression 
for elastic settlement of rigid annular footing resting on 
semi -infinite soil. These factors are reported to yield 
results close to those obtained from field data. Since 
field conditions would always involve plastic yield, these 
factors can be looked upon to include the effect of plastic 
yield or load level indirectly. For comparing their values 
with other solutions a factor G- is introduced as below which 
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reflects the effect of depth of layer. 



(3.25) 


where p is the settlement for footing resting on soil 
of finite depth and Pe is the elastic settlement for 
semi -infinite soil under the same footing and loading. 
Some of the values of 0 which are taken directly from 
the Tables given by Egorov et al, (1979) are quoted in 
Table 3.4 along with worked out values based on elastic 
solutions by Poulos (1968) and those obtained from the 
present work using elasto -plastic analysis* 


TABLE 3.4 : VALUES OE C (Eq. 3.25) FOR TIPIGAL OASES 
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The corrections given by Egorov et al, are lower than 
in the elastic case. 0 worked out from the present work 
depends on Poisson's ratio and ratio and the value 

in any case is lower than that in elastic case. The lower 
values of G indicating higher settlements; can be attributed 
to elasto-plastic settlement which occurs in field or which 
is predicted from the present analysis. Since it is more 
convenient to use directly the charts of Pig. 3.10, factor G 
has not been used to e.ypress the effect of depth of layer 
except for the comparison discussed above. 

Prom the design charts of Pig. 3.10, it can be 
seen that the effect of h/a is more pronounced for lower 
values of n, for lower values of qav/% higher values 
of Poisson's ratio. For example, considering h/a = 0 and 4, 
the values of Ig/Tj^^ (indicating ratio of settlement for 
semi-infinite soil to that for finite depth) for Poisson’s 
ratio = ,5 and q^^/q^ = ,3 are about 1,38 and 1,25 for n = 0 
and n=,8 respectively; the value of Ig/^ ^'or n=0, Poisson's 
ratio =,5 and qg_y/q.y.= ,6 is 1.34 and for n=0, q^^^/qy. =.3 and 
Poisson's ratio =0, it is 1.21, For higher values of n, 
the effect of depth of layer although less pronounced, 
amounts to overall less settlements. Hence an annular 
footing resting on a finite depth of soil will show very 
small settlements. 
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From the design charts (Fig. 3*10) it can also be 
observed that as far as overall load- settlement response 
is concerned, the effect of having less depth of layer is 
similar to the effect due to increasing n. In both the 
cases, the settlements are less, the load-settlement 
curves are nearly linear except in the higher ranges of 
load level and a somewhat sudden failure occurs as load 
approaches failure load. This behaviour for lower values 
of h/a is explained as follows. 

The contact pressures tend to be like unifom as 
h/a decreases. Initially, under increasing load, the yield 
zone progresses slowly. At the same time the contact 
pressures in the elastic zone increase and become more 
uniform, A stage comes when the contact pressures in the 
elastic region becomes almost uniform and only slightly 
less than the yield stress, A small increase in load at 
this stage is enough to cause a rapid advance of yield zone 
and a consequent high settlement. 

Fig. 3. 'll shows the progress of yield zone from 
some topical cases. 

The value of q for footings resting on finite depth 
of soil may be taken as the same for footing on semi -inf inite 
soil for h/a > 1.42. For lower values of h/a, qy. appropriate 
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for the required value of h/a is to be used. An approximate 
correction to qy for depth of layer can be obtained by 
using Mandel and Salecon's (1969) values of . These 
values are for strip footings on clay, but their mutual 
proportions may be assumed to be same for annular footings, 
^ith this assumption, the correction factors, based on 
Mandel and Salecon's values are shown in Table 3.5. 


TABLE 3.5 : APPROXIliA.TE OORREOTIOIB TO q,, EOR DEPTH OP 

y 

LAIER OP SOIL 


h/a 

>1.42 1 

.667 

.5 

.4 

.335 

.25 

.2 

Correction 

Pactor 

1 1.02 

1.11 

1.21 

1.30 

1.40 

1.59 

1.78 

The 

only change 

, due to the 

correction for 

depth 



mentioned above, in the use of design charts (Pig. 3.10) 
will be because of the corrected value of q^.. The revised 
value of qg_^/qy is to be used. Due to this correction, 
the settlement at any load will work out less in comparison 
to that obtained without correction except for h/a > 1.4^. 
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3.5 C0M3LUSI0I© 

^ simple elasto -plastic analysis using available 
elastic solutions and a modified procedure has been 
developed to obtain the entire load-settlement response 
of centrally loaded smooth rigid ring footings resting on 
or embedded in a semi-infinite soil mass or resting on a 
finite depth of compressible layer of soil underlain by a 
rough rigid base. It is seen that the results when compared 
with available solutions show good agreement with experimental 
as well as PEM results. Oomputational economy is achieved 
mainly due to the absence of soil mass discretization and 
due to condensation of flexibility matrix as plastic yielding, 
progresses. 

The following are the conclusions, drawn from the 
present study, 

(1) Due to the absence of purely elastic conditions at 
any level of loading, the load-settlement response 
is nonlinear throughout. 

(2) The settlement is implicitly related to the level 
of loading with respect to the ultimate load, 

(5) For the same load annular footings show lesser 
settlement in comparison to a circular footing, 

(4) The plastic yield starts near the edge/edges and 
proceeds towards centre under increasing load, 



66 


causing redistribution of contact pressure. Under 
increasing loads, the contact pressure tends to 
become uniform and this is likely to limit the 
induced bending moments in the footing loaded 
uniformly. 

(5) The response of an embedded footing in a semi- 
infinite soil is qualitatively very similar to 
that of surface footing resting on semi -infinite 
soil. It has been therefore possible to provide 
simple correction factors, 

( 6 ) The load-settlement response tends to become less 
non-liijear for narrower rings and/or shallower depth 
of soil layer. Also in these cases the settlement 
increases with load level at a slow rate but it 
increases suddently near to the ultimate load. This 
needs attention from a practical point of view, 

A design procedure has been outlined for predicting 

settlements in all the three cases considered. 
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FIG*3*4 TYPICAL LOAD SETTLEMENT CURVES 



WIDTH RATIO w 

FIG-3-5 EQUIVALENT CIRCIII AR ANH AMMlil AP 





u 


ns 0 



.3-7 TYPICAL CONTACT PRESSURES BENEATH RIGID FOOTING 
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FIG‘3*8 CONTACT PRESSURE BENEATH 

RIGID RING ON SEMHNFINITE SOIL 
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FI6-3-9Q EMBEDDED FOOTING-COMPARISON 
OF LOAD SETTLEMENT RESPONSE 
WITH SURFACE FOOTING 



FIG3«b EMBEDDED FOOTING -COMPARISON 
OF ELASTIC SETTLEMENTS OF 
CIRCULAR FOOTING 







Sjv'qy 


2 



M]*J0V3D> 


o 


FIG. 3-10 FINITE DEPTH SETTLEMENT FACTORS 
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FI0 3-10 FINITE DEPTH SETTLEMENT FACTORS 
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FIG 3-10 FINITE DEPTH SETTLEMENT FACTORS 
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FIG.3H TYPICAL CONTACT PRESSURES BENEATH 
RING FOOTING ON FINITE LAYER 



CHAPTER 4 


ELAST0-P31ASTIG SETTLEMEITT AMLTSIS OP GEMIRAXLY LOADED 
' AHMJLAR FOOTINGS OP FINITE RIGIDITY 

4.1 INTRODUOTION 

The footing slabs used as foundations for heavy 
structures have to be strong enough and therefore can be 
classified as rigid. But rigidity of a footing relative 
to the rigidity of soil is an important consideration for 
the analysis of settlements-mauciraum and diff erential-and 
of bending moments and shear forces in the footing. Also 
in certain cases like storage tanks at or below ground 
level, the bottom slab or plate resting directly on the 
soil can be thin. It is therefore necessary to study the 
effect of relative rigidity of a footing. With increasing 
flexibility of the footing, the contact pressure distribution 
tends to be uniform and this aspect is important from the 
point of view of the elasto-plastic anaQ.ysis and several 
design aspects. 

Since the study of finite flexibility of footing 
involves interaction between soil settlements and plate 
displacements, a study of bending moment and shear is also 
included along with the elasto-plastic analysis of 
settlement. The nature and positioning of the imposed 



86 


load on the footing has a significant effect and is also 
made a part of the study, Another specific feature of 
the study of footings with finite rigidity is the 
possibility of loss of contact of the tooting from the 
soil under certain values of relative rigidity and 
certain loading configurations. 

One of the earliest important studies of flexural 
rigidity in connection with smooth circular footing is 
found in the work of Borowicka (1939). In this study, 
the load is in the form of a central concentrated load 
and the footing rests on semi-infinite soil. For uniformly 
loaded circular footing, the solution is given by Brown 
(1969) for semi -infinite soil as well as for finite depth 
of soil, The bending moments are studied for rigid circular 
plates on finite layer by Egorov and Serebryanyi (1963) 
and Poulos (1968), All the solutions mentioned so far 
are based on the assumption of purely elastic behaviour 
of soil and these can be referred to in the compilation 
by Poulos and Davis (1974). Some of the studies based on 
elastic theory but considering loss of contact are by 
Oheung and Hag (1968), Svec (1974) r G-ladwell and Iyer 
( 1974 ) Yalliapan et al. (1975), Scott (1978) . 

The study of smooth annular footings of finite 
rigidity from the point of view of settlement, bending 
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moments and shear forces in the footing taking elasto- 
plastic behavioTor of the soil into account is not found 
in the literature. 

The use of thin plate theory to deal with the 
footing-soil interaction is in vogue and is presumed 
to be acceptable if the aim is an overall study of 
interaction behaviour and major structural design aspects 
of the footing and not the details of inplane stresses 
etc,, in the footing. 

The coupling of the plate-bending problem with 
that of soil deformations is achieved by inverting one 
of the plate rigidity or soil flexibility matrices and 
just adding if only elastic behaviour is involved. 

Since in the present work, elasto-plastic behaviour is 
considered, a somewhat complicated procedure had to 
be adopted, but since the size of the problem in\ terms 
of the size of the matrix to be solved is quite small, 
the computations are still within highly economic limits. 

4.2 ASSUMPTIONS AND STATEMENT OP THE PROBLEM 
4.2,1 Assumptions 

The assumptions contained in Section 3.2,1 (A) (a) and 


(b) are applicable. 
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For the purpose of this Chapter, the following 
assumptions are also applicable. 

(i) The footing is in the form of an annular ring 
(a circle being a particular case), 

(ii) The footing is loaded vertically and axisymmetrically 
and has finite rigidity, 

(iii) The footing is treated as a thin plate* 

4. 2.2 Statement of the Problem 

The statements contained in Section 3.2,2 and 
Figs. 3.1 and 3.2 are valid except that in this case 
the footing has finite rigidity and that the vertical 
loading may be positioned in any axisymmetric configuration. 

Since several combinations and parameters are 
involved, only some cases of load distributions have been 
studied from which at least a qualitative conclusion 
of behaviour under any configuration of loading could be 
drawn. The following load configurations are studied 
(Fig. 4.1). 

(l) Uniformly distributed load (U.U.l,) over the entire 
area, 

U.B.Ii. over the outermost one-ninth width of the 
ring (or circle). 


( 2 ) 



(3) U.B.L, over the innermost one-nineth width of 
the ring (or over a circle of one-nineth the 
outer radius in case of a circle), 

( 4 ) U.D.L. over the central one-nineth width of 
ring or circle. 


4.3 AMLYSIS 

4.3.'* The Soil Flexibility Matrix 


As for rigid footing (Chapter 3)j the footing is 
discretized into a number of small elements and subelements 
The settlement, of an element i due to contact pressure 
p^ on element 3 is expressed by 

1 - 

s 

p = p . a - — - — I . 

^3 ^3 Ec 


's 


ID 


(4.1) 


Here subscript s has been used to indicate 
parameters pertaining to soil. 

The above Sq. 4 .I in matrix form, for all the 
elements of the footing, condensed to the size m x m , where 
m is the number of radial elements, can be written as 

{p} = (1-Vg) [I] {p} (4.2) 

1 

where {p} = -- {p> is the non-dimensional settlement 

a 

1 

vector, {p} = “ {p} is the non-dimensional contact 

E 

s 

pressure vector and [l] is the influence coefficient matrix 



90 


The condensation of [l] matrix into raxm size is possible 
because of axisjnnmetry. 

For convenience, Eq, 4,2 is cast into the following 

form 

{p> = [zl (p) (4.3) 

where [z \ represents non-dimensional soil flexibility 
matrix of size mxm. 

4.3.2 Footing Rigidity Matrix 

For an axis 3 mimetric ring shaped thin plate, the 
governing plate equations are (Timoshenlco and Krieger, 

1959) : 


B 

d'^w 

2 d^vy 

1 

,2 
d w 

1 

dw 


( + 
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dw 


P 
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o 

— ) 
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dr3 

r 

dr^ 

r^ 

dr 


where B = Bpt^/I2(1- ’^p ) is the flexural rigidity of 
the plate (footing), t is the thickness of the plate, Ep 
is the Young's Modulus of the plate, Vp is the Poisson's 
ratio of the plate , w is the deflection of an element 
situated at a radial distance r from the centre, q is the 
intensity of load on the element, p is the contact pressure 
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On the element, M and F are the bending moment and shear 
force per unit length respectively. 


Using Taylor's expansions and central difference 
scheme, Eq, 4.4 is converted into finite difference form, 
for node i (Fig, 4.2) as 

^ 2 ^i 

7—4. O-Yi) + (-4. + 2y^ _ - -J-) 

(8ri) 

+Wj^(6 + 2Yi)+ Wi+1 (-4-2yi-y| + ~) 

+ Wj_^(l+y^)) = - Pi (4.9) 

where y. = 6r/ri , Ti being radial distance of r element, 
Eq, 4,9 can be written, in matrix form, as 


D fj» 

_ L^i i [w._2 w^_4 1 


{ q - P } 


( 4 . 10 ) 


where ai 


^ " ^i 3 

2 n 

~ 4 + 2yi - yi - ^ 


d. 

1 


6 + Zyl 

2 

-4 - 2y^ - y^ t ~ 


and 


1 + yi . 



For convenience, the m radial elements are 

denoted by numbers 3, 4 m+2, 1, 2, m+3 and m+4 are 

the fictitious nodes added (Fig. 4.2). Bq. 4.10 embodies m 
equations for m nodes but would include 4 unknown deflections 
at the fictitious nodes 1, 2, m+3 and m+4. lo eliminate 
therse, the following (free edge )boundary condition for the 
footing (plate) are used, i.e. 


1) Moments at the inner and outer edges of ring 
footing are zero. 

2) Shear forces at the inner and outer edges of 
ring footing are zero. 


Considering the inner edge, (Pig. 4.2 one 
the boundary conditions as 


M, 


'A 


%- 1 = - B ( 

1 d 

1 = B ( - 


^2 
d w 

4* 


dw 

dr2 


r 

dr 


1 

d w 

1 

-- 


I IWI 


r3 

r 

dr2 



) = 0 
dw 


) = 0 


can. write) 


(4.11) 

( 4 . 12 ) 


Using central difference, these equations are converted into 
the following finite difference matrix forms. 

1(1/2) (-1/2- (-1/2+ ''pyaHl/2) \ ’'i-1 "l "1+1 1* 

= [O 0 0 0 l’' (4.13) 

|(-i^ - )(3- ^ )(-3- ^ -ya ) r ) I • 

^i-1 "^i "^i+ll^ =[0 0 0 of 


(4.14) 



where = Sr/rg^, Xg^ being radial distance of A from 
centre. 

Similar expressions are written for the outer edge. 

]^or a circular footing the above boundary conditions 
for the inner edge are replaced by the following: 

1 ) Because of symmetry about the centre, the deflection 

at the fictitious node 2 is mirror image of deflection 
at 3 i.e. W2 = w^. 

2 ) Similarly w-| = w^. 

Thus using the four equations obtained from boundary 
conditions, w-j , W2> w^^^ and are eliminated from 

Eq. 4.10 leaving m equations to be solved. 

For convenience, Eq, 4*10 after incorporating the 
boundary conditidns, is cast into the following form 

[?] {w} ={q-p} = {q} - {p} (4.15) 

- 1 

where {w} = — {w} represents non-dimensionaj. deflections 

^ 1 

vector, { q - p } = -5-- { Q. - P } is the nondime nsional net 

•“S 

pressure vector and [p] represents nondimensional plate 
rigidity matrix of size mxm . 

4.3.3 Coupling of Soil and Plate Responses :Elastie Oase 
Inverting Eq, '4,5 
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For compatibility of vertical displacements of 
plate and soil, the deflections of plate should be equal to 
settlements of the soil, hence 

{w} = {p} (4.17) 

Using Eqs, 4.15 through 4.17, 

[Z + P] {P } = {q> (4.18) 

Phis equation is solved to get settlements. Contact pressures 
are obtained from Eq, 4.16, 

4.3,4 Coupling of Soil and Plate :Elasto -Plastic Case 

For the elements where the contact pressure exceeds 
yield stress q^. of the soil, say k in number, the displacement 
compatibility is not to be satisfied and at these elements the 
contact pressure is held at qy. It is necessary to remove 
the equations corresponding to yielded elements from Eqs. 

4.15 and 4.16 before coupling is done. Considering Eq. 4.16, 
the rows corresponding to yielded elements, are eliminated, 
one by one, and in doing so, using the known valqe q , the 
[z]matrix as well as the right hand side of the equation 
get modified. let {p-| } indicate the modification to the 
r.h.s. vector. After removing the k rows corresponding 
to yielded elements, 

[z'-”'] {p > = {p } _ {p^ } (4.19) 
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where is now a matrix of size (m-k) x (m-k), 

The effect of contact pressures acting on yielded elements 
is now contained in the vector {p^ }. 

Similarly, Eq. 4.15 is treated to obtain 

[P'j {w} = {q} - {p} - {P 2 } (4.20) 

where {p2^ indicates the effect of yielded elements and 
[P'}, is the modified form of matrix, now reduced to 
the size (m-k) x (m-k). 

The settlements { p} in Eq, 4.19 and the plate 
deflections {w} in Eq. 4.20 can now be equated since they 
belong to elements of elastic zone only. Therefore, from 

Eqs, 4.19 and 4.20 one gets, 

1 

[z‘" +p']{p} = {q} - {p-,} - {Pq}- (4.21) 

Eq, 4,21 written for compatibility in elastic region only 
(see Chapter 3) is solved to get the settlements (=plate 
deflections) in the elastic region. The contact pressures 
in elastic region are obtained by using Eq, 4.19. The contact 
pressures in the yielded region are all equal to qy. The 
plate deflections in the region where soil has yielded are 
obtained by solving Eq, 4.15 for {w> , 

4.3.5 Separation or Tension Cut-off 

In footings with finite rigidity, situations can 
arise, such as when the load is concentrated at the centre of 
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a circular footing, creating separation of the footing 
from the soil in some region. This occurs since soil can 
not take tension. To account for this, first one can find, 
from elastic solution, elements where contact pressures 
become.- negative. At these elements say k' in number, 
the contact pressure value is held at zero. The compati- 
bility of displacements at these elements is not satisfied 
and hence the procedure follows exactly on the lines of 
yielded elements. For the case of elasto-plastic analysis 
with separation, Sqs. 4.19 and 4.20 are obtained by excluding 
the k' elements so that the size of the matrices obtained 
becomes (m-k-k') x (m-k-k*). The rest of the procedure 
remains the same. 

4.3.6 Bending Moments and Shear in the Footing 

While plate deflections are obtained directly from 
the solution of Eqs.4.18 or 4.21, the deflections at the 
fictitious points are obtained from the appropriate 
expressions used for boundary conditions as given in 
Section 4.5.2 above. 

Finally Bqs. 4.jU^and 4.42' in Finite Difference form 
as given below are used to get bending moments and shear 
forces at any node (centre of element) i of the ring footing. 
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M. 




■ ^2 wj^_i-2wi+(i+^ yiK+n^ 

:^3(- j ’'i-2+(i+yi+ 5 r^Wi-i-syi 

+(-l+yi- I y?) j Wi+g) 


(4.22) 




(4.23) 


where - 6r/r^ , being radial distance of node i 
from centre. 

For a circular footing, the bending moment at the 

centre is the only important quantity, for which instead of 

Bq. 4.22 the following equation is used 

I)(1 + Vp) 


Mo 


(w^-Wj) 


(6r)2 

The bending moments referred to above so far are 
radial bending moments. The tangential bending moments 
are easily obtained from the following relations 

1 dw 


(4*24) 


-H ( V 


,2 

a w 


P dr2 


r dr 


) 


(4.25) 


Eq. 4.25 in Finite Difference form for node i of a ring 
footing will be 


D 


M 


T - “ 
■^1 


(('' + ^ )’»i_1-2''p w^+( V - (4.26) 


(6r) 

where yj_ = 6r/rj|_ as before 

At the centre of a circular footing, the radial and 
tangential moments are equal. . 
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MTq = (4.27) 

The values of radial or tangential bending moments 
will be maximum at certain points and further, the higher of 
the two will be of particular interest. The absolute 
maximum values of bending moment and shear force have 
been denoted by Mq . and respectively, 

4.4 RESULTS AID DISCUSSION 

4.4.1 General 

The parameters used in the present study are the 
Poisson’s ratios for soil and footing, ( Vg and v^) and the 
relative rigidity K defined as 

K = C1- \ (+.28) 

£lg 

All the other required elastic derived quantities 
are expressed in terms of '’s’ "^p 

For brevity, and are assigned the values 0.5 
and 0.25 respectively in major portion of the work included 
herein. 

During the numerical solution in each case, the 
settlements, contact pressures, bending moments, and shear 
forces at all the elements are calculated. However, for 
the purpose of presentation, more stress is laid on certain 
quantities as mentioned below, in which a geotechnical or 
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a structural engineer may be interested most. 

The maximum settlement is of particular interest 
and settlement refers to maximum settlement unless otherwise 
stated. For this purpose of settlement, it is expressed as 

P= (4.29) 

*S 

where is the maximum settlement influence coefficient 
for footing Of finite rigidity. The differential deflection 
is also a quantity of importance and it is expressed in the 
dimensionless form as 

Differential Deflection Factor 

Maximum Differential Deflection x 

= 5 ^ (4.30) 

( 1 - V ) a 

^ S'' Hay 

The results of the present work corresponding to 
elastic portion are compared in Figs. 4.3(a) to (g) v/ith 
available elastic solutions for circular footings . An 
agreement within a range of about 3.5 percent is found to 
exist in all the cases. 

The various aspects of the results obtained from the 
present work are discussed below in an order convenient for 
comparative study. 

4.4.2 Footings of Finite Rigidity Resting on Semi-Infinite Soil 

4,4 .2,1 Effect of plastic yielding of soil on bending 
moments and shear forces 

Except for an extremely flexible footing, the contact 
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Pressures obtained by elastic theory are very high near the 
edges of the footing (Brown, 1969). The soil near and under 
the edges is bound to yield. Under increasing loads, more 
and more portion of soil under the footing yields and in 
the limit the entire soil in contact with the footing 
yields so that the contact pressure becomes unifom and 
equal to qy- the yield stress of the soil. Thus the 
magnitude as well as the distribution of contact pressure 
depends on the relative level of loading. This gives rise 
to very significant changes in moments and shear forces in 
the footing under varying load. In some cases it would be 
highly erroneous to calculate these quantities on the basis 
of elastic theory alone, neglecting yielding of soil. 

(a) Uniformly loaded footing 

Pig. 4, 4(a) shows the value of critical (absolute 
maximum) moment, in a circular/ring footing loaded uniformly, 
against values of q^^^/q^ (load level). This moment is 
expressed in non-dimensional form as a . The value 

of non-dimensional moment reduces with increasing load and 
equals zero at failure (q^^/qy= 1). 

A similar trend is seen for the critical shear force 
expressed in nondime ns ional form as plotted against 

load level qav/q.y (5‘ig« 4.4(b)). The changing nature of 
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moment and shear with load level can be brought out better 
by plotting the nondimensional moment and shear in terms 
of and since in these terms q.y and a are 

constants for a given soil and geometry of footing. 

Fig. 4.4 shoTOthe nondimensional moment and shear 
plotted versus load level (q^y/ly) for a unifomly loaded 
annular footing. It is seen that the nondimensional moment 
for a relatively rigid circular footing (K = 10) increases 
from zero to about ,037 at qg_^/qy = .65 and then reduces 
to become finally zero at = 1, Ihe corresponding 

maximum value for K = .1 is about .011 at Ig-v^V = .87 and 
for n = .4 and K = .1 it is about .0065 at qg_^/qy=.8. A 
similar trend is seen for nondimensional shear, 

TiVhen the loading intensity is less, the contact 
pressure distribution is more like elastic, with very high 
values near the edges. This distribution tends to give 
higher values of moment and shear but the load itself is small 
and so the values of or are small. With 

increasing load, they increase but the contact pressure 
distribution also changes. With more and more region being 
plastic under increasing load, the contact pressure tends 
to become unifomn and equal to qy everywhere. Therefore, 
in the limit, the footing is loaded with qy and the contact 
pressure becomes also equal to qy resulting in no moment or 
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shear anywhere in the footing. In between the extreme 
conditions associated with no load and failure load, maximum 
values are obtained. These are the design values of 
nondimensional moment and shear, Pig, 4.6 shows the design 
moment values indicated by and plotted as dimensionless 
quantity ]%/qy ^ versus the radius ratio, n, for various 
values of relative rigidity , K. In this figure, the curves 
for h/a = °° apply to the case of footing resting on semi- 
infinite soil. Similarly the design shear can be obtained 
from the curves for h/a = °° in Pig. 4,7. 

The critical bending moment obtained from elastic 

theory for a circular footing with K = 10 and Vg=.5, 

2 

=,25 would be given by MQ/qg_^ a =.07. In elastic 
treatment of the problem, this value is independent of the 
load level and hence the maximum value of which is the 
design value is given by = qy a , Prom elasto- 

2 

plastic analysis for this data is (Pig, 4.6), M^=.037 qya , 
which is only about 53 percent of that based on elastic 
theory. Table 4.1 gives some values of the design moment 
and shear as percentages of those based on' elastic theory 
alone. 
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TABLE 4.1 : DESIG-B MOMENT AND SHEAR VALUES BASED ON ELASTO- 
PLASTIG ANALYSIS EXPRESSED AS PERCENTAGE OF 
VALUES BASED ON ELASTIC THEORY 

(Loading-unifoim, g=.5 or 0, ''p =.25,li/a;= °°) 

Design value based on elasto-plastic analysis 
as percent of that based on elastic theory 
alone 

n 

Bending Moment Shear Force 



K=10 

K=1 

K=.1 

K=10 

K=1 

K=.1 

0.0 

52.9 

60.0 

78.6 

56.4 

62.7 

69.7 

0.4 

47.6 

51 .8 

92.9 

51.8 

55.1 

59.2 


From this Table and from Figs. 4.5 and 4.6 can be seen 
that a design based on elasto-plastic analysis will be highly 
economical, more so far K > 1 which is the more likely range 
of K values in practice. In the conventional current design 
practice, the elasto-plastic behaviour is not considered, 
although it is quite significant as shown in the present study. 

From Figs. 4.6 and 4.7 it can also be seen that the 
design nondimensional moment and shear values are less for 
lower values of K and higher values of n. For example, the 
values of M^/qy. a^ for n=.4 are .0295 and .0215 for K=10 
and K=1 respectively and for n=.8 they are .014 and .005 for 
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K=10 and K=1 respectively. 

In all the cases it has been found that out of the 
radial and tangential moments, the latter is critical and 
that it occurs at the inner edge (centre in case of a circle). 
Except for circular footing (n=0 ), in which case both these 
moments are maximum and equal, at the centre, the maximum 
radial moments are less than the maximum tangential moments. 

This is evident from the typical values shown in Pig. 4.8. 

It has been found that the maximum radial moments are much 
smaller than the maximum tangential moments for higher values 
of n. It can be said in general that the maximum tangential 
moment is critical in all the cases. 

It has been found that, for the case under consideration, 
the critical shear occurs near mid-width of the footing (near 
half -radius in case of circular footing), 

(b) loading at the inner edge of ring footing 

Por the purpose of comparative study of various 
loading configurations the footing is divided into annular 
rings of equal width. The programme developed for this thesis 
has provision for loading these rings in any desired proportion. 
However for brevity, the 5 cases mentioned in Section 4.2,2 
are discussed herein. 
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Fig* 4.9(a) shows the effect of width of uniformly 
loaded area on the maximum dimensionless shear and moment, 
the loading being on the innermost portion of a ring or 
a circular footing. The values shown in this figure are 
lav/% ‘V'alue of about ,45 but a similar trend is found 
to exist under other load levels values). As the 

width of loaded area reduces, the maximum dimensionless 
moment and shear increase at increasing rate. For example, 
the maximum moment in terms of M^/qa,^ a is about .07, ,32 
and ,8 for load covering 100 percent, 50 percent and 10 percent 
radius respectively of a circular footing with ■ K= 10 . 

A situation of very narrow loaded area at the inner edge or 
centre should therefore be avoided. 

Fig. 4.9(b) shows the effect of positioning of 
loaded annular area of width equal to one-tenth of outer 
radius, on maximum dimerisionless moment and shear. For this 
figure qa_^/qy. is about .4»but values for other .load levels are 
found to show similar trend. The maximum moment and shear 
have small negative values for loading at the outer edge. 

These values reduce to zero and then increase on the positive 
side as the loading moves from the outer edge to the inner 
edge. The values of F^Aav example, are about -. 15 , 

.75 and 3,8 for circular footing with K=10 with loading at 
the outer edge, centre of width (half radius ) and inner 
edge respectively. 
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In both the cases mentioned above, the effects of 
load position on maximum moment and shear are less for 
higher values of n and a little less for lower values of K 
(Fig. 4.9 ). 

An extensive study of maximum moment and shear 
under elasto-plastic conditions of soil revealed that they 
remain fairly constant at all load levels e.g. Fig. 4.10. 

This behaviour is explained as below. 

Under smaller loads, the contact pressure is much 
higher near the edges, whereas under higher loads, because 
of plastic yield of soil, the contact pressure tends to 
be uniform. The effects of these two factors oppose each 
other in changing the maximum bending moment and shear force 
and hence the change remains small (within 10 percent), 

A study of a large number of cases has led to the 
conclusion that the footing under this type of loading could 
be either designed on the basis of elastic solution or 
on the basis of complete yield condition of soil amounting 
to uniform contact pressure of intensity The latter is 

recommended being easier to deal with. 

The critical shear force is calculated easily at 

any radial distance r as equal to the unbalanced force outside 

this radius, divided by the perimeter (2'n:r) at this section 
(vide Fig. 4.11). The force due to the contact pressure is 
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calculated easily on the basis of uniform stress of intensity 
q.y( yield stress of soil). The critical shear occurs at the 
outer edge of the loading and this would be the design shear. 
The design moment occurs, as tangential bending moment, at 
the inner edge (centre in case of circular footings) and 
can be calculated by assuming uniform contact pressure jqty 
and then using, and superinposing if necessary, standard 
elastic thin plate solutions (e,g, Timoshenko Krieger ( 1959 )), 


The design procedure suggested above, amounts to 
replacing soil by a uniform pressure q^. and the moment and 
shear become independent of K. 

(c) Other positions of annular loading 


An extensive study under elasto -plastic conditions 
revealed that the critical moment and shear remain fairly 
constant or slowly increase as load increases. Therefore, 
the design procedure given in (b) above is satisfactory. 

For loading on a thin annular ring, it can be assumed, without 


much loss of accuracy, that this loading is replaced by 
an annular line support. For example, the design moment for a 
ring footing loaded on a thin annular ring at the outer edge 
can be found out fairly accurately from the fonnula 


(Timshenko- Krieger, 1959) 
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The above solution is for a thin annular plate 
simply supported at the outer edge and loaded uniformly 
with intensity 

For annular loading at the outer edge ,the maximum 
shear occurs at the inner edge of the loading ap-d the 
maximum moment (tangential) occurs at the inner edge, both 
being negative. 

For annular loading at the centre of width of a 
ring footing (at half radius for circular footing) , the 
maximum shear occurs at the outer edge of the loading. 

The radial and the tangential moments are maximum and 
nearly equal under the loading. 

In all types of loadings, the design moment and 
shear force values reduce with increasing value of n. This 
aspect is discussed in detail later. 

4. 4.2. 2 Blast 0 -plastic settlement :-Bf feet of relative rigidity 

The settlement is expressed as (Eq. 4.29) 

'lav®' j 

p = ■‘■f 

B 

Except for an infinitely rigid footing, the 
settlement under a footing is non-uniform and refers to the 
maximum settlement under a footing unless otherwise specified. 
The settlement has been generally found to be maximum at 
the inner edge (centre in case of circle) except in some 
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cases like annular loading at the outer edge in which case 
it is maximum under the load, 

For footings with finite rigidity, the differential 
settlement is an important design consideration and it is 
indicated by (Differential Deflection Factor) given in 
Eq. 4.30. 

In Fig. 4.12 are given the plots of and D^ for 
two values of K (K = 10 and 1) and for uniformly loaded 
ring footing with n = 0, .4 and .8; the curves corresponding 
to h/a = °° being applicable for the case under consideration. 
The settlement (indicated by I^) is more for lower values 
of K in the lower range of values of q^y^qy. For example, 
at values of for n = 0 are 1,0 and 1,1 

for K = 10 and 1 respectively (Fig, 4,12), This difference 
is small and it is due to sagging of a flexible footing 
near the centre. 

In case of a flexible footing, the contact pressure 
is uniform- like and hence under increasing load, the 
pressures increase somewhat equally without causing yielding, 
giving rise to very slow increase in If. The curve for If 
therefore goes somewhat straight until failure develops 
over appreciable area after which settlement increases at 
a fast rate. On the contrary, for a comparatively rigid 
footing, because of highly non-uniform contact pressure, 



yielding spreads even under small load increments, giving 
rise to a drooping curve for The settlement in this 

case increases in a uniform manner* This explains the 
difference in nature of curves for K = 10 and K = 1 in 
Pig . 4-.12, This difference in behaviour due to K is 
less dominant for higher values of n which can also be 
seen from Pig. 4- 12, This is mainly because of already 
lower values for higher values of n. 

The differential settlement, indicated by 
is however found to be highly sensitive to K (Pig. 4.12), 

Por example, considering the load level qg_^/qy=.5, the 
value of for n =.4 is .0165 for K = 10 as against *11 
for K =5 1 , for n = 0 it is ,0185 for K = 10 as against 
.15 for K = 1, and for n=.8, it is ,0025 for K = 10 as 
against .01 for K = 1 . As the load approaches the failure 
load, the differential settlement under a uniformly loaded 
footing tends to vanish since a situation of uniform 
contact pressure equal and opposite to uniform loading 
on the footing is approached. The differential settlement 
under a uniformly loaded footing is therefore most critical 
under working loads. 

Pigs, 4.13 and 4.14 show the plots of and for 
three types of annular loadings, for two values of K (10 and 1 ) 
and two values of n (0 and .8). Apart from n, the maximum 
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settlement is controlled by the load position, relative 
rigidity K, and the load level indicated by lav^^' 

Referring to Fig. 4..I3 ( n= 0 case) and considering the 
effect of K on it is seen that for loading at the outer 
edge, is more for lower value of K and the difference 
increases rapidly beyond > ,3. For example, at 

q.av/'ly ^f for K = 1 as against 2,1 for K = 10, 

In the other two cases of loading, is less for K = 10 
than for K = 1 upto lav^'ly ~ vice versa beyond 

lav'^^y ^ * 5 * O^^e of the effects of load position on 

is already mentioned above (loading at the outer edge), 
l£ is less sensitive to load position, at any load level, 
for higher value of K as seen from Fig, 4,13 where the 
corresponding curves are close to each other. 

Prom Fig. 4.13 it is seen that is low and fairly 
constant for all load positions and levels for higher value 
of K. Hcwever it can be seen that B^ is highly sensitive to 

load position for lower values of K, For K = 1 , B^ at qa-^/qy=. 4 - 
is about 4.7, 2,1 and 1.5 for loading at the inner edge, 
centre of width and outer edge respectively. The corresponding 
values for K = 10 are , 6 , .3 and . 2 , This shows that 
loading at the innermost portion is critical from the point 
of view of differential settlement, more so for lower values 
of K. 
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Fig. 4.14 is similar to Pig, 4.13 except that it 
is for n=,8 instead of n = 0, It is seen from this Pigure 
(Pig, 4 . 1 4 ) that is sensative to load position for lower 
value of K (K=l) and the rest of the curves show less 
sensitivity at all load levels to load position or K. It 
is also seen that in this case (n=.8), I)^ is minimum for 
load position at centre of width. 

Considering I^ and at all load levels it can 
be said that the position near the centre of width of ring 
footing (near half radius for circular footing) of an 
annular loading is better than any of the other positions 
of load. This applies particularly, for relatively more 
flexible footings. 

4 . 4.2, 3 Effect of radius ratio n 

Prom Pigs, 4.3 to 4.6 and 4.8 to 4.14, it can be 
seen that the maximum or design moment and shear values 
decrease with increasing values of n, the reduction being 
faster as n approaches unity, Por the sake of comparison, 
some design values for n=,8, obtained from the present work 
are quoted in Table 4.2 in terms of percentage of corresponding 
values, for a circular footing after applying the correction 
(Section 3.4.3 and Pig. 3.5). The comparison stands for 
equivalent footings i.e. for equal total loading on the ring 
and circular footing and equal area for ring and circulax 
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footing. The design shear for loading o'bliGr than uniform 
loading is calculated as explained in section 4.,4:,2.l(b) and 
Pig. 4.11. 

TABLE 4.2 : TIPIOAL LBSIdE MOMENT AND SHEAR VALUES POR A 

RINO POOTING WITH n=.8 EXPRESSED AS PERCENTAGE 
OP G0RRESP0iE)ING VALUES POR CIRCULIR FOOTING, 
UNDER EQUAL AREA AND TOTAL LOAD CONDITIONS 


Quantity Loading K Percentage 


Design shear 

Uniform 

10 

10 


Uniform 

1 

13 


A 

Any value 

20 


B 

Any value 

U 


C 

Any value 

36 

Be sign moment 

Uniform 

10 

100 


U nif orm 

1 

25 


G 

Any value 

36 


Jlnnular loading of one tenth width of ring 
footing at inner edge (loading within .1 a 
of circular footing). 

Annular line loading at centre of width of ring 
(at .5a for circular footing). 

Annular line loading at outer edge. 


Note : A : 

B : 

C : 
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From the above Table, it can be seen that the 
design moment and shear values are considerably less for 
the ring footing in almost all the cases. 

From Fig. 4. 1 3 and Fig. 4.14 and Table 4.2 it can 
be seen that the settlement of footing is less for higher 
values of n for any K, position of load and load level 
q.av'^'ly after applying correction for equivalence 

of footings. The differential settlement, represented by 
Djp in these figures, is also seen to be highly sensative 
to n, being very much less as n approaches unity. For 
a rational comparison, the quanity D^/unit width of footing 
should be compared for equal total load. To account for 
this, for hing footing/D^ for circular footing is to be 
multiplied by Y( 1 -n )/(l-n) (Fig, 3 . 5 ). If this is done, 
it is found that the real advantage is obtained if n is 
greater than about ,6. Referring to Fig. 4,11 it can be 
found Out, after applying the correction, that 

Dj/unit width for n=.8 

= ,154 = 15.4 percent 

D^/unit width for n=0 

for equal total load with K = 10 and ^or K = 1, 

this percentage works out to 20 percent. This is a 
substantial advantage in favour of ring footing. 

Thus ring footings are better suited than circular 
footings for carrying the same total load. This is true 
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from the point of view of moment and shear in the footing, 
as well as settlement and differential settlement, 

4.4.3 Footing Resting on a Soil layer of Finite Depth 
Underlain by Rough Rigid Base 

As for the case of infinitely rigid footing in 
the previous Chapter the only change due to finite depth comes 
in through the use of the appropriate elastic solution. 

The entire discussion given in Section 4 . 4,2 above 
for semi -infinite soil holds good, in general, for the case 
under consideration. The discussion below is therefore 
aimed at only a comparative study involving the effect of 
depth h, of the layer. 

The design moment and shear values are shown in 
Figs. 4.4 and 4.5 for two values of h/a, three values of K 
against n for uniformly loaded footing. It is seen than 
those values are lower for lower values of h. For example, 
the values of design moment as percentage of those for 
h/a = °° are, for n = 0 , K = 10 , about 87 percent and 69 percent 
for h/a = 2 and 1 respectively. These percentages v/ork out 
to about 80 percent and 50 percent for K = 1, Thus for 
lower values of h the design moment and shear values are 
less. This effect is less prominent for higher values of n 
(narrow ring footings). The design values drop considerably 
for very low values of depth of layer. 
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Under loadings other than uniform loading, an 
extensive study of maximum moment and shear has revealed 
results similar to those for h/a = «>. Therefore, the 
design values be better based on uniform contact pressure 
qy. Thus for the loadings considered (other than uniform) 
the design moment and shear values would be independent 
of K and depth ratio h/a also. iPor finding the design 
values, the procedure given in Section 4.4.2.1(b) and Fig. 

4.11 is to be followed. 

The settlement and differential deflection are 
expressed as in Eq. 4.29 and Eq. 4.30 respectively, Fig. 4 . 1 5 and 
Fig. 4.1 6 show • the plots for and for uniform loading 

for various parameters and two values of h/a. The curves 
for h/a = 0 ° shov/n in Figs. 4 .l 3 and 4.14 can be used 
for comparison. From these figures qt can be seen that 

is appreciably affected by h/a. For example, for qg_^/qy=.5, 
for h/a = 1 expressed as percentage of value for h/a;=°° 
is about 33 percent for all the values of K and n shown 
in the figure, (This percentage is however not constant for 
all values of K and n at all other values of qg^^/qy). It 
is observed that the effect of h/a on is significant and 
it remains so at all load levels. For uniform loading the 
differential settlement represented by is also highly 
sensitive to h/a. For example, B^ for n=0 and qg_^qy=.4 
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expressed as percentage of for h/a = °° is about 87 percent 
for h/a = 2 and 55 percent for h/a = 1 . The corresponding 
values of for n = ,4 are 81 percent for h/a = 2 and 
45 percent for h/a =1. B^for higher values of n which is 
already low even for h/a = is further reduced very 
much with decreasing values of h/a. 

Pigs, 4.15 and 4. 16 show plots for and for 
two values of h/a, the loading being at the inner edge on 
one-nineth width. The curves for uniform loading are 
also shown for the sake of comparison. Since loading on 
one nineth width presents considerably different situations 
for vaxious values of n, a direct comparison with respect 
to n is not advisable. Broadly it can be said, referring 
to Pigs. 4.15 and 4.16, that is less sensative to loading 
configuration and load level for higher values of K, lower 
values of n and lower values of h/a. is high for 
loading at the inner edge even for low values of h/a, 
except for n values of about .8 and above. 

In short, the effect of radius ratio n as discussed 
for h/a;= °° is similar for finite values of h/a. As mentioned 
earlier, the design nondimensional values of moment and 
shear under uniform loading are lower for higher values 
of n. Under other loading positions, h/a does not affect 
these design values (since the design is to be based on 
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uniform contact pressure equal to q ) and the discussion 
for h/a = °° holds good. Both and are reduced, the 
latter considerably with increasing values of n, and 
even considering equivalence, narrow ring footings will 
give very small values for design moment and shear, as 
well as settlements and differential settlements, in 
comparison to circular footing. All these values reduce 
further for lower values of h/a. 

All the design values discussed above depend 
on q^/and since qy for finite values of h/a, is to be 

corrected as explained in Section 3.4.5 for thickness of 
layer will be lower than those for h/a = 

4,4.3 Effect of Embedment of Footing in Semi-Infinite Soil 

For this case Mndlin's solution has been used in 
place of Boussinesq's solution. The behaviour agrees in 
general with that of a surface footing except for some 
differences which are discussed below. 

Fig, 4.17 shows the bent portions of the footing 
near the centre of a uniformly loaded circular footing 
embedded at depths of .3, .6 and 1 times the outer radius 
(c/a =,3, .6 and 1). It is seen that the curvature is 
sharper for c/a =.6 than for c/a =,3 or c/a = 1, (Although 
settlements are less for higher values of c/a). This type of 
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behaviour is found to exist for all values of K, n and 
lav^^y* bending moments are very sensitive to the 

deformed shape of the footing (Bqs. 4.5, 4.25). It is 
therefore found that the dimensionless design moment 
(Mg^/ly ) i'S maximum for a certain value of for 

c/a values between a range of .5 to .6 for uniformly loaded 
embedded footings* Pig. 4.18 gives these majximum galues for 
K = 10, 5 and 1 and n = 0, .4 aud .8, It is seen that 
these values increase as c/a increases, become maximum 
between c/a = .3 to c/a = .6 and then decrease at a slower 
rate as c/a increases, beir^ almost constant for c/a > 2, 

All the above discussion holds good for dimensionless 
design shear (P^j^/ciy a) also and the maximum values are 
shown in Fig, 4.18. As pointed out in the earlier Chapter, 
the value of q^y appropriate to c/a should be used. 

For loadings other than uniform loading, it is 
again found that the design moment and shear should be 
found from the basis of uniform contact pressure qy, on the 
lines given in Section 4.4.2.1(b) and Fig, 4.11. Thus for 
such loadings, the design moment and shear values are 
independent of depth of embedment, 

I^ and follow the same overall pattern as for 
surface footing, the values being lesser with increasiiig 
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"Value of depth, of embedment. However, the curves, if 
plotted as for earlier cases do not show a simple pattern 
of variation with ly* It is suggested that the 
correction be applied to values found for surface footings 
as given in Chapter 3 for rigid footings. 

4,5 OONOLUSIOHS 

(1) A study of axisymmetrically loaded annular footings 
of finite rigidity resting at or below the surface 
of a semi -infinite soil or resting on a compressible 
soil layer of finite depth has been made using thin 
plate theory for the footing. The Bending Moments and 
Shear Forces in the footing and the settlements are 
studied under various loading configurations and 
under elasto-plastic soil condition and with tension 
cut-off wherever applicable. 

(2) It is found that differential settlement and 
design moment and shear under uniform loading 
are highly sensitive to relative rigidity K of 
footing, Maximam and differential settlements are 
higher and design moment and shear are lower for 
lower values of K. 

(3) It is found that in all the cases, a ring footing, 

^ especially with inner radius more than about ,6 

times the outer radius gives lesser design bending 



moment and shear force and settlements and are 
advantageous over a circular footing of the same 
area and carrying the same total load. 

Under uniform load, it is found that the footing 
need be designed for values much lower than those 
obtained from elastic theory alone*. Design charts 
for some values of parameters are presented* 

Under annular loadings, it is found that the design 
moment and- .shear in all the cases should be found 
by assuming uniform contact pressure equal to yield 
stress of soil. 

For the loadings considered, the locations of 
maximum moment and shear are investigated. Under 
uniform loading the maximum shear occurs near the 
centre of width of ring footing or near half radius 
of circular footing. The maximum moment occurs at 
the inner edge (or the centre) except for annular 
loading at half width in which case the maxim-um 
moment as well as shear occur near/under the load. 

Out of the maximum radial and tangential bending 
moments, the latter are found to be critical in all 
the loading cases considered. 

The design moment, shear and settlements are analler 
for lower values of depth of layer or for higher value 
of depth of embedment. 



The settlement of embedded footing can be 
approximately fdund by applying correction procedure 
given for rigid footing (Chapter 3). 

A design procedure has been given for evaluating 
the design moment and shear in the footing and the 
maximum and differential settlements of footings of 
finite rigidity resting on or in semi-infinite soil 
or on soil layer of finite thickness. 
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TYPE-u TYPE-q 


I 



TYPE-b TYPE-c 

u - Uniformly distributed load 

Q - Annular load at inner edge 

b - Annular load at centre of width 

c - Annular load at outer edge 
w - Width of ring 

FIG.4-1 TYPES OF LOADING CONFIGURATION 
ON FOOTING OF FINITE RIGIDITY 
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(a) Elements in the neighbourhood of i 



Fictitious elements l,2,m+3,m4-4 


(b) Numbering of plate elements 


EDGE 



Fictitious ^ , 
elements 


(c) Elements in the neighbourhood 
of edge 


FIG. 4-2 PLATE DISCRETIZATION 
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FIG.4-3b CONTACT PRESSURE BENEATH FIG* 43 c DIFFERENTIAL DEFLECTION 

CIRCULAR RAFT yp = .3 IN CIRCULAR RAFT yp:^-2 
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FIG.4'3d MAXIMUM MOMENT IN CIRCULAR 
RAFT >J>=-3 

r — BROWN (1969) 

-- - PRESENT WORK 



r/a 


F!e.4-3e BENDING MOMENT DISTRIBUTIONS 
IN CIRCULAR RAFT Vp=*3 






EGOROV ET AL 

(1963b 

POULOS(1968) 
PRESENT WORK 



FIG AOf RADIAL BENDING MOMENT DISTRIBUTION 
IN UNIFORMLY LOADED RIGID CIRCULAR 
RAFT. >>g=0-4, Vp=0-25,h/a=1 



FIG.4'3g RADIAL BENDING MOMENT DISTRIBUTION 
FOR CONCENTRATED LOAD WITHIN r=0-1a 
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FIG.4-4 MAXIMUM (Critical) BENDING MOMENT 
AND SHEAR AT VARIOUS LOAD LEVELS 
IN A UNIFORMLY LOADED FOOTING 






FIG.4-5(a) DEPENDENCE OF CRITICAL B-M. ON 
LOAD LEVEL 
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FIG'4-6 DIMENSIONLESS DESIGN MOMENT FOR 
UNIFORMLY LOADED RING FOOTINGS 








s is 



FIG. 4-8 COMPARISON OF TYPICAL MAXIMUM VALUE: 
OF RADIAL AND TANGENTIAL BENDING 
MOMENTS 
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Inner loaded width 7<> 

F10.4-9(a) EFFECT OF WIDTH OF UNIFORMLY 
LOADED AREA (innermost portion) 

ON MAXIMUM SHEAR AND MOMENT 



h/osco 

yp=-25 
— K -1 


\ nzO 

V 


A 


1 2 3 4 5 6 7 8 9 10 

Loaded annular ring (number 1 for 
innermost ring) 

FIG-49(b) EFFECT OF POSITION OF 
SINGLE LOADED ANNULAR 




FI6.4-10 MAXIMUM MOMENT AND SHEAR VARIATION 
WITH LOAD LEVEL FOR INNERMOST ONE 
NINETH WIDTH LOADED 



FIG’4-11 METHOD OF CALCULATING DESIGN SHEAR 
FOR RING FOOTING LOADED AT INNER EDGE 




FIG-4-12 DESIGN CHARTS FOR SETTLEMENT 

INFLUENCE FACTOR AND DIFFERENTIAL 

deflection FACTOR FOR UNIFORMLY 
LOADED RING FOOTING 

(a) K=10,n=0 (b) K=10,n=-4 




FtG.4-12 DESIGN CHARTS FOR SETTLEMENT 

INFLUENCE FACTOR AND DIFFERENTIAL 
DEFLECTION FACTOR FOR UNIFORMLY 
loaded RING FOOTING 
(c) K= 1Q,,n=-8 (d) K= ^ 






qav/<^y 


Flfs-M2 DESIGN CHARTS FOR SETTLEMENT INFLUENCE 
FACTOR AND DIFFERENTIAL DEFLECTION 
FACTOR FOR UNIFORMLY LOADED RING FOOTING 

(e) K = 1,nr-4 (f) K = 1,n=-8 




FIG-4-13 MAXIMUM SETTLEMENT FACTOR AND 
DIFFERENTIAL DEFLECTION FACTOR 
FOR VARIOUS LOADINGS ON FOOTING 
ON SEMI-INFINITE SOIL 
Loading on one nineth width 
Inner edge-Q,d ^Centre of width-b,e 
Outer edge-c,fiK=10 for a,b,c jK=1 or d,e,f 
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fio.4-u maximum settlement factor and 

DIFFERENTIAL DEFLECTION FACTOR 
FOR VARIOUS LOADINGS ON 

FOOTING ON SEMI-INFINITE SOIL 
Loading on one-nineth width : 

Inner edge-a^dj Centre of width-b^e 

Outer edge- c,f iK=: 10 for a,b,CjK = 1 for 
d,e,f 






u 







“Uniformly distributed load 

“One nineth width at inner edge loaded 
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FIG.4-17 SETTLEMENT OF EMBEDDED CIRCULAR 
UNIFORMLY LOADED FOOTING RELATIVE 

TO CENTRE OF FOOTING 
( rrdistance from the centre) 




144 



4-18 VARIATION OF DESIGN MOMENT AND SHEAR 

WITH DEPTH OF EMBEDMENT FOR UNIFORMLY 
loaded FOOTING 



OHAPTER 5 


EIAS TO -PLASTIC AHALYSIS OE CEITRALLY 
LOADED PILED RIGID ANMAR FOOTINGS 

5.1 INTRODUCTION 

In the design of the foundation for a heavy structure 
on a deep deposit of clay or for a medium structure on a 
deep deposit of soft clay, it may be found that excessive 
settlement and not bearing capacity failure is the deciding 
design criterian. In order to restrict the settlement 
it is advisable to provide a few short piles. In such 
cases, traditionally the piles are designed to cej^ry all 
the load, which is not logical. In fact in the pile-raft 
system, the aim is to take best advantage of the load 
bearing capacity of the raft and the settlement restricting 
capacity of piles. An integrated analysis, taking into 
account the load sharing between the raft and piles and 
the extent of plastic yield in each component of the system 
at all load levels, is justified for predicting the behaviour 
of the system. 

Pile raft foundations have been reported (Poulos 
and Davis, 1977) to have been successfully uaed in Mexico 
(Zeevaert,l957), and London (Hooper, 1973). 



A brief description of the various approaches used 
for analy Gin,','; the ])ilo-raft systems has been given by Poulos 
and Uavi, G (iPBO), wliich includes the comparative study made 
by Bro.wn et al, (1975). The method of 'Plate on piles and 
continuum' (Hain, 1975) was assumed by 'Brcwn et .al,(1975) as 
reference solution since it involved the least approximation 
and the elaciticity based method by Poulos and Davis was 
reported to be comparatively the most satisfactory from 
settlenuitit point of view amonf^ the methods considered 
(whicli excluded f iuiit e-element analyses). In simplified 
finite element anal.yses (Hooper, 1973, Desai et al.l974), the 
piles liave boon replaced, in an approximately axially- 
symmeti’ic problem, by an equivalent continuous annulus. 
However, choosing an appropriate stiffness for the annulus 
simulating'’: the piles is difficult, A complete three- 
dimensional analysis can give a satisfactory solution, but 
the cost involved, if a generalised or parametric study is 
aimed at would be definitely prohibitive, A method which 
considers the interaction at all load levels between piles, 
raft and soil; the nature and extent of plastic yield at 
all the contact points; is able to predict directly the load 
settlement behaviour throughout the' entire load range upto 
failure and which is not computationally expensive, is 
therefore needed. An attempt is made in this study to develope 
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such a iiicthod , altliotigli some assumptions are made regarding 
the I'igidity of the raft and piles. 

For a proper load transfer from above the footing to 
the piles, it is necessary and advantageous to have the 
f ootin^'^ of relatively high rigidity. If the footing is 
quite flexible, a load, if it is located above a pile, 
will not be shared by the footing and if it is not above 
the pile , it will not be transmitted to the pile. In 
this study, therefore, the footing is idealised to be fully 
rigid. Also, relative to the heavy loading and/or soft soil 
for which the pile-raft combination is intended, the piles 
can also bo treated, in the first instance, as rigid. Hence 
the piles .'ire assumed incompressible. If necessary , correct ions 
for compreosibility of piles (and flexibility of footing) 
can be applied as outlined by Poulos and Davis (1980), 

The advantiiige of annular footings over circular ' 
footing with regard to settlements is brought out in the 
earlier Chapters and if they are provided with even a 
small number of short piles, the settlement is likely to be 
reduced considerably. _ Hence piled annular footings are 
ideally suited for many heavy tower Hke structures. An 
integrated elasto-plastic analysis for piled annular footings 
is not available presently in literature. An integrated 
analysis becomes all the more necessary when the piles are 
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a.ri*3Jig0d ill inoiro "tlian one concenliric cicle and would be 
having different loads and plastic slips. 

Piled annular footings are expected also to exhibit 
excellent resistance to eccentric loading. However, this 
aspect has not been tackled in the present work. 

Basically, numerical integration of available elastic 
solutions has been employed in this study. Hence the size of 
thti problem, for the axisymmetric cases considered in this 
Chapter, In terms of size of matrix is equal to the number 
of annular rings in w?iich the footing is discretized plus the 
number of elements in which dissimilar piles are discretized. 
Therefore the size of the problem still remains relatively 
small, and the computer programme developed for the purpose 
of this thesis can bo used to analyse any specific case 
at a very low computational cost, 

5 . 2 ASSUllPTIOHS iU® STaTEiffim OE THE PROBLEM 

5.2,1 Assumptions 

In addition to all the assumptions contained in 
Section • 5 . 2 . 1 (Chapter 3), the following assumptions are 
also rmulu for the purpose of the study in this Chapter. 

( 1 ) The soil is semi -infinite in extent. 

(2) The raft rests on the piles and the surface of the 
soil and is perfectly rigid and smooth. 


(3) The piles are incompressible. 

(4) The pile-raft system and the loading form an 
jsLX isyiimetr ic jarrcmgement . 

( 5 ) -Ihe shear strength of the soil is constant with 
depth. 

(6) The inclusion effect of piles is neglected. 

( 7 ) Displacement compatibility is satisfied only in 

the vertical direction. 

5 , 2,2 Statement of the Problem 

A rigid rinii; footing of external radius 'a 'and 
internal radius 'na* rests simply on incompressible piles 
and on the surface of a semi-inf init e soil mass (Fig, 5.'!). 
The pile-raft system and the central loading form an 
axisymrautric arrangement. All the piles in a system are 
astnunod, for the sake of simplicity, to have the same 
size (Length- L, diameter = d). 

5 . 3 ANALYSIS 


5.3.'! General 

In addition to the discretization of the footing, 
the piles are also discretized (Fig, 5.2(a)). Each pile 
is divided into a suitable number of elemen-us of equal 
Imipbh plus an element representing the base of the pile. 
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"I 

whor© ( P ) {p) is the non-dimensional displacement 

a 

1 

vector, { q } ~ { q} is the non-dimensional contact stress 

vector and [I] is the influence coefficient matrix, {q} 
includes contact pressures at the bottom of the footing, 
the strusDos along the periphery of the piles and the 
stresses at the bases of the piles. 


5,3.2 Influence Coefficients 

The settlements or displacements of various elements 
are found out as below. The stress on an element, is assumed 
to act aa a concentrated force at the geometric centre of 
the element except for finding the influence of an element 
on itself in vrtiich case, the mathematically integrated 
formulae are used. 

(a) Inflmmce of footing elements on footing elements 

Although Boussinesq's solution is applicable, 
Mindlin's solution with depth equated to zero is used for 
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the salce of generalisation (Poulos and Bavis, I 974 ). In 
general, the settlement at a point in the soil mass due to 


a 


point load P (Pig. 5.2(b)) is given by 


^ (1 + v) ^ 3 _ 4 v 


) 


R 


1 


+ 8(l-v )-(3-4v) + 

Ro 


p d ' 

(3--4v) (a+c) -2cz 602(2+0) ^ 

+ ) 


l>3 






( 3 - 0 )' 


R^ 


(5.3) 


where z ia the det)th of the point under consideration 
below the surl’acse, c is the depth of the point below the 
surface where the jjoint load is acting and R-\ and R 2 are 
distances, as given below, in terms of c, 2 , and the 


hori; 5 ont;al distnnoe between the two points under consideration, 
r (Pig. 5.2(b)), 
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(5.4) 


where dA..: is the area of the element 3 . 

il 

Por the particular case under consideration, 


c = 0 I 

2=0 I 

= R 2 “ ^ t 


(5.5) 




I'lie influence coeff icients are obtained by using 
ISq. 5.::5 for various elements in the form of Eq. 5.1 and 
finally of Eq. 5.2. 

For 1=3, Eq. 5.3 is not valid and the Mndlin's 
solution integrated over a rectangular area is used assuming 
the small element to be rectangular. 'The settlement at the 
centre of the element is obtained using this solution 
(E’oulos ami Davis , 1974) as 
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(5.8) 
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wliGX’B X tind y "tlio l6ngt;h,s of shorter and longer sides 
of the area respectively. 

(b) Irifluence of footing element on pile element 

For the whole pile-raft system, it is rational 
to consider the centre of a pile element as the centre 
of the cylindrical element of the pile. Eqs. 5.3 and 5.4 
are used in which, for the particular case under consideration, 
c « 0 

z '-S depth of centre of pile element, 
r - horizontal distance between centre of the 
footing element and the centre of pile, 

(c) Influunco of pile element on element of another pile 

■^Qs. 5.3 and 5.4 are used in which, for the particular 
case under consideration, 

z ~ depth of centre of pile element on which 
influence is being found out, 
c = depth of centre of pile element whose influence 
is bein{f found out, 

r — horizontal distance between centres of piles 

6.Jk ~ Tcd . dll where 

d = diameter of pile and 

dL == length (height) of pile element,. 

dA = TLd^/4 , for the base element. 
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(d) Ijif. lufcjnce ot pile 6161116111; on footing element 

Bqs. 5.3 and 5.4 are used in which, for the 
particular case under consideration, 
a -- 0 

c = depth of centre of pile element, 
r = horizontal distance between centres of pile 
and footini5 elements 

dA wi.-dL as in (c) above 

2 

dA -- TOi‘"/4 for tip element of pile. 

(e) Influence of pile element (other than base) on pile 
eleminit of the same pile 

For this case, the pile element surface is divided 
into subelcjmonts (Fig. 5.2(c)) and the concentrated force 
on the subolemonts has been integrated mathematically to 
got the influence of the pile element on an another element 
of the same pile, as 

(1+ V )q=d 

P ^ ^ 

^ HK ( 1 - V ) 

a % 

- 4(1 - V )ln tan ( £ 4^ ^ 

4^5 

- (3-4 ) sin P ~ ^ 



8 (l-v) In tan( 2 4 ^ 


(5.9) 
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2(3-0) ^ 2 (z+c) 

whore a =- tan” , P = tan” — - 

d d 

tlL dL 

0^ c — , Oo “ c + -TT- > 


O ’ ^2 


c ai 


2 = dex)tii of centre of pile element on which 

influence is being found out, 
id cJL' are depth to centre and length (height ) of pile 


elument whose influence is being found out respectively, 

(f) J'ld’luence of base element of pile on an element of 
the sjAinc pile 


For this case, the base element circle is divided 
into suboleraents (Fig, 5.2(d)) and the concentrated force 
on the subelememt has been integrated mathematically to get 


the influence as 



!il.. [C3-4v) + 8(1-v^) TB 

4l!)(1-v) 

- (3-4v) ys - (z-c)^/ Y'A-(3-4 v) (z+c) /fB 

2 .^/2 

+ 2cz/'fB--2cz(c+z) B"^'^ ] 


(5.10) 


where A = (o-z) + r'^ 

/ x2 2 ^ 

B = (c+z) + r ; 

r = d/2 for the upper limit j 

r = 0 for the lower limit j 

c ~ depth of pile base below sui'face 1 — length of pif 
z =s depth of centre of pile element on which I 

influence is being found out. 
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(g) Influence of base of pile on itself 

fhio is obtained by using the same formula derived 
in Section 5-5.1 (f) with z = g= L, but multiplying by the 
factor %/A- to account for the rigidity of the base assuming 
that this factor true for surface area is approximately 
applicable for embedded area also (Poulos, 1977). 

5.3.3 Pinal Form of Equations to be Solved 

Since the footing is assumed as rigid and the 
piles as iriconipi'esslble, the settlement will be equal 
to say, Pq* Hence Eq. 5.2 can be written as 

{IJ} .. (1- v^) [l], {q} (5.11) 

where p ~ {U} is unit vector. 

w KJ 

For overall equilibrium, 

Z q^ = Q (5.12) 

where Q is the external load and the summation is taken 
over all the elements of the footing and the piles. 

5.3.4 Elastic Case ' 

Equations 5.11 and 5.12 can now be solved to obtain 
unknown stresses {q} as well as Pq’ 

5.3.5 Elasto-Plastic Case 

Even under very small load, the contact pressures 
near the edges of the footing and the tips of the pile may 
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reach a very high value giving rise to yielding of the soil. 
Under higher loads, more and more soil elements yield 
and as per the assumption made, the stress is to be held at 
the yield value of the soil, for such elements. 


In case of a pile-raft system, the yield value 


of the soil for the footing and the pile will be different. 

In the present study, the yield stress along the pile is 
assumed to remain constant. The yield stresses are as below . 

(a) For the footing = 5.69 c^ 


(b) 

(c ) 


For the pile shaft 
For the pile base 


'lyb “ ^ °u“ 


5.69 

9^yf 

5.69 


for a = 1,0 


Under increasing load, as soon as the value of 
stress at any element reaches the corresponding yield value 
as given above, it is held constant at that value and 
partitioning and partial inversion of, modified influence 
matrix is done as given in Chapter 5. 


The settlement is expressed in the following form 


where I is the settlement influence factor for the pile- 

IP 

raft, Q is the total load on the system and d is the diameter 


of a reference pile. 
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A quantity of particular interest is the load 
sharing:; between the footing and the various sets of piles. 

A set of piles for this purpose is one having piles that 
carry equal loads due to their concentric placement and 
such a set will be denoted as a 'pile-ring*. 

The piles in different rings will not only carry 
different loads, but they will be at different stages of 
failure under any load. Therefore interaction factors found 
out on the basis of equal loads on two piles are not 
applicable. Hence for a more rigorous analysis an integrated 
pile-raft analysis is necessary and the same is adopted 
in this work. 

5.4 WiiiULTS Al© DISCUSS lOI 
5.4.1 fkm^ral 

For a combined pile raft system, the load level is 
indicated rationally by the ratio of total load on the footing 
to the total ultimate load, (and not by qa.Y'^%- 
case of footing alone). The ratio Qy/Q indicates overall 

factor of safety. 

calculated on the basis of yield-stress being 
reached at all the points of contact between the soil and 
the pile-rirft system is given by 

0^ - TtA'qyf + mTtdl qyg + UybA 


(5.14) 
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Where ra is the number of piles, A» is the net area of the 
footin^,^ (a' = A-rnmi /4) , lyf* lyg and are yield stresses 
as defined in Section 5*3.5 and A is the area of the footing. 

Although the computer programme dereloped in this 
thesis allows for unequal diameters and lengths of piles, 
the dicinui'ter and lerif^th of all the piles in a system are 
assumed, for the pui'pose of simplicity, as equal. 

The load shr.3red by a pile ring expressed as a 
percentage of bhe total load is denoted by 'p'. The value 
of is assumed as .5 throughout. 

It was sufficient to divide the pile into a small 
number of elements (Poulos and Davis, 1977). Depending on 
the length, each pile was divided into 5 to 20 elements. 

5.4.2 Single Floating Pile in Semi-Infinite Soil 

Some results have been obtained for the purpose 
of studying the settlement response of piles under various 
load levels upto failure. Fig. 5.3 shows the case of a 
single pile, the load-settlement curves for which are drawn 
as P/ c^d*" versus PE/c^d, where P is the load on the pile, 
c^ is the ultimate shaft resistance (assuming a = l) and 
d is the diameter of the pile. It can be seen that, 

(i) !I!he load-settlement curves are almost linear 

upto about 55 percent, 68 percent and 75 percent 
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o,t ultimate locid for L/d values of 9.2, 18 and 25 
respectively. Thus the linear- like response is 
over a larger range of load level for higher values 
of L/d. 

(ii) Beyond the linear- like range, the curves become 

nonlinear at increasing rate as the load increases. 

In Pig. ‘3.3 , the load-settlement cur^fes obtained 
from field tests for a stiff clay by O'Neill and Reese 
( 197 ^) foi' two viilues of L/d are shown. For drawing these 
curves the values of E and c^ are obtained by back-figuring. 
It is assumed that the soil is semi-infinite in extent, A 
curve obtained from PBM (Besai, 1977 )for the case of L/d 
= 9.7 is also shown in which the back-figured values for 
the corresponding field test are again used. These curves 
also show a high non-linearity beyond the initial linear 
portion, and a rapidly increasing settlement as failure 
is approached, A large number of field tests such as those 
by Whitaker and Oooke (1966) and Burland et al. (1966) show 
the same trend, although the tests are in different types 
of soils. Recently Boek Madsen and Lagoni (I 98 I) have 
reported from a field test on pile in whit-e chalk that the 
settlement varies in cubic proportion to the load. The 
present analysis gives results conforming to these trends 
and the agreement as seen from Pig. 5.3 with the results 



from field test and from PEM, Is satisfactory, thus ralidating 
the analysis, k load-settlement curve for L/d = 25 based 
on the method of Poulos and Davis (I 98 O) is also shown in 
Pig, 5 . 3 . it agrees well with the curve obtained from the 
present imalysis upto a load level of about 76 percent of 
ultimate load. 

5 . 4,3 Effect of Radius Ratio n on Settlement Response of 
Pile-Ri Footing? System 

i?i{?:ure 3.4 showy the load settlement responses of 
three footijvjs with n =- 0, .6 and , 8 . The outer radii 
of the tlireo f ootlif'ns are selected in such a way that the 
area of the footings for all the cases remains the same. 

Each of the footin^i is provided with a ring of piles 
consiotini*: of six identical piles and hence the footings can 
be treated as equivalent systems. Prom Pig. 5.4 he seen 

that the settlement response improves with n under all load 
levels (denoted pj Q • For example for Q/Qu=* 2 the 
values of Ip (which indicates settlement) are about ,036j 
.029 and ,025 for n =0, .6 and .8 respectively and for 
QA1^“.5 , the values of Ip are about .087, .0 68 and .05 
for n-0, ,6 and ,8 respectively; thus higher n being more 
advantageous for higher values of It is better to 

adopt a value of n as high as possible sub^t'^t to considerations 
of availability of space in plan and structural design of 
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the ring* In the seme figure (Fig, 5.4) » the settlement 
response expresvsed in terms of Ip for an equivalent footing 
with no piles has been shown; equivalence being attained by 
designing the unpiled footing also for the same total ultimate 
load 0^ as for piled footing. Even though a high value of n 
is chosen (n=.8), it is seen that the settlement of unpiled 
footing? is much more than that of piled footing. For example, 
Ip foi‘ piled footing with n=.8 and six piles of diameter .15a 
and L/d value of 25 is about 19 percent of that of unpiled 
footiiif'i at Q/CL=*2, this percentage being about 53 percent 
at Q/a^=.7. 

The percentage pile load (relative to total load 
on the system), p, is shown in Fig, 5.4 for three equivalent 
systems. It is seen that p is higher, at all load levels 
except at failure, for higher values of n. For example, at 
Q/Qu = .5| p is about 50 percent for n = 0 and about 63 percent 
for n=.8. Thus for higher values of n, the participation 
and utilization of piles is better. It is also seen that 
in comparison to working load range, the pile load would be 
overestimated if based on elastic analysis and under- 
estimated if based on failure load. In any case it is 
illogical to design the piles for the entire load. 
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5.4.4 Effect of Positioning of Pile-Ring 

Pi^pAre 5.5 shows the effect »of positioning a 
pile ring ,on the settlement (indicated by I ) and the load 
shared by piles ( indicated as percentage p) for two cases 
of n ( n=0 and .5)* For n = 0 case, d the diameter of pile 
is 0,075 a and there are 6 piles with L/d =25. For n=.5 
case, d = #075 a and there are 6 piles with I/d = 30, 

The position of pile ring is defined by ’s’, (Fig, 5.1), 
such that 

Rixdius of pile ring- Inner radius of ring 

® Width of ring(=outer radius-inner radius) 
r„- na 

= -J (5.15) 

(1-n) a 

The settlement of the system (Fig, 5.5) in the 
lower range of load level (Q/Q^) is a, little lower for 
higher value of s but in the higher range of Q/Q^, it is 
lower for smaller value of s (Pile-ring near inner edge). 
This difference in settlements for various values of s, is 
very low for higher values of n except near failure. 

The minimum settlement in the range of Q/Q^ values 
greater than ,25 occurs roughly for s =.75 and .5 for n=0 
and ,5 respectively. It can be said that if the pile 
ring is placed at a position dividing the footing area 
inside and dutside into equal regions, minimum settlement 
is obtained (although p is not a maximum). For narrow ring 
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f octings "til© pil© icing can. b© placsd ab "tlis csniic© of 
width of ring. This would be advantageous from the point 
of view of positioning the load also near the centre of 
the width and causing minimum design bending moment and shear 
in the footing. 

In the same Figure (Fig. 5.5), the values of p 
(percentage load shared by piles) are plotted for the case 
of n~,5, for 3 values of s (s=.1, .5 and ,9). These curves 
join together at failure but at low values of Q/Q^^ the 
difference is considerable for various values of s. For 
example, at Q/'.J^~.3, p values are about 42 percent, 50 percent 
and 59 percent for s = *1 , .5 and .9 respectively. However 
the settlements at this Q/Q^^ are almost same for 3 values 
of s. It appears that the load sharing between the piles 
and the footing and the distribution of contact stresses 
in various cases produce some balancing effects on settlement 
for various values of s. 

5.4.5 Combined Effect of Length and limber of Piles in . 

Nearly Equivalent Systems 

One can thinic of alternative systems of piled 
ring footings, keeping the size of the footing, the diameter 
of piles and the number of piles times L/d the same. Such 
systems will be nearly equivalent since their ultimate 
load bearing capacity will be approximately equal. 
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Figure 5.6 shows the load settlement response of 
two sets of such systems. For the sake of general comparison 
the dominance or otherwise of piles can be represented by 
percentage pile load at failure, denoted by p^. The values 
of p^ are about 76 percent and 4.6 percent for the n=.8 
systems and n=. 5 systems respectively. For lower value of 
p^, the settlement response for the two equivalent systems 
(L/a =10, 3 piles and l/d=5, 6 piles) is almost idLentical, 
whereas for the higher value of p^, the Ip for L/d = 40, 
and 3 piles is lower than that for L/d = 20 and 6 piles, 
beiri^,^ about 67 percent at Q/Q^^ =.3 and 50 percent at 
q/Q^=, 7. Since the cost of piles per unit length increases 
with increasing? depth, a cost analysis will be required to 
decide which of the equivalent alternatives is economical 
in the lori^? run for the sets of equivalent systems having high 
pile loads. For lower values of pile load one can say that 
it would be preferable to use an equivalent alternative having 
more number of shorter piles than having less number of 
longer piles provided the piles are not very near to each 
other. 

5.4.6 Ring Footing with Two or More Pile Rows (Rings) 

It may be necessary to provide piles in 2 or. more 
rows if the number of piles desired to be provided is large 
causing too less spacing between piles. It may also be 
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convenient to provide the piles in more than one row 
even though they could be accommodated in a single ring. 

The piles in different rings will carry different loads 
and will be at different stages of failure. 

Figure 5.7 shows the effect of positioning piles 
in two rows. The pile radius is .06a where a is the outer 
radius of the footing. For each system, there are 12 piles 
with L/d =5* Ihe three systems are equivalent. From this 
figure it is seen that the best positioning of the two 
pile rings (rows) is at the inner and outer edge. In the 
same figure , the response for 3 pile rings with 4 piles 
each, two beini^ placed near the edges and one at the centre 
of width (in staggered angular location) is shown. It 
shows that this arrangement is, in an overall way, better 
than that of placing 2 pile rings, with 6 piles each, near 
the edges. 

Some of the values of load carried by piles are 
given in Table 5.1 for the above cases. 

From the Table 5.1 , it can be seen that, the 

load shared by the pile ring near the inner edge is less 
and that for pile ring near the outer edge is more in all 
the cases, except at failure. 
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Ti3LE 5.1 : LOAD SHARING BY PILE RINGS 


System 

Number 

Pile-Iiing 

Position* 

No. of 
Piles in 
the Ring 

s 

Load Shared (Percentage 
of total lo^) 

At Q/C) 

= .2 

At Q/Q^ 
= .67 

At Q/ 

= 1 

1 

Outer 

6 

.9 

28,5 

19.1 

13.5 


Centre 

6 

.5 

17.5 

16.9 

13.5 

2 

Outer 

6 

.9 

29.5 

19.3 

13.5 


Inner 

6 

.1 

14.5 

15.7 

13.5 

3 

Centre 

6 

.5 

18,8 

17.0 

13.5 


Inner 

6 

.1 

14.2 

15.0 

13.5 

4 

Outer 

4 

.9 

19.0 

17.0 

9.0 


Centre 

4 

.5 

11.8 

11.8 

9.0 


Inner 

4 

.1 

7.8 

8.2 

9.0 


Provision of more pile rin^^with staggered arrangement 

is the best arrangement. 

5.5 OONOLUSIONS 

(1) An integrated elasto -plastic analysis of a rigid ring 
footing supplemented by incompressible piles in one 
or more annular rows or rings, is presented, 

(2) Providing even a few short piles below a footing reduces 
the settlements considerably. Ihe b.^nefit of the use 
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of the combination of advantages of the piles and 
the raft will be lost if either of the component 
dominates. 

( 3 ) It is highly advantageous to increase n (keeping 
the area same) of a piled rijog footing. 

(4) To provide more number of short piles or less number 
of lon^j piles (keeping number of piles x L/d=constant ) 
makes little difference from settlement point of view 
unless the piles share a high percentage of load (more 
than about 50 percent). Hence the former alternative 
may be preferable. For pile dominating systems a 
cost analysis will be necessary. 

(5) It is preferable, from the point of view of overall 
settlement response, to place the piles in as many 
number of pile rings as possible (subject to 
possibility of accommodating in the width of ring 
footing). In such cases, angular staggering of 
piles in various rings in essential. 

( 6 ) Charts for settlement response as provided in the 
present study can be used for the purpose of design* 
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RING FOOTING 
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IS fbuhd out 


(a) PILE DISCRETIZATION (b) MINDLIN PROBLEM 


•B 



Ic) INFLUENCE OF PILE ELEMENT (d) INFLUENCE OF BASE I 
ON ELEMENT OF SAME PILE ELEMENT OF SAME Pll 

FIG.5-2 APPLICATION OF MINDLIN’S SOLUTION 
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0 4- 8 12 16 


f E /Cud 

FIG.5‘3 LOAD SETTLEMENT RESPONSE FOR 
SINGLE FLOATING PILE 





All piles L/<1=25 
Pile dia = -15a 

FIG.5-4 SETTLEMENT FACTOR AND PERCENTAGE PILE LOAD FOR 
EQUIVALENT PILE RAFT SYSTEMS 
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Al-n*0,S:^'5 1 

A2-n»0,S=-75 ' 

A3* n -0 1 S" *875 
B1-n=:5,S^*1 ] 

B2-ns5,S»*5 ? 

B3-nts5jSs*9 J 


6 Piles 
d=*025a 
L/cl=25 

6 Piles 
d-*075Q 


FIG* 5*5 EFFECT OF POSITIONING OF PILE RING ON 
SETTLEMENT AND PILE LOAD SHARE 
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FI0.5*6 COMPARISON OF NEARLY EQUIVALENT 
SYSTEMS OF PILED RING FOOTINGS 
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•08 
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*16 

0 2 *4 *6 *8 1*0 

Q/Qu 

FIG.5-7 EFFECT OF POSITIONING OF PILE 
RINGS 
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ELAS3!0-^"*LAS!I!IC AIALYSIS OE SETTLSIISIE} OP 
BQGEHTRIGALLY LOADED RIG-ED AMULAR POOTIRGS 

6.1 IMDRODIJGTIOR 

Annular or circular footings which are well suited 
for tall tower like structures are generally centrally loaded. 
Howe'^rer forces s\ich as wind load may cause eccentric loading 
on the footing. Because of the tall nature of the structure, 
a high eccentricity under such load can not be over-ruled 
at least occassionally. Also for such structures, the 
footing has to be strong (rigid) enough. Hence it is 
appropriate to study the settlement response of eccentrically 
loaded footings which are highly rigid, 

I'he e last 0 -plastic analysis of centrally loaded 
footings has been made in Chapter 3 and the same is extended 
in this Chapter for eccentric loading. Because of eccentric 
loading, rotation is caused. This, and the possibility of 
lift-off under high eccentricity on the tension side of the 
footing are the two important aspects introduced because of 
eccentricity of loading. 

Ihren though soil can be treated generally as an 
elastic continuum, it would be erroneous to assume that it 
could take tension. This aspect of tension cut-off has to 



177 


be therefore taken into account in any rational and realistic 
analysis of settlement. 

For a rigid circular footing the elastic solutions 
for moment loadings are given by Borowicka (1943) and Sneddon 
( 1946 ). The solution given by Yegorov and Eichiporovich 
(1961), although approximate, is of particular interest since 
it is applicable to a footing O.n a finite layer of soil. 

Both these solutions do not take into consideration the 
plastic yielding of soil or the lift-off or separation, 

Ihere is no solution available presently for eccentrically 
loaded rigid annular footing taking these factors into 
account. 

In the present study, the horizontal forces, if 
any, on the footing are not considered. These forces are 
not expected to contribute appreciably to the vertical 
settlement and the effect on the rotation is considered in 
the form of eccentricity. 

6.2 ASSUMPTIOIS AM) STATEMFIT OF THE PROBLEM 
6.2.1 Assumptions 

In addition to the assumptions contained in 
Section 3.2,1 (A) and (B) (Chapter 3) , the following 
additional, assumptions are.made f or ■the purpose : of 
study in this Chapter, 
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(1) The external force acting on the footing is vertical 
but not central. 

(2) V/herever tensile stresses tend to develope at 

the contact face between the footing and the soil, 
the stress is held at zero value, 

6*2.2 Statement of the Problem 

A smooth rigid annular footing of outer radius 
*a' and inner radius 'na' carries a vertical load Q at an 
eccentricity ’e ' (Pig. 6.1). 

Three types of problems as in Section 3.2.2 
(Chapter 3) are considered, the only difference being that 
the load is eccentric. 

6.3 AiNlAhYSIS 

6.3.1 General 

The analytical procedure outlined in Chapter 3 
and the portion of Chapter 4 involving tension cut-off 
are generally applicable. The elastic solutions mentioned 
in earlier Chapters are all applicable, 

6.3.2 Elastic Case 

The footing is discretized into small elements 
(Pig. 3.1) and the stress q.^ over a small element is assumed 
uniform and acting at the centre of the element as a point 
load, except for the case when a formula is available for 



finding the influence of an element on itself. 


The settlement at the centre of element i 

due to stress on element, j , is expressed (Sq. 3.4) as 


1~ v' 


q^ a 


^13 


( 6 . 1 ) 


13 J JJ 

This equation in matrix form for all the elements 
is written as ^q, 3.5), 


{p > 


(1- v*^) [I] {q} 


( 6 . 2 ) 


where {p} and {q> are non-dimensional settlement and 
contact stress vectors respectively. 


Since the footing is assumed to be rigid, the 
settlement at any point at a distance x from the centre 
can bo expressed as 

p j= a X + k (6.3) 

where x = x/a , ic is the non-dimensional vertical 
displacement of the centre of the footing and cc is the 
rotation of the footing about its centre. 

The settlement-contact stress relationship in 
matrix form (Eq. 6,2) along with Eq. 6.3 can be written as 

(1- /) [I] li) = (0 S + 6 ) 

For overall equilibrium of forces and moments 


on the footing, 
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X- dii.^ ~ Q ( 6 . 5 ) 

1 dA^x^ — Q, e (6,6) 

where dA^ is the area of element 3 , and is the distance 
of centre of element j from the centre of the footing 
measured in the direction of eccentricity, and Q is the 
vertical force acting on the footing at eccentricity e. 

If m is the total number of elements into which 
the footin^^ is divided, m+2 equations contained in ( 6 . 4 ) 
to (6.6) above are solved to yield m values of contact 
pressures and also the two unknowns a and k which give the 
settlements. The solution can be written in the following 
form 

Q 1-v^ 

p ~ ( 0x + 6) (6.7) 

3* ■ E 

where ^ and 6 are rotation and settlement influence factors 
respectively. 

6.3,3 Elastic Case with Tension Out-off 

If the eccentricity of the load is high enough, the 
solution of Eqs. 6.4 to 6.6 above would indicate some negative 
values of contact pressure on the other side of eccentricity. 
Since soil can not take tension, such a solution would stand 
invalid. To avoid this situation, the contact pressures 
over all such elements are assigned zero value (Fig. 6.2), 
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and "ths ©equations are resolved, keeping in mind ■feh.a't 
compa'fcibili'ty of displacements for these elements need 
not be satisfied. It may turn out that after the redistri- 
bution of contact pressures obtained after the solution, 
some other elements may indicate negative contact pressures, 
for which the procedure is to be repeated until satisfactory 
solution with no tensile contact stresses is obtained. 

6.3*4 Blast o-Plastic Case Without Tension Out -off 

Bven under a small load, the contact stresses 
near the edges of a rigid footing are extremely high and 
the soil under some of the elements of the footing would 
yield. Let suffixes e and y indicate elastic and plastically 
yielded values respectively, Eq, 6.2 is converted, for 
convenience of the analysis, into partitioned form as below 
(same as Eq, 3.10, Chapter 3). 

t = (1.v2) (6.8) 

Py % 

In the above equation, q^ indicates the value of 
yield stress of soil to which, as per assumption made in 
this analysis, the contact stress is held (Big, 6,2) after 
yielding takes place. (qy=nondimensional yield stress =q/E). 

Si.nce compatibility of displacements is satisfied in 
the elastic region only, Eq, 6.8 above is written as 
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{p^)= [I.] 

V 

wb.er6 * 1 is the modified coefficient matrix corresponding 
to the rows of elastic elements only ( Sq. 3.11 ). 

The equilibrium equations 6.5 and 6.6 can be 
re-written, after dropping the suffix j for brevity, as 

E q^.dA + Eq^.dA - Q (6.10) 

i: q^.dA .x+iq^.dA.x = Q.e (6.11) 

where subscripts e and y denote elements in elastic and 
yielded ssonos respectively* 

Equations 6.9 through 6.11 can now be solved, along 
with Eq, 6.3 to get the values of contact pressure in the 
elastic region (the value being known and equal to q^. in 

yielded region) and the unknown settlement components in 

the form of central settlement and rotation. It may turn 
out that the solution thus obtained contains contact pressure 
values higher than yield stress q^. In that case, the 
procedure is to be repeated until satisfactoiy solution with 
contact pressure not exceeding q^ anywhere is obtained. 

6.3.5 Elasto-Plastic Case with Tension Out-Off 

Under higher values of eccentricity and at any 
loading (even small) both lift-off and plastic yielding of 
soil would take place. After applying the conditions, namely, 
for any clement 
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h “ 'ly ’ ^ 1y 

(6.12) 

and 

Q.Q "" 0 y if* K 0 > 

(6.13) 


it may turn out, after the first trial solution using the 
techniques given in Section 6.4.3 and Sedtion 6.4.4 above, that 

(i) some contact pressures exceed and /or , (ii) some 
contact pressures are negative. 

‘Ihe contact pressure values at these elements are 
made to take up the values as in iiqs. 6.12 and 6.13 above 
and the solution is obtained aigain. If the required conditions 
are not found to be satisfied, the procedure is repeated until 
the conditions of Jiqs. 6.12 and 6,13 are nowhere violated. 

6.4 RBSIJLTS Am DISCUSS lOIT 

6,4,1 Minimisation of Computational Effort 

At the outset, only half of the footing symmetrical 
about the di;ameter on which the resultant external load acts 
is considered. 

The restribution of contact pressures at the beginning 
of each trial solution, along with the involvement of two 
components of settlement and two restrictions on the contact 
pressure values is expected to run into a long iterative 
procedure. However, it was found that the zone of lift-off 
(or separation) marches with an almost straight front towards 



the centre of the footing in all the cases considered. 
Thexefore it was possible to adopt the following procedure 
and reduce the iteractions normally to 2 to 3 only. 

Equation 6.4 is split as below, without loss of 
generality • 


(1- v^) 

[I] 

{ C-j 

q+C2l}={0x + ^ } 

(6.14) 

(1- v^) 

[I], 


= {X} 

(6.15) 

(I^v2) 

[1 1 

^^4 

11 

(6.16) 


where , {!]'} is unit vector and c-j , C 2 , and c^ are 
constants such that 

c-^ = 0 ; c,| = C 2 6 and c^ + Cj = 1 (6.17) 

For convenience, a beginning is made with an 
assumed extent of tension cut, -off or zone of separation. 

For the plastically yielded zone, the procedure given in 
section 6,4.4 is adopted. In case some negative values of 
contact pressure occur, c-j and C 2 are simply chosen in 
such a way that, while Eqs. 6.17 are satisfied, the recombining 
of Eqs. 6.15 and 6,16 would give all positive contact pressures 
and exactly zero pressures on the boundary of the separation 
zone. Partitioning and elimination of the yielded and 
separation zones from the computations reduces the matrix 
size successivel 3 r under higher loads and/or eccentricities. 





A computer run for the entire range of loading 
and eccentricities under given geometry of footing took on 
DEC 10 system o^ily about 2 minutes of CPU time. 

It may be- mentioned that, even allowing for the 
diametric symmetry, a PEl solution of the three-dimensional 
nature as involved in the present problem, even under simplified 
elasto— plastic or non-linear elastic analysis involves 
enormous computations in terms of computer memory as well as 
OKI time. 

6.4,2 Footing Resting on the Surface of Semi-Infinite Soil 
6. 4.^.1 Cruneral 

Pigu.rG 6.5 shows t-'/pical load settlement responses 
of centrallv and eccentrically loaded footings. The particular 
case of e = 0 of the present analysis is seen to tally well 
with experimental and PSM .results given by Desai (1968) 

(curves 1 to 3» Pig. 6.3) • I^ the same figure, some load 
settlement curves for typical values of radius ratio n and 
eccentricity e ( e = e/a ) are also shown (Curves 4 to 8 
Pig. 6.3). For e other than zero, the settlement refers 
to the maximum settlement of the footing (which occurs, 
at the outer edge of the footing). It is seen that the 
settlement response for n=,9 remains nearly linear upto 
about 85 percent, 70 percent and 60 percent of failure load 





for e = 0, .2 and .3 respectively and then shows a tendency 
of rapid increase to failure. The settlement response for 
n=0 , on the other hand shows a steadily increasing trend. 

For comparing the performance under equal load, the settlement 
obtained under the same load level for a ring footing with 
a certain value of n is to be multiplied by the Squivalenoe 
Factor (Section 3.5.3)) given by1/Y(l-n ), Using this and 
the curves of Fig, 6. 3 , it is seen that the settlement of 
ring footing with n =.9 , at 50 percent of ultimate load, 
is about 44 percent, 28 percent and 24 percent of that for 
n=0 (carrying; the same load and having the same area) for 
e = 0, ,2 and .3 respectively. This shows that the settlements 
of a rin ^5 footing are much less than those of a circular 
footing at aJLl eccentricities and load levels. 

From the elasto -plastic analysis it is found that 
the settlemcjnt components indicated by central settlement 
factor 6 and rotation ^ depend on eccentricity, radius 
ratio and load level. 

pf , b = f (n , e , q^y/Uy ) (6.18) 

The results of the present investigation are 
presented in compact form to serve as design charts (Figs, 6. 4 
to 6,7), In these charts, the values of 6 and 0 are given 
against i (»e/a) for various values of n and qy/lav* 



settlement at any point x from the centre measured in the 
direction of eccentricity can be found out from Eq. 6.7 , 
namely, 

Q 1 - M ^ 

p = — — (0x + 6) 

a E 

where x = x/a and 0 and 6 are to be read directly from the 
design charts (Eigs, 6.i to 6.7). These charts can be used 
in a number of ways. For example, for a given size of a 
ring footing, the permissible load and/or eccentricity for 
a desired ratio qy/q^y (Factor of Safety) can be worked out 
using the charts and the. above equation. 

The revsults from the present analysis, in the lower 
range of load level, nearly represent elastic behaviour and 
the value of 6 in this range is found to be equal to .5 for 
circular footing (same as in Chapter 3) and the value of 
<0/e for circular footing allowing tension is found to be 
.75 which tally with solutions by Borowicka (1943) and 
Yegorov and Nichiporo.vich (1961). 

A discussion about the settlement response is made 
below on the basis of the design charts (Figs. 6.4 to 6-.7), 
with reference to various parameters, It ought to be mentioned 
here, that since the parameters are not explicitly operative 
the discussion is only meant for broad comparative study. 



6.4.2. 2 Sffect of radius ratio n 


Since the values of 0 and 6 are to be used in 
£q. 6.7, a compaj:ison at the same load Q and the same load 
intensity between a circular and a ring footing can be 
made by applying the Equivalence Factor as below. 

The outer radius for a ring footing equivalent to a 
circular footing (same area) is given by (Section 3.5.3) 

b = a/Y'(l-n^) 


Therefore, the settlement of ring footing equivalent 

to circul£yo footing of radius a (carrying the same load and 

having equal area) can be written as 
Q 

b- 


(l-v^) X 

^ + 6 ) 
E D 


2 

where b = a/Y(l-u ) . 

The maximum settlement is given by (x=b) 


max 


Q ) 


(0 + 6 ) 


(6.19) 


( 6 . 20 ) 


This can be written as 


Q (1-v ) 


max 


E 


^^eq ^eq^ 


where 0eq " 0Y(1-n^) and = 6f(1-n ). 


( 6 , 21 ) 


Thus the 0 and 6 for an annular footing are to be 

multiplied by Y(l-n^) to get equivalent values for circular 


footing. 


oome of "bhe 6 qui valent values worked out using tlie 
above relations and tlie design charts are siven in Table 6,1 
assuming that the values of non-dimensional eccentricity 
(e/a or e/b)are the same, 

TABL-fi 6.1 : BQUIV.ALEI'JT GEMTRAL SSTTIEIffiET AM) ROTATIOM 
GOEFEIGIEETS (qy/q^v = 5) 


n 

b/a 

e/b 

Equivalent values 

.. ^eq+ ^eq 

/eq 

^eq 

0.0 

1 .00 

.1 

0.20 

0.60 

0.80 



.3 

0.40 

' 0.63 

1.03 



.5 

1.70 

1.10 

2.80 

0.8 

1.67 

.1 

0.06 

0.39 

0.45 



.3 

0,20 

0.42 

0.62 



.5 

0.70 

0.40 

1.10 


From the Table it can be seen that the central 
settlement, rctation and maximum settlement, indicated by 
6eq_, 0eq respectively are much less for n=,8 

than for n=0. This differende will be more if e instead of 
e is taken same for the ring and the circular footings* 

A study Cf design charts (Eigs. 6-. 4 to 6,7) shows 
that ring footints are advantageous over circular footings 


from the consideration of settlement and rotation for i 
values of e/a and q„/q___. 

6. 4- 2. 5 Effect of eccentricity 

Prom the charts for / in Pigs. 6.4 to 6.7), it' is 
observed that 0 (rotation coefficient) increases with e/a 
at any load level, as expected. The increase in 0 with e/a 
is very rapid after a certain eccentricity for low load 
levels and it is so for even for a small value of 
eccentricity for higher load levels. Both these types of 
behaviour are explained on the basis of elasto -plastic 
nature of soil. The contact pressures which are already 
very high at the edges increase further on the compression 
side under eccentric load even if the load is small in 
magnitude. This induces larger yielding of soil under 
increasing eccentricity causing more rotation. In addition 
to this, tension cut-off is created on the tension side 
which also accentuates rotation. Under loads close to 
failure, a large portion of the soil remains yielded due to 
which the rotation is sensitive to even small changes in 
eccentricity, 

Por the same reason involving plastic yielding 
of soil, the central deflection also increases with 
eccentricity as seen from Pigs. 6.4 to 6.7. However, this 
alone does not give the correct idea of the movement of the 
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foo'fcing, since "blie seiilementis must be looked upon as the 
sum of and 6, Thus the following quantities ane of 
particular interest. 

Maximum settlement coefficient = 0+6 (x/ap=1) (6.22) 

Minimum settlement coefficient + 6 (x/a=-l) ( 6 . 23 ) 

Differential settlement coefficient =20 ( 6 . 24 ) 

Distance of point of rotation from the centre 

= ^ a , ( 6 . 25 ) 

These quantities can be obtained from the design 
charts. It can be found that the centre of rotation moves 
from infinity (no rotation) to smaller distances from the 
centre under increasing values of e/a. 

Some typical values obtained from Eqs. 6.21 to 6.25 
above are given in Table 6.2 after converting to equivalent 
values as given in Section 6.4.2. 2. 

Prom Table 6.2 it is also seen that at large values 
of eccentricity and/or load levels (qg_y/qy), even a physical 
lift -off -the-ground can take place (indicated by negative 
values of minimum settlement). 
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TABLE 6.2 : TYPICAL 7ALUES 0? MXII/UM AI® MliTOUM ODLPPIGIEI® OP 

SETTLEMlDlfE AM) DISTANCE OP POINT OP ROTATION 


n 

ly/ lav 

e/a 

CO 

^eci 

Maximum 

Settle- 

ment 

Coeffici- 
ent = 
^eq’^^eq 

Minimum 
Settle- 
ment 
Coeffi- 
cient = 

“^eq"^^eq 

Distance 
of Point of 
Rotation 
from Centre 
= (6e^/52^eq)a- 

0 

4 

0.0 

0.00 

0.51 

0.51 

0.51 ■ 

oo 



0.3 

0.23 

0.57 

0.80 

0.34 

2.480a 



0.5 

1.00 

0.75 

1.75 

-0.25 

0.750a 

0 

3 

0.0 

0.00 

0.52 

0.52 

0.52 

OO 



0.3 

0.40 

0.63 

1.03 

0.23 

1.575a 



0.5 

1 .70 

1.10 

2.80 

-0.60 

0.295a 

, 8 

3 

0.0 

0.00 

0.55 

0.55 

0.55 

oo 



0.3 

0.20 

0.42 

0.62 

0.22 

2.100a 



0.5 

0,70 

0.40 

1.10 

-0. 30 

0.570a 


6. 4. 2.4 Marching of yielded and tension cut-off zones 

As a part of the solution, the contact pressures and 
the extents of yielded zone and tension cut-off zone (if any) 
are also obtained for any stage of loading and eccentricity. 
Typical yielded and lift-off zones are also sketched in 

i 

Pig. 6,8. The nature of typical contact pressure distribution 
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und©ir yiGlcIsd. ajid sspsxsi'tjLoii coiidi'tioiis is shown, in Pig,6«9» 
Wh©!"©©.© "til© finsi yi©ld6d son© devslops nsan ih.© ©irtii*© 
outer edge of a circular footing under central loading, 
it develops only on the side of eccentricity for eccentric 
loading. It then proceeds towards the centre with increasing 
load. For an eccentrically loaded ring footing the first 
yield starts at the inner and outer edges on the side of 
eccentricity* 

The tension cut-off zone starts on the tension 
side and under increasing load or eccentricity it marches 
towards the centre, with almost a straight front. Fig. 6.10 
shows the relationship between ly/lav eccentricity at 
inpending tension cut-off. It is seen that tension cut-off 
starts, at value of 5 at e/a values of about 0.34 

and 0.39 for n=0 and 0.8 respectively. For any load level, 
the tension cut-off starts at higher values of e/a for 
footings with higher values of n (ring footings with narrow 
width) which is in favour of ring footings. Another interesting 
observation from this figure (Fig. 6.10) is that for values 
of q /lav enough (infinity in the limit )which indicates 

close to elastic conditions, the e/a necessary to cause • 
tension cut-off is seen to be about .39 for n=.8, whereas, 
when ©last 0 -plastic analysis is made, it is found that tension 
cut-off can be caused at e/a values of about .385, .375, .36 and 
.2 at ly/q^.^ ratios of 2.5, 2, 1,5 and 1.25 respectively. 
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Thus tension cut-off is likely to occur at values of 
eccentricities lower than those arrived at from elastic 
theory. 

6.4.3 Effect of Embedment of Eooting 

For the purpose of this study, the embedment 
amounts to application of net load Q on the footing 
embedded at a depth c below the surface of semi— infinite 
soil which is basically treated as elastic medium and the 
Mindlin’s solution is applied instead of Boussinesq’s 
solution. 

Figure 6,8 shows a tsrpical contact pressure 
distribution for c/a = 0 and c/a = 5 (The contact pressures 
are almost same for the two cases). The contact pressure 
distribution as well as the behaviour of the embedded 
footing is found to be very similar to that of a surface 
footing at all load levels and hence it has been found 
possible to correlate it to that of surface footing after 
applying correction factors. The procedure for applying 
the coCTection fantors is given in Section 3.4.4 end the 
correction factors to be applied for depth and for radius 
ratio and v for central loading are given in Table 3.2 

and Table 3.3. 

For eccentric loading, the settlement of a rigid 
footing, at any point at a distance x from the centre of the 
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footing luoasurod in tho direction of eccentricity is 
expressed as 
Q 

^ a * IT" ^c^^n % W 9^' ^ ) (6.26) 


where Q 


a 



C,,,0, c= 


X = 
6 ', 0 ’= 


load on the footing 

External radius of footing 

Elastic constants of soil 

Correction factor for depth of embedment 

(I’able 3.2, Chapter 3) 

Settlement correction factors for n and v 
(fable 3.5, Chapter 3). 

Rotation correction factors for n and ^ 
x/a 

Centraf settlement, and rotation coefficients 
ob-cained for surface footing (Eigs. 4.4 to 4.7), 
but for appropriate value at depth c as 
explained in Section 3.4.4. 


The values of R^. and R^ are given in Table 6.3. 

The effect of n on Rj^ is not much ' as apparent 
from Table 6.3; however the effect of n on other parameters 
such as 0 is significant. 

The effect of depth of embedment alone on rotation is 
of the same proportion as for central settlement and is 
included through the correction factor Cq. 
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The contact pressure distribution, changes in it 
at various load levels and the nature and the extent of 
yielded and separated ?.ones for embedded footings axe 
similar to those for surface footing. Also the discussion 
for surface footings applies to the case of embedded footings 
except for the effect of the correction factors. The load 
level and eccentricities at which tension cut-off is caused 
is the same as for surface footings (see the curves for h/ap^ 
in Fig. 6.10), However here again, since the reference 
itself is higher, tension cut-off in case of embedded footings 
will be caused at higher load/eccentricity levels in 
comparasion to surface footings. 

6. 4-, 4 Footing Eesting on soil Layer of pinite Depth underlain 
by Rough Rigid Base 

6, 4. 4,1 General 

Solutions for this case are obtained by using the 
appropriate elastic solutions for finite layer, I 

The settlement in this case., is expressed as ; 

2 

p = — ^ (a V 4 -) : 

8 # 

where the subscript h is employed to indicate finite depth | 
and a is the correction factor for effect of Foisson’^s ratio 
It is found that it is not necessary to apply correction to $4 i 

Also, a s= 1 for '^ = .5 



6 . 4 . 4 . 2 Settlement and rotation response for Poisson's 
ratio = .5 

Figures 6.11 to 6.16 give tlie values of $ 2 |j^and 
for various parameters and for Poisson's ratio =.5 (a= 1 ). 

The nature of the curves is generally similar to that for 
semi-infinite soil except for the following: 

(i) Depth of layer has a considerable effect on the 

settlement and rotation of the footing, Por example, 
if 6 j;iare to be restricted to say ,5 each at 

a load level then for a circular footing, 

e/a has be less than about 0 . 01 , 0,25 and 0 . 3 I 
for h/a = 0 ° , 2 and 1 respectively. For the same 
data, if e/a is .2 , qy/lav "values are 5, 1 .C and 
1.5 for h/a = °°, 2 and 1 respectively (Pigs. 6.4 , 

6.11 and 6 . 12 ). 

(ii) The effect of higher values of n is similar to that 
for the case of semi— inf inite soil; however, this 
effect and the effect of depth of layer are 
additive giving rise to much less settlement and 
rotation for a narrow ring resting on soil of 
finite depth. Por example, the (4. value, at <ly/lav“^ 
and e/a = .2 , is about .57 for h/a;=-, n=0. The value 
is ,3 for h/a = °° , n=,8 and it is .18 for h/a =1, 

n=, 8 (Pigs. 6.4, 6.6 and 6.16). 



6.4. 4*3 JifXect of and oorrection for Poisson's raiio 


Figure 6 .17 shows t3rpical values of 0^ obtained 
directly for h/a = 2, n=0 and Poisson's ratio =.5 and 0. 

It is seen that the curves for Poisson's ratio of ,5 and 0 
are nearly the same for the case considered, for two values 
of qy/ lay* Phis trend is found in general and therefore 
it is not found necessary to apply any correction to on 
account of Poisson's ratio. Also in this Figure (Fig.6-,17) 
are shown typical curves for obtained directly as well 
as obtained by appllying the correction a (Sq. 6.27). Ihe 


value of a is obtained as follows: 

It, (1-(.5)^) 

a = X 5 

% 5 (1.. v^) 


( 6 . 28 ) 


where v= Poisson's ratio under consideration. 


I. = Value of settlement influence coefficient obtained 

ily 

from design chari^s for centrally loaded rigid 
annular footings (Figs. 3.10, Chapter 3) against 


required value of Poisson's ratio and for the 
required values of n and q^y/iy. 


% = As above for Poisson's ratio =.5 

» 5 

a = 1 for V = .5 


(6.29) 


It is seen from Fig. 6.17 that the values of 6^^ 
obtained for v = o by applying the correction as above 
tally fairly well with those obtained directly for \)= 0 



similar treiid is found to exist generally. Hence, tlie 
simplified correction procedure as above is suggested, 

V/hile finding out the value of in all the 

cases of footing resting on finite layer of soil, if the 
corrected value of q^. for depth as given in Section 3.4.5 is 
used, the value of corrected load level will be 

reduced on account of this correction and all the design 
values of settlement, moment etc. will also reduce, 

6.5 GONGIiUSIOffi 

(1 ) A simplified elasto-plastic analysis has been 

conducted to obtain the entire response of settlements 
and rotations of rigid annular footings (circular 
footing being a particular case) under vertical 
eccontric loading. 

(2) The contact pressures at the edges of the rigid 
footing which are very high under central loading 
are further increased on the compression side if 
the load is eccentric and this gives rise to plastic 
yielding of soil even at lower load levels and to 
non—lineso^ sottlement /rotation characteristics with 
increasing load and/or eccentricities, 

( 3 ) Design procedure and charts are presented for the 
three cases considered. 



Tension, cut-off occurs at eccentricity values lower 
than those predicted by elastic theory when elasto— 
plastic behaviour of soil is taken into consideration. 
A.nnular footings are found to be better suited than 
circular ones from the consideration of not only 
settlements but rotation also under eccentric 
loading. 

The settlements and rotations are less for footings 
resting on finite layer of soil and embedded footings 
in comparison to surface footings on semi -infinite 
soil. 
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F(6.6‘t ECCENTRICALLY LOADED RIGID ANNULAR 

FOOTING: DEFINITION SKETCH (Problem type A) 



T -Tension cut-oft zone 

E-Elastic zone 
y-yiclded zone 

F 10 6*2 CONT ACT STRESSES 



Mail, surfoce dt$placcfi»cf}t in x 
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Applied vertical pressure psi0i&9KN/ir^} 



FIG.6-3 typical LOAD SETTLEMENT CURVES 
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FIG. 6-4(a) ROTATION ECCENTRICITY CURVES 















FIG.6-5 (b) SETTLEMENT ECCENTRICITY CURVES 
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SETTLEMENT ECCENTRICITY CURVES 
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FIG. 6-9 TYPICAL CONTACT PRESSURES UNDER RIGID 
CIRCULAR FOOTING UNDER ECCENTRIC LOAD 
FOR c/a=0 AND c/at:5 



FIG. 6*10 LOAD ECCENTRICITY RELATION AT TENSION 

CUT-OP 





















FIG.6*17 EFFECT OF V ON SETTLEMENT FACTORS 
FOR ECCENTRICALLY LOADED RIGID 

ANNULAR FOOTING ON FINITE DEPTH OF 
SOIL 



CHAPTER 7 


OONSOLIDATIOH SETTLEMENT OE ANNULAR EOOTINGrS 

7.1 INTRODUCTION 

7.1.1 General 

In addition to immediate or final settlements, the 
time rate of settlement is an important aspect requiring 
investigation for the purpose of predicting the behaviour 
of a footing with respect to time and for the purpose of 
design. 

Although footings will have finite rigidity, it is 
convenient to study the extreme eases of fully flexible and 
fully rigid footings and these are the cases which are 
considered in this study. For a flexible or rigid footing, 
an aspect which needs due attention is the effect of period 
of construction on consolidation. Except for soils of 
very low permeability, a considerable amount of consolidation 
settlement takes place during the construction period itself. 
A study of this aspect is included in the present work with 
reference to annular footings. 

Many tj^es of soils, due to their formation process, 
are more permeable in the horizontal direction than in the 
vertical direction. Thus the effect of anisotropy also . 



needs consideration. 


7.1.2 Consolidation Under Flexible Annular Footings 

For the general three-dimensional consolidation 
analysis, two major types of approaches are available. 

The pseudo-three-dimensional approach (Terzaghi-Rendulic) 
assumes three-dimensional flow but one dimensional settlement 
caused due to flow of fluid only. This approach is comparatively 
simple and has been successfully used, to predict time -rate 
of settlements. For predicting the overall time -settlement 
response, this approach, after corrections to time scale etc,, 
has been used for three-dimensional problems and it has been 
reported (Davis, 1969) that the results are fairly accurate 
when compared to the more rigorous approach (Biot), In this 
approach a simple Finite Difference numerical scheme can be 
used to solve the governing differential equation. The Terzaghi- 
Rendulic approach has been used by Madhav and Vitkar (1981) 
for flexible circular footings with special reference to 
period of construction and anisotropy of soil with respect 
to permeability. In the present work, the same approach is 
used for the analysis of consolidation of flexible annular 
footings. 

In the other approach (Biot, 1941)> "bh.® flow of fluid 
and the elastic displacements of the soil are treated together 
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as a coupled problem, The governing differential equation 

is complex enough and can be solved in closed form only 

under a number of simplifying assumptions; else the help 

of numerical techniques such as Einite Element Method is 

to be sought. Once the soil is discretized and EEM is 

applied, the pore pressures and displacements in the soil 

mass are obtainable at any time. The discretization in 

time domain involves the choice of a weightage coefficient 

which decides the variation of pore pressures between two 

time steps. Due to the discretization of soil mass and 

the necessity to modify and solve a large matrix, high 

computational costs are involved in this procedure. However 

for a study of pore pressures and three-dimensional displacements 

in a specific problem, the approach with PSM solution is 

the only choice. The manifestation of Mandel-Cryer effect 

is possible due to the coupling of flow and elastic deformation 

problems in the Biot's approach. This effect is not 

indicated in the Terzaghi-Rendulic approach in which the I 

coupling is not there. However the difference between the 

two theories is never very large as far as the overall 

settlement response is concerned (Murray, 1978). i 

For the consolidation under a completely flexible 
circular raft quite a few solutions are available based on i 

the sinqple diffusion theory (e.g. G-ibson and Iiumb ( 1953), 



Davis and Poulos (1972)) as well as those based on Biof s 
theory (e.g. Be Jong (1957), Mandel (1957, 1961), Paria 
(1957), McWamee and Gibson (i960), Gibson et al.(l970)). 
Madhav and Vitkar’s (1980) solution takes into account 'the 
effects of construction time and anisotropy. However, 
no solution is found presently for consolidation under 
annular footing considering construction time and anisotropy 
with respect to permeability. 

TSnough evidence already exist (Davis, 1969) to show 
that there is good overall agreement between the results 
based on Biot's theory and Terzaghi-Rendulic theory in case 
of flexible circular (or strip) footings, 

7.1.3 Consolidation Under Completely Rigid Annular Pootings 

The problem of consolidation under rigid footing 
becomes much more involved because of the conditional rigid 
settlement under the footing, 

Analjrtical solutions for consolidation under rigid 
circular raft have been obtained by Agbezuge and Deresiewicz 
(1975) and Chiarella and Booker (1975). However these are 
for the condition of send. -infinite soil. No solution is 
found presently for the consolidation of finite layer under 
rigid annular (or even circular) footing taking into 
consideration construction time and anisotropy with respect 
to permeability. 





The only available approach for this type of problem 
appears to be the PEM based on Biot’s theory. Although it 
involves relatively higher computational cost, the contact 
pressure distribution, pore pressure distribution,deformation 
and stresses in the soil mass etc, can all be obtained. by 
advancing the solution with respect to time. However, it is 
possible (as shown in the present work) to use the simpler 
Tersaghi-Rendulic approach and finite difference scheme 
of solution to tackle this type of problem, provided one 
is not interested in the spatial distribution of stresses 
and strains in the soil. The overall behaviour of consolidation 
settlement and the changes in contact pressure distribution 
which are the important aspects for the purpose of design 
can still be investigated. The limitations of the Terzaghi- 
Rendulic approach are noted. 

The major considerations involved in dealing vn.th 
consolidation under rigid footings are the requirement of 
uniform settlement and the consequent changes in contact 
pressure as consolidation proceeds and the changes in 
pore pressures which, in turn, are caused due to changes 
in contact pressures, in addition to those due to the process 
of consolidation, Ydiile these considerations are automatically 
handled in Biot's theory due to the coupling of flow and 
elastic deformation problems, they are to be indirectly 
introduced if the Terzaghi-Rendulic approach is used. 
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Firstly it inust be recognized that the changes in 
contact pressure are mainly, if not entirely, due to the 
requirement of uniform settlement at all the times, under 
the rigid footing. This involves the response of the 
consolidating soil mass as a whole * The consolidation has 
to continue conditionally i.e, the Finite Difference scheme 
has to become implicit. Secondly, the changes in .contact 
pressure have to introduce corresponding changes in pore 
pressures. It is found (D’Appolonia et al,,l971) that 
the pore pressures axe predicted well as equal to octahedral 
total stresses cjilculated from the theory of elasticity. 
The calculation of changes in pore pressure due to changes 
in contact pressure at each time step would therefore 
involve elastic solution considering the elastic parameters 
of the soil and the geqmetry of the problem. The changes in 
pore pressure have beep, calculated in the present work on 
this basis and therefor? a coupling, although indirect, 
has been achieved between the flow and the elastic deformation 
problems. Thus the cha^nges in pore pressures due to 
consolidation, the changes in contact pressure and the 
consequent changes in octahedral total stresses in the soil 
mass and changes in pore pressiires due to changes in contact 
pressure get implicitly coupled. Although Manoel .Gryer effec-t 

total stresses do take 
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place corresponding to changes in contact pressures. Due to all 
this f this approach, with the indirect coupling employed 
in the present work, is presumed to predict the behaviour 
of the rigid footing fairly accurately, as far as settlement 
and changes in contact pressure are concerned, 

7,2 CONSOLIDATION SETTLEIi[El\fT 0? FLEXIBLE ANNULAR FOOTING 

7.2.1 Assumptions and Statement of the Problem 

7 . 2 . 1 . 1 Assumptions 

(1) The consolidating layer of soil is homogeneous and 
of uniform thickness, 

(2) The consolidation is governed by Terzaghi-Rendulic 
theory. 

(3) ' The loading can be instantaneous or time-dependent, 

and in case of time-dependent loading, its variation 
is linear with time. 

(4) The soil is anisotropic only with respect to 
permeability, 

(5 ) The loading is uniform giving rise to axissnmaetric 
problem. 

7. 2. 1.2 Statement of the problem 

A uniformly distributed load in the form of an annular 
ring of outer radius 'a' and inner radius 'na' is applied 
on a homogeneous layer of thickness *h* underlain by a rigid 

base (Fig, 7. 1(a)). 
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The load is applied in a linearly increasing manner 
(lig. 7 » 2 (b )) during ihe construction period t^> which in 
non-dimensional form is expressed as 
c t 

(7.1) 

where is the coefficient of consolidation of the soil, 

T ^2 = 0 refers to the case of instantaneous loading. 


The following types of drainage conditions at the 
boundaries are considered; 

(a) IFIB - Impervious footing and base 

(b) IPPB - Impervious footing, pervious base 

(c) PTIB - Pervious top, impervious base 

(d ) PTPB - Pervious top and base. 



where r indicates radial (horizontal) direction, z indicates 
vertical direction, and k and c^ are the coefficient of 
permeability and consolidation respectively. 


7.2,2 Analysis 

For a homogeneous isotropic elastic soil, the 
general differential equation for 3 -B consolidation is 




5u 


(7.5) 


2 ? 2 

a U a U i ae 

— C~ ( k + ~"W + ~ p ) ““ 

^ ax^ ay2 dz^ 3 at 

where, 6 is the total bulk stress (=a^+0 + a ), 

0 ^ is the three dimensional coefficient of consolidation 
and u is the excess pore pressure. 


The Terzaghi-Rendulic theory amounts to assuming 


^ - 0, which is approximately valid for all practical 

at 

purposes and is fairly accurate as far as overall consoli- 
dation process is concerned, (Davis, 1969), With this 
assumption, Eq. 7.3 written for axisymmetric case reduces to 


1 1 du a^u 


(7.4) 


The average degree of pore pressure distribution, Up, 
on any selected vertical line is obtained, after solving 
Eq. 7.4, from the expression, 

/u(z,o)dz-/u(z,t)dz 

U = — — ^ (7.5) 

^ /u(z,o)dz 

where t refers to the instant under consideration and zero 
refers to reference time. The average degree of settlement, 

Ug, Is approximately equal to Up for 3-D case if the time 
scale is adjusted in such a way that the time factor is 

defined as (Poulos and Davis, 1972; Christian- , 

19T2). 

(7.6) 

T - — 

^ ■ ■ ■ 

Where ci is one dimensional coefficient of consolidation. 



equation for anisotropic 


The governing differential 
ax i symmetric case is, 

^ 1 1 8u ,2^ 

3t r2'0r2 *r "a; '> * 


In the present study, it is necessary to know the 
initial excess pore pressures induced in the soil mass under 
a given loading on the surface. These can be obtained from 
the Henkel’s empirical relation (D’Appolonia et al.(l971). 


+ f T 


(7.8) 


where f = a coefficient depending on the conpressibility 
of soil. 


A study by the author (Table 7.1) for consolidation 
under a strip footing showed that the effect of taking 
the term f into consideration is quite less and it 

reduces as consolidation advances. If consolidation of 
two closely spaced strip footings is considered, a situation 
somewhat similar to ring footing, the effect of f is 

slightly more but still quite small. 

TABhE 7.1 : EFFECT OF ON GOHSOLIMTIOH UHDER 

FIBXIBIB STRIP IiOADING 

r "TT-- --t— - 1 --— — * — *1 ' 

b/h gA Percentage degree of consolid 
tion 

f=0 . f=0 f=.2 



In the evaluation of 


pore pressures beneath an 


embankment, D'Appolonia et al. (1971) also observed that the 
pure pressures are predicted well by the relation, based 
on three dimensional theory of elasticity, between pore 
pressures and cfoct Therefore, the initial pore 

pressures were obtained from the relation 

0 


where 6 is the total bulk stress. The values of bulk 
stresses are easily obtained by numerical integration 
of available solutions fix)m theory of elasticity (Poulos 
and Davis, 1974) • Thus pore pressure (or the change in 
it) at any point in the soil mass can be expressed as the 
product of a constant coefficient and the load (or the 
load increment). 


7.2.3 Solution Procedure 

For solving the governing differential equation (7.7) > 
it is written, using the usual Taylor's expansions and the 
approximations, in the Finite Difference form as 

6z 2 

^i,D,t+5t " ^ ^ ^ ^^i,j+1,t'*’ 

6z 





^^i+1 ,it'‘' ^i- 




2u, 


( 7 . 10 ) 


where , p 


(7.11) 


c^6t 

(62)2 

On the axis of symmetry (r=0), Sq, 7.7 in 

Finite Difference form becomes, 

6 s 2 


u 


+ 


"VZ 




(7.12) 


The time and time steps are given by 


- ^vz"^ Ovz^"*^ <2vz ^ \ 

T = -1^ ; 5.1 = ^ p ^ _ZZ 

a 3. ^vr 


The soil mass is represented by a grid (Fig. 7. 1(a)) 
made up of square units (6r = 6z), A lateral extent to 
represent infinite soil mass with sufficient accuracy 
is chosen and the consolidation is marched forward through 
one time step using Fq, 7.12 and Fq. 7.13 and the given 
drainage conditions on the boundaries. The average load 
increment during the next time step is now applied so that 
additional pore pressures will be generated, which can be 
expressed as, 




( 7 . 14 ) 


where b. ^ indicate coefficients of pore pressure generation 

XfJ 

and 6 p is the change in the imposed load ( 6 p = 0 for 
instantaneous loading and bj_^j depend only on the geometry 


of the problem). The next time step is taken up after 
adding the pore pressures so generated to the existing 
pore pressures and the process is continued until consoli- 
dation is almost complete. 

The finite difference grid and the time steps 

are so chosen that numerical stability and accuracy are 

achieved at reasonable computational cost. A value of 
2 

p = Cyp 6 t/( 62 ) = ,16 was used and found satisfactory, 

A lateral extent of 5 times the outer radius of the ring 
was found to be sufficient to represent soil of infinite 
lateral extent. 

7 , 2.4 Results and Discussion 

7 , 2. 4.1 General 

Although the present work is oriented to solution 
of problems involving annular footings, some results were 
taken for circular footing and they were found to tall 3 '‘ 
well with those of Pou^os (1967). (Typical curves are 
shown in Pig, 7.2). 

The settlement under the flexible footing is not 
uniform and the degree of settlement in all the ca-se.s herein 
refers to the settlement below the annular footing. 

The effects of various parameters are discussed 


below in brief. 



7. 2, 4 . 2 Effect of drainage boundary conditions 


Figure 7.3 shows for = 0, n=.75 and h/a =4 
the effect of drainage boundary conditions on the rate 
of. settlement. It is seen that the drainage boundary 
conditions have a significant effect on the rate of 
settlement (indicated by the degree of settlement U„ 
or degree of pore pressure dissipation U ). IFIB and 

Jr 

PTPB represent, as expected, the two extreme cases. For 

the particular case shown in Fig. 7,3, the degree of 

settlement, expressed as percentage, for PTPB case is 

2 

about 29, 61 and 96 percent for T = c-it/a =, 02 , ,2 and 
2 respectively, whereas the corresponding percentages for 
IFPB case are 6, 49 and 87 percent. The curves for IFPB 
and PTIB are found to lie between the two extreme cases, 

A similar trend exists for all values of T^, n, h/a but 
the difference in rate of settlement due to drainage 
conditions is less for higher values of T^, n and h/a. 

IFPB case is similar to IFIB case in the early 
time range but approaches PTPB case towards the end of 
consolidation, PTIB is similar to PTPB at small times but 
approaches IBIB towards the end of consolidation, i.e, in 
the beginning of consolidation the permeability condition 
of the footing or the top surface dominates and towards 
the end of consolidation the permeability of the base 



cases are intermediate 


dominates. Since the IFPB and PTIB 
between the extreme IPIB and PTPB cases, the latter are 
hereafter taken up for discussion unless otherwise mentioned, 

7. ^•4. 3 Effect of radius ratio n 

figures 7.4(a-) and 7 . 4 (b) bring out the effect of n 
on time rate of consolidation settlement, Prom these f’igures 
it can be seen that a small opening at the centre of a 
circular footing (n less than about , 25 ) has negligible 
effect on the rate of consolidation. As n increases, the 
consolidation tends to be faster, more so in case of IPIB 
case than PTPB case. An n approaches a large value closer 
to 1 , the ring acts like a narrow strip footing and the 
consolidation is very rapid. For example, the degree of 
consolidation expressed as percentage, achieved at T = 0.1 
is about 41 . 5 , 51 sm.d 66 percent for n=,5, .75 and .9 
respectively for PTPB case (h/a?:4> ^^=0), The corresponding 
Figures for IPIB case are 21, 54 and 58.5 percent. 

The final maximum elastic settlement of a flexible 
ring footing with =, 9 , h/a =4 (obtained by the author by 
numerically integrating tabulated values for point load 
(Poulos and Davis, 1974) is of the order of 27 percent of 
that for an equivalent circular footing (both the footings 
being of equal area and carrying equal load), Por IPIB 
case, and n = .9 , Ug from Pig. 7 . 4 (a) is about .13 and .585 


for T = 0.01 and 0.1 respectively. !Phe corresponding 
figures for n = 0 are ,005 and .12 . Hence the 
settlements at T = .01 ,0.1 and for n = .9 will be 
respectively about 7 , 1.3 and .2? times that for an 
equivalent circular footing. This shows the following 
two consequences which need attention from a practical 
point of view, 

(i) ring footings, the consolidation settlement 
in the early stages is going to be faster in 
comparison to an equivalent circular footing. 

(ii) The settlement of a ring footing during the later 
major span of time will be much less in comparison 
to that of an equivalent circular footing, 

7,2, 4.4 Effect of Period of construction 

Figure 7.5 shows the time-degree of settlement curves 
for h/a = 4 and 2, IF IB and PTPB conditions, n=,25, .5 , 

.75 and. 9 and = .1, .5 and 2. The effect of thne 

of construction on the nature of the time-settlement curves 
is appreciable, A higher value of T^j which means slower 
rate of construction induces settlement at a slower rate. 
However, during the long time of construction a considerable 
amount of settlement takes place. This aspect has been 
brought out clearly in Fig. 7.6 in which the degree of 
consolidation achieved at the end of construction has been 



237 


plotted versus From Fig. 7 . 6 , for PIPE and ii/a=2 

it cun be seen tiiat if ttie construction time = c^t^/a^ 
is equal to .1 , about 60 percent consolidation would be 
over at tb.e end of construction of a ring footing witb. 
n = ,9. Some of the values are given in Table 7.2 for the 
purpose of comparison at a glance, 

TABLE 7.2 : DEGREE OF OOESOLIDATION ACHIEVED AT THE EHD 
OF GONSTRUGTIOH 


h/a 

Drainage 

Condition 

T 

c 

Percentage Degree of 
Consolidation at T,^ 




n=.25 

n=.5 

n=,9 

4 

IFIB 

0.1 

5 

10 

38 



1.0 

46 

51 

72 

2 

IF IB 

0.1 

6 

12 

45 



1.0 

55 

57 

76 

2 

PTPB 

0.1 

55 

57 

60 ' 



1.0 

80 

81 

89 


From Fig, 7. 6 and Table 7.2 it is seen that the 
consolidation achieved during construction time is more 
for higher vaAues of n and lower values of h/a* 

Considering once again. Fig. , 7.5, it is seen that 
for higher values of Tq although consolidation is slower 





initially, it becomes rapid some time before the end of 

cent ruction (depending on Tq) and continues to show a 

soiTiowliE't rdt© of consolidation, For Bxample, for li/ap=4, 

H— ,5? I1‘ IB , tli 0 consolidation is advancod from 20 parcont 

to 80 percent in time T of about 1,1 and 1.4 for 1=0 and 

c 

Tp = 0,5 respectively. From Fig. 7.5 it can also be seen 
that the difference in time-settlement response in the 
extreme drainage cases IFIB and PTPB is lower for higher 
values of n, which means that narrow ring footings are 
less sensitive to drainage conditions, more so for higher 
values of h/a and higher values of In fact for very 

high values of (greater than 5 for isotropic case), 
the curves for all values of h/a , n and all types of 
drainage conditions appear to merge into a single curve. 

The charts as in Fig. 7.5 can be directly used 
to estimate the degree of settlement at any time for 
various values of depth of layer, outer and imier radii 
of loaded area and construction period T^. xilso charts, 
of Fig. 7.6 can be used directly to predict the degree of 
settlement at the end of construction. 

7. 2. 4. 5 Effect of depth of consolidating layer 

Figure 7.7 briiags out the effect of depth ratio 
h/a on the rate of consolidation settlement, '^Then the 
depth of the consolidating layer is very small in comparison 
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to the dimensions of the loaded area, the consolidation 
except for IB IB case, is nearly one dimensional and the time 
required to attain a particular degree of consolidation is 
nearly in proportion to the square of the depth. However 
as the depth ratio h/a increase, the corresponding increase in 
time for attaining a particular degree of consolidation is 
lesser because of the three dimensional effect. For example, 
referring to Fig. 7.7, it can be seen that for n=.5 and 
PIPE condition, the time for 80 percent consolidation with 
h/a = 4 is about 2.58 times and for h/a = 10 about 6.45 
times that for h/a = 2; whereas, these figures for simple 
one dimensional consolidation would have been 4 times and 
25 times respectively. The effect of h/a on rate of settlement 
is slightly lesser for higher values of n i.e. ring footings 
are less sensitive to h/a ratio than circular ones, 

7, 2, 4. 6 Effect of Anisotropy 

Phe effect of anisotropy with respect to permeability, 
as defined in Bq. 7.2 , has been studied for instantaneous as 
well as linearly built up load. Figs. 7.8 shows typical 
curves for Uo versus log T for m = 1 (isotropic case), 10 
and 100, The time-settlement characteristics change 
significantly as m increases from 1 to 10, this change being 
more, for higher values or n. For values of m beyond about 
10, the curves tend to become somewhat parallel to each other 



240 


for Tq 0. However, it has not been possible to suggest 
a simple relationship to correlate the curves for vacious 
values of m even beyond m = 10, 

Figure 7,8 also shows the Ug— log T cujrves for m =1, 

10, 100 and for = 0, and .1, It can be seen that the 
curves for various values of m tend to come closer to each 
other as increases. 

The curves such as in Pig. 7.8 can be directly 
used to estimate degree of settlement at any time for 
given values of h/a, n, and m. 

7.3 OOKOLIDATIOI SEmEMEOT OP RIGID AOTULAR POOIIIG 

f ^ 

7*3.1 Assumptions and Statement of the Problem 
7. 3. 1.1 Assumptions 

The assumptions contained in Section 7.2.1, 1 except the 
last one about the loading are valid. In this case the 
assumption in connection with the loading (replacing (5) in 
7.2. 1,1) is as below. 

The loading is in the form of a central load on 
a smooth rigid annular footing, giving rise to an axsymmetric 
problem, 

7. 3.1'. 2 Statement of the problem 

It is the same as in Section 7.3.2, 1 except that 
the footing is a smooth rigid footing carryiiig an axxsjrmmetric 
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loading of any nature (Figs. 7,3). 

7 . 3.2 Analysis 

The governing differential equation for the 
consolidation process and the method of calculating pore 
pressures from elastic theory remain the same as for 
flexible footing and hence the analysis for the flexible 
footing remains valid. However two major factors required 
to be considered for the consolidation under rigid footing 
are : 

( 1 ) The pore pressures generated due to external loading 
require revision at each time step since contact 
pressures under rigid footing keep on changing with 
t irae . 

(2) The consolidation process can not be explicitly 
marched forward in time since it has to incorporate 
the condition of uniform settlement under the 
rigid footing. 

These major factors are incorporated in the solution 
procedure, 

7.3.3 Solution Procedure 

The consolidation under rigid footing has to 

continue in such a way that the settlement under the 
footing has to be uniform. To achieve this the contact 
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pressures must vary satisfying the following conditions: 

( ) Redisijribution of contact pressures at eacli time 
step sliould take place so as to obtain a uniform 
change in degree of consolidation on vertical lines 
below the rigid footing under the combined effect of 
load increments (if any) and consolidation. 

(2) The overall chcuige in contact pressure must be 
equal to the change in external loading during 
a time step. 

( 3 ) The changes in contact pressure added up during 
the forward marching in time must yield the final 

contact pressure consistant with that obtained 
from elastic theory alone. 

In crrder to couple the rigid-settlement of the rigid 
footing with the consolidation, the ring is divided into 
small annular rings (say, *i' in number) and on each 
small ring, the contact pressure is assumed to be uniform 
(Pig. 7 . 9 ). Let 6pjj. indicate the average change in contact 
pressure for the k"^^ small ring, over a time step 6T, Let 
bi^ . . indicate the coefficient of pore pressure generation 

such that at any nodal point (i,3) soil, the pore 

pressure generated due to incremental contact pressures 

on all the small rings is given by . 

1 


(7.15) 



243 


Let t indicate the instant just after the 
instajitaneous pore pressure generation and just at the 
b ©ginning o.f b * coiisolida.tiionL Th© por© prGssur© ai 

t is given by 

^i,j,t"^ = ’^i,j,t ^^i,j (7.16) 


Since changes in contact pressures and in pore 
pressures are introduced due to the rigid-set tiement 
phenomenon in between the consolidation process the 
consolidation Eqs. 7.10 and 7.12 in finite difference form, 
are firstly required to be modified using Eq. 7.16 above 
to get pore prc 3 ssures at t"*" + 6t instead of at t + 6t. 

Bqs. 7.10 and 7.12 after rearranging, will therefore be of 
the following form: 


j,t'' + 6t 


®i 


+ 0, 


1 

s 

k^1 


6pic 


(7.17) 


where corresponding to Eq. 7.10, 
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6z 2 6z 

Kp = p (~) _ p 

2j,6r 


K, 
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Similarly, values of and can be written corresponding 
to Eq, 7.12, 


Bie common change of degree of consolidation 

(assumed equal to chcuige of degree of settlement) along 

any vertical line below the rigid footing is given by 

/u,- (3,o)dz - /u^Czjt'*') dz , 

6U = i (7.18) 

f'u.j_(z,o) dz 

where Uj^ indicates pore pressures on i"^^ vertical line. 

The integrations are carried out numerically so 
that if 's* is the number of nodes on each vertical line, 
and if indicates the coefficient of numerical integration 

V 

■feh 

for 3 node on the vertical line, 



/u^(z,o) 

CO w ^ 

It 

K.j.o ^ S}'> 

(7. 

.19) 

and 

/ujL(z,t'*') 

s 

= 1 
3='' 


(7. 

.20) 


In case of time dependent loading, the zero in Bq. 7.19 
and in the numerator of 'Bq. 7.18 is to be replaced by the time 
at the beginning of the time step. 



Indicating x by and combining Bqs. 7.15 

■throufrh -fici. 7 .20, the following form of equation is obtained 

°i 1,1 ®I>k - V ® (7.21) 

f 

Where are known quantities in terms of coefficients of 
pore-pressure generation and coefficients of numerical 
integration and are known quantities in terms of pore 
pressures at the beginning of a time step. 


The above equations are written for '1* number of 
selected vertical lines below the rigid footing. One more 
equation is written for overall equilibrium of the footing 


1 

Z Aj_ 6p]- 
*1 ** 


6W 


( 7 . 22 ) 


where Aj^ is the area of k small ring, 6i|j. is the 
corresponding^ change in contact pressure and 6W is the 
change (if any) in the external load during a time-step 
under consideration, 

Eqs. 7.21 and 7,22 are combined and written in 
the matrix form as, 


■11 


h2 


<^21 ^22 


o1 

g2 


« • '♦ • 


• • * • * 


• • • • • 


<^11 


L‘^1 


^^12 * 

A2 . 


U, 


ol 


0 




, 

"ei^ 


6p2 


®2 

1 

■ 

• • • 

• • • 

1 = < 

' 

• •' « 

i . 


6pi 

I 

i 

1 

1 

l®i 

1 

—• f 

6U 


i 6ff V 


(7.23) 



In the abOYe Eq. 7.23 , the matrix coefficients 
are derived from 0; and are therefore calculated only once. 
The K.11,.,U. vector is derived from known pore -pressures at 
the be/^innij;!^.^ of a time step and hence E(]_, 7 •25 is to be 
solved at each time step to get values of changes in 
contact pressure and change in degree of consolidation 
(settlement) of the rigid footing. The incremental pore 
pressures axe evaluated using Eq. 7,15 and are added, 
after applyin^r a vreightage coefficient a as required in 
coupled pi'oblems involving time (Besai and Ohristian, 1977) 
to the original pore pressures and consolidation is advanced 
through a time-step. The process is continued until 
consolidation is almost complete, Eig. 7.10 shows the 
flow -chart for the computations involved. 

7.3,4 Results and Discussion 

7. 3. 4. "I General 

The computer programme developed for this Chapter of 
the thesis is a general one covering flexible as well as 
rigid footings and for the latter, provision is made to 
generate a matrix, triangularize it and use it at the 
beginning of each time step. The size of the matrix is 
quite small being equal to the number of rings into which 
the footing is subdivided plus one (Eq. 7.23), 
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The soil is discretized as for flexible footing 
and as mentioned in Section 7. 2.4.1 (]?ig. 7.9)^ 

Ihe increase in contact pressures near the edge 
is critical since it may introduce increase in bending 
moment hence the study and discussion is oriented to 
signify this aspect. The increase is found to be maximum 
for XT IB case of drainage and for intantaneous loading 
^ discussion refers to these cases unless 

otherwise stated. 

The discussion made in Sections 7.2. 4.2 to 7. 2. 4. 5 
applies broadly to rigid footings also. 

7 . 3 . 4 .‘- Ghfmges in contact pressure distribution 

The ratio of contact pressure p to average contact 
pressure (total load divided by area of footing) is 
plotted versus dimensionless time factor T in SUgs, 7 J 1 to 
7.14. In Fig. 7. 'll, the contact pressure at a value of 
T = ,0225 is shown for a ring footing (h/sF: 2 ) and at 
T -- ,08 for a circular footing (h/a = 2 ), It these time 
factors, the contact pressure shows an increase at the 
outer eind inner edges of the ring footing and at the outer 
edge of the circular footing. In both the cases, the 
contact pressure away from the edge /edges has decreased. 

This change in the distribution is unfavourable from the 
point of view of moments induced in the footing, Figs, 7 . 12 , 7 .l 3 
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and 7.14 show the peaks of r‘n■n+af^+ 

I contact pressure aweraged over 

the one litth annular portion near the outer edge of the 

footing, in Fig. 7.12(a) the effect of h/a is brought out. 

Xt is ouon t.}iat the inaximuni (peak) values of contact 

pressure are higher for larger values of h/a, the values 

of p/po for n = 0 , =0 being about 1 . 9 , I .84 and 1.61 

for h/a = 8, 4 and 1 respectively. They work out to an 

increase of about I 4 , 12 and 9.5 percent over the initial 

contact pressures (which are already high near the outer 

edge), it is also seen that the peaks are attained at 

about 0.’ = 0.08 for the cases considered. 

'Dho peak p/Pq values obtained from Fig. 7.12(b) 
for h/a =.^1, n=0, = 0 are about 1.61 , 1.59, 1 .58 and 1.56 

for IFIB, X’TIB, IFPB and PTPB boundary drainage conditions 
respectively. This shows that IPIB is the most critical 
condition, 

Pifpire 7.15 shows the changes in contact pressures 
for Tq = 0, n = 0, IFIB, for the central ring and for an 
annular portion near the half radius of the circle. The 
average over one fifth width is plotted in this figure also. 
Whei‘oas the edge contact pressures increase, attain a peak 
and decrease with time, the contact pressures away from the 
outer edge decrease, attain a minimum and then increase. 
This happens because the total load on the footii:^ remains 
constant. It is seen from this figure that the minima do 
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not ocour simultaneously with the peaks of outer edge 
pressures. For example, for hA =2, the peak of outer edge 
prosniro occurs at about 1 = o .08 (Fig. 7.12) v^hereas 
the ininiiiiuni values at the central ring and at the half radius 
reach et about 1 =0.2 and T = 0,1 respectively, This means 
that the most unfavourable contact pressure distribution 


may or may not coincide with the peaking of outer edge 
contact pccsrnuru. Prom the same figure (Pig. 7.13) it is 
soun tliat the contact pressures near the centre show a 
maxirmuti decruaou and that too for higher values of h/a. 
For uxample, for h/a = 8, the contact pressure near the 
centre drops from about .5 to .39 times average pressure 
on the footing. 


Figure 7.14 shows the effect of radius ratio n on 
the changcas in edge contact pressures. The peak p/Pq values 
are about 1.915, 1 .856 and 1 .796 for n=.5 , .2 and 0 
respectively, ThfeSQ work out to an increase 01 about 13.3, 
11,S , imd 10.9 percent over initial values of edge contact 
prciscure for n = .5, .2 and 0 respectively. This shows that 
the increase is more for higher values of n. However for 
a ring footing;? the bending moments are already less than 
those for a circulai' footing. Hence the change in contact 
pressure will be more critical for circular footings than 
for ring footings. 





Bor an actual footing, the increase in contact 
prcrv'ure will not be so much as indicated in the preceeding 
par.-vuvaph!.; because of several reasons: 

(i) llie lootiiig will have finite rigidity and hence 
the valiie of increase in edge contact pressure wi ll 
lie inbetween that of perfectly rigid footing as 
analysed in this chapter and zero for fully flexible 
foot Ing. (Also the initial edge contact p^ressure itself 
will be lower for flexible footings than for rigid 
footiru^s). 

(ii) 'Jhc peal, construction time being greater than zero, 
the peak values will drop (This aspect is discussed 
later). 

Fipiire 7-15 shows the comparison between contact 
pressure changes for strip footing obtained from the 
present analysis applied to strip footing and from Biot's 
theory usin'; the computer programme given by Verruijt (1977). 
The agreement is fairly good. 

7, 3,4.3 ii;:ffect of infinite rigidity in comparison to 

infinite flexibility on time rate of settlement 

The present study shows that the consolidation 
..nder rigid ring footirg with UIB drainage condition would 
be faister than that under a unifomly loaded flexible footing 



for the sane drainage condition except tov^ards later stages 
of CO ,..'0 Li.d-itioii vvhen it becomes slightly slovier (curves for 
11 O)* I- n .5 the difference is not much, fhe reason 

for faster consolidation under rigid ring footings (circular 
as a pax'ticular case) is that the high edge pressures, 
which also increase initially, increase the flow near the 
edges. The preferential outward radial flow under an 
axisymmetrin loadinig is further accelerated due to high 
edge pri.'ssuro at the outer edge. 

In the cases shown in Figs. 7.16 (a) and 7, 16(b) 
the dirrorenco in the degree of settlement for rigid and 
floxiblu footings for the case of = 0 is appreciable 
only in the earlier stages ( T less than about .4), and 
this difreronco is less for higher values of h/a. In 
Fig. 7.16(b) the curves for = .5 are also shown from 
which it is aeon that for this the difference in the 
time rate .of settlements for rigid and flexible footings 
is quite small as compared to that for Tq = 0. At hign 
values of (> 2) the behaviour of flexible and rigid 
footijign with respect to time rate of settleiii.ent ,is almost 
identical, 

A coinijarison with available result for a rigid 
circular footing on serai-iafinite soil (Agbezuge and 
Bore siewic z, 1 975 ) with the result from the present study 
for deep layer of soil is shown in Fig. 7.2. The agreement is 

sec'n to be good. 
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7 . 3 . 4.4 Effect of radius ratio n 

C/O ribolidat ion. is faster for higher values of n in. 
u way vB.cy biiiiileir to that for flexible footing (discussed 
in Sect.ion f.2,4.3). Ihe peaks in contact pressure show 
a higher percentage increase over initial values with 
higher values of n (Pig. 7. "^4 and discussion in Section 
7 *2, 4.2), Pig. 7.17 also brin^out this aspect, Por 

C) , the percentage increase in edge contact pressure 
(taken over one fifth width of ring) is 115.0, 111.9 and 
110.9 for n .5 , .2 and 0 respectively. Therefore the 
effect of increase in edge pressures on bending moments in 
the footing will be slightly more for higher values of n. 
However, as mentioned earlier, the moments are already 
lower for higher values of n. (For example for relative 
rigidity of 10, the nondiraensional design moments (M^/q.ya ) 
for uniformly loaded footings for h/a = 2 are about . 0205 , 
.029 and .032 for n = .5 , .2 and 0 respectively (Chapter 4). 

7 . 3 . 4. 5 Dependence of peak contact pressure on construction 

time Tq 

Figure 7,17 shows the effect of on peak values. 
The figures in parantheses express the increase, in 
percentage, of peak contact pressure value, averaged over 
one fifth width of footing, over the initial contact 
pressure. The increase is less for higher values of T^, 
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as expect<.Kl, but the increase in peak contact pressure 
drops si./piiL’icjantly for higher values of n. For example, 
the increase drops from 110.9 to 109.8 for n = 0 and from 
113.0 to 108.7 for n = ,5 as increases from 0 to .2. 
liiince will be greater than zero in actual practice, 
the increase in contact pressures will not be as severe 
as for - 0 and much less for higher values of n, Phis 
again goes in favour of providing footings with higher 
values of n, 

7.3»4.6 ih'fect of anisotropy with respect to permeability 

Fii'hU’o 7.18 shows the effect of anisotropy on the 
time rate of aott lenient behaviour for a particular case 
of h/a -- 2, n =■ ,5 , IFIB. The curves for flexible footing 
are also shovm (dotted lines) for comparison. It is seen 
that consolidation is faster, for higher values of m 
( ra k^/kj, = Gvr/^vz^ expected. It is seen that the 
cui'ves for ra value more than 10 are somewhat parallel 
to each other. However it is not possible to suggest a 
simple correlation. The difference in behaviour of rigid 
and flexible footings shows the same trend as for isotropic 
case i.c,, the consolidation for rigid annular footing 
being faster than that for uniformly loaded flexible footingj 
in the earlier period of consolidation. 
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i'i^nire 7.1 9 shows the effect of anisotropy 
on perU-: \r, •flues of contact pressure. It is seen that 
the inot'easo in contact pressure at peaks (shown in 
paranthuBis) is reduced with higher values of m. However, 
the difrcjrenoe is snail and it could be said that anisotropy 
has little effect on peak contact pressures. 

7,4 COMCJLIJSIOHS 

(1) The consolidation settlement under flexible and 
rif^id aamular footings has been investigated 
usirit Ter iaaghi-Rendulic theory, For rigid footings, 
the flow due to consolidation is linlced up with 
elastic behaviour of soil mass through pore pressure 
generation and the conditional uniform settlement. 

( 2 ) Plots for time-settlement response for various 
parjamoters including construction time and degree 
of anisotropy are presented which could be used 
foi’ the design of annular footings. 

(3) The consolidation under a flexible annular footing 
is -Taster for higher values of n especially in the 
earlier stages of consolidation. 

(4) Considerable consolidation settlement takes place 
by the end of the construction period, especially 
for higher values of n and lower values of depth 
of consolidating layer. 



iii.ng foo’biti^'js ar6 less sensitive to h./a ratio than 
ciirtjular ones in regard to time rate of settlement, 
An:inoi,ropy with respect to permeability has a 
signi, Lcimt effect on the time rate of settlement, 

Ihe contact pressure changes during the consolidation 
settlement of a perfectly rigid annular (or circular) 
roatiyis!: are not insignificant, Por ring footings, 
the peak percentage increase in edge contact pressures 
is a.lightly higher than for circular footings. 

Thu ])eak contact pressures decrease considerably 
if time of construction is taken into account, 
(Jorisolldation under rigid annular (or circular) 
footini;a is faster than that under flexible footings 
in the initial stages ( T less than about ,4). 

Later the difference is quite less. 
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FIG.7-6 FLEXIBLE FOOTING-DEGREE OF 
CONSOLIDATION AT THE END OF 
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T=c,t/a2 

FIG-7-17 EFFECT OF CONSTRUCTION TIME ON 
PEAK PRESSURES 
(Figures in brackets indicate 7o peak 
edge pressure over initial averaged over 
one fifth of ring) 
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CHAPTER 8 


SUMIVIARY AM) GOICJIUSIOMS 


8.1 SUMimRY 


The aim of this thesis is to make a comprehensive 
study of the perfOCTance of annular footings as foundations, 
with reference to the load-settlement and time-settlement 
response. Although annular footings are highly suitable as 
foundation for structures like water tanks and towers and 
although millions of such structures are constructed every 
year all over the world, the potential of the use of annular 
footings does not seem to have been explored fully, probably 
because of lack of simpler means of predicting the behaviour 
at all loads levels and at all times. This thesis is an 
attempt to provide, for annular footings, a broad parametric 
study, and some design data for direct use, Ihile circular 
f ootingi is studied as a particular case of annular footing , 
a critical comparison between the two types of foundations is 
included with reference to all the aspects covered in the 
thesis, 

Except for extremely flexible footings, the contact 
pressures near the edges of an annular (or circular) footing 
are very high even under small loads, giving rise to plastic 
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yielding. Plastic yielding can not therefore be neglected 
in any rational analysis. So also the tension cut-off 
which occurs especially under eccentric loading, due to 
the inability of the soil to take tension. These aspects 
are given utmost importance in the thesis. 

The load-settlement response has been studied for the 
following cases of annular footings (circular footings as 
a particular case). 

(1) Rigid centrally loaded footings* 

(2) Axi symmetrically loaded footings of finite rigidity 
with typical loading configurations, 

(3) Rigid centrally loaded footings supplemented by 
incompressible piles (Piled Ring Poundations) , 

(4) Rigid footings eccentrically loaded. 

Except for case (5) above, the following t 3 ^es of problems 
are considered. 

(a) Footing resting on the surface of a semi -infinite 
soil. 

(b) Footing embedded in the semi -infinite soil. 

(c) Footing resting on the surface of a soil layer 
of finite depth underlain by rough rigid base, 

A piled annular footing is assumed to rest at the 
surface of a semi -infinite mass of soil. 



278 


In all the above cases, the footing and piles are 
discretized into small elements, the stresses on the elements 
are assumed to be uniform and idealised to act as point loads 
wherever appropriate at the centres of elements. The 
available elastic solutions (such as by Boussinesq, Mindlin) 
are employed to obtain the influence coefficients of the 
flexibility matrix connecting the stresses and displacements. 
Wherever stresses tend to exceed the yield stress value 
of the soil, they are held at the yield stress value 'and 
wherever stresses tend to become negative, they are held at 
zero value. The influence coefficient matrix is partitioned, 
the rows corresponding to yielded zone are removed using 
the limiting values as mentioned and it is modified to 
represent stress-displacement relations in elastic zone 
only. Solutions are obtained after incorporating overall 
equilibrium conditions of the footing. An entire load- 
settlement response, distribution of contact pressures, and 
the progress of yield zone under the footing are obtained as 
solution at all levels of loading and eccentricity if any. 

In case of eccentric loading, the settlement response is 
obtained in terms of central settlement and rotation about 
the centre of the footing. 

The effects of embedment and the effect of depth 
of finite layer are studied and compared with reference 
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to the condition of semi-infinite soil, in all the cases 
except piled footings. For piled footings, the effect of 
length to diameter ratio and number of piles is studied in 
combination with the footing size, ihe effect of positioning 
of piles below the footing in one or more concentric circles 
is also studied. 

The rate of consolidation settlement is studied for 
fully flexible and fully rigid ainnular footings resting on 
a saturated soil layer. Various drainage boundary conditions 
are included for the purpose of the study. The effect of 
period of construction, in which the load is assumed to 
build up linearly, is studied and the amount of degree of 
consolidation during the construction period is quantified 
for annular footings. Terzaghi-Rendulic theory is used with 
time scale correction to formulate a direct finite difference 
scheme. The same theory is used for rigid footings but 
using a scheme which becomes implicit because of conditional 
uniform settlement under the rigid footing and v^hich involves 
indirect coupling of flow-elastic deformation problems 
because of generation of pore pressures under changing contact 
pressures in terms of total stresses elastically induced in 
the soil mass. Changes of contact pressures with time under 
the rigid annular footing are also obtained. 
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In the study of consolidation settlement, the 
effect of anisotropy resulting from higher permeability 
in the horizontal direction in comparison to that in the 
vertical direction is also studied briefly, 

8.2 GOMJIUSIOIB 

(a) With the help of a simplified elasto-plastic analysis, 
it has been possible to obtain load-settlement 
response of annular footings over the entire range 

of loading upto failure. The load-settlement response 
when compared with some available results from 
experiments as well as those obtained from FEM, tally 
well. 

(b) The load-settlement response is non-linear throughout 
the range of loading. At lower load levels, it is 
nearly linear and as load approaches the failure 
load, it is highly non-linear and shows a rapid 
increase of settlement (and/or rotation in case of 
eccentric loading). 

(c) For the same load, annular footings indicate less 
settlement in comparison to circular footings, 

(d) In contrast to circular footing, an annular footing 
shows a nearly linear load settlement response upto a 
considerable load level, after which settlement 
increasely rapidly. 



The plastic yield of soil starts near the outer 
edge of a circular footing and advances towards the 
centre a,s load increases. In case of annular footing, 
it starts near the outer as well as inner edges and 
progresses towards the centre of the v^^idth of the 
footing, this being faster from the outer edge. As 
load approaches failure load, the contact pressure 
tends to become uniform. 

The settlement of embedded footing is less than that 
of a surface footing, but the nature of contact 
pressure distribution and the progress of yield are 
similar in the two cases. It has been possible 
to apply corrections to the settlement obtained 
for surface footing to get the settlement of embedded 
footing at any load level. 

The settlement of a footing resting on finite depth 
of soil is less than - that resting on semi-infinite soil 
but no direct correlation could be found since the 
nature of contact pressure distribution itself is 
different in the two cases. The performance in 
terms of settlement for the two cases is qualitatively 
similar. Annular footings are advantageous in this 
case also in comparison to circular footings from 
settlement point of view, at all load levels. 



For a centrally loaded footing of finite rigidity 
it is found that whereas the overall load settlement 
response is not very sensitive to relative rigidity K , 
the differential settlement, and the bending moments 
and shear in the footing are sensitive to K, For 
lower values of K (more flexible footings) the 
differential settlements are higher and the design 
moment and shear in the footing are lower. These 
three quantities and also settlements are lower for 
embedded footings as well as for surface footings on 
finite depth of soil in comparison to footings on 
semi-infinite soil. 

For footings of finite rigidity loaded uniformly, 
it is found that the footings should be designed for 
values of moment and shear much lov;er that those 
given by elastic theory. 

For footings of finite rigidity it is found that 
the worst position of load is near the centre or 
inner edge, from the points of view of settlements 
and design moment and shear in the footing. This 
type of loading may also introduce separation of 
footing from the soil near the outer edge, if the 
footing is relatively flexible. 



^or a loading in the form of an annular ring, it is 
found that the design moment and shear in a footing 
of finite rigidity could be calculated on the basis 
of uniform contact pressure equal to yield stress 
of soil, 

A ring (or circular) footing can be advantageously 
supplemented by providing a small number of short 
piles, due to which the settlement of the footing 
reduces considerably. If the piles do not dominate, 
the bearing capacity of the annular footing is 
utilised advantageously, A judicius combination 
of advantages of a raft and piles can yield a very 
economic alternative. If a single ring of piles is 
provided it should be positioned' in such a way that 
it divides the area approximately into equal portions. 
For ring footings it should be near the centre of 
the width. This is found to give minimum settlement. 
It may be found more convenient from the points of 
view of positioning of the loads on the footing and 
of structural design of the footing. If more than one 
row of piles is provided an angularly staggered 
arrangement is preferable and the response is better. 
For eccentrically loaded rigid annular footings, it is 
found that the rotations are also less for annular 



footings, carrying the same load, in comparison to 
circular footings. It is found that eccentricity 
also adds to yielding of soil on the compression 
side and. introduces separation or lift-off on the tension 
side. Separation can occur at eccentricities less than 
those worked out from elastic theory, because of 
plastic yielding of the soil, 

lerzaghi-Rendulic theory has been employed, along 
with time scale corrections to the axisynmetric 
problems df consolidation under flexible and rigid 
annular footings. It is found that consolidation 
settlement of annular footings is faster for higher 
values of n (narrower footings). The major portion 
of consolidation settlement in case of annular footings 
will occur in the earlier stages. 

It is found that considerable amount of consolidation 
settlement would occur during the period of construction 
itself. This aspect has been quantified for annular 
footings. 

It is found that anisotropy with respect to 
permeability has a marked effect on time rate of 
settlement of annular footings. 

For consolidation under rigid footings, Terzaghi-Rendulic 
theory only is used. However the finite difference 
scheme becJ^mie implicit because of conditional rigid 
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settlement. Elastic theory is used to obtain 
pore pressures due to changes in contact pressirre 
with time due to which an indirect coupling of 
flow and elastic deformation problems is achieved. 

The time rate of settlement and contact pressure 
changes are studied, 

(r) In aj.1 the cases, charts are presented which would 
directly serve as design aids for the range of 
parameters considered, 

(s) In all the problems, it is found that the computational 
effort and cost involved is considerably less although 
an entensive parametric study and generation of 
considerable amount of design data were involved, 

8,3 SUGaiJSTIOES FOR FURTHER STUDY 

(1) Application of the simplified elasto-plastic 

settlement analysis used in this thesis to the 
following cases: 

(a) Footings of any shape. 

(b) Surface footings of any shape resting on Gibson 
soil (using available elastic solutions (Brown 
and Gibson, 1979)), 

(c) The case of layered soils using simplified 
equivalent layer approach (Odemark, 1 949 ) or 
tabulated elastic solutions(e.g. compilations by 
Poulos and Davis, 1974). 
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( 2 ) 


( 3 ) 


( 4 ) 

( 5 ) 

( 6 ) 


( 7 ) 


(d) Single pile-square raft unit and pile-raft 
systems of usual configurations. 

Experimental verification through laboratory 
and field models for annular footings with and 
without piles. 

Using available elastic solutions for horizontal 
loading and including the effect of adhesion between 
the footing and the soil. 

Studying the effect of compressibility of piles 
and of adhesion coefficient a. 

Study of piled ring footings under eccentric loading. 
Use of the computer programme developed for this 
thesis for consolidation under multiple strip footings 
and eccentric loading. 

Study of consolidation under construction loading 
sequence other than linear. 




287 


RBJ’ERSaaES 


Agbezuge, L.K. and Deresiewicz, H* (1975)» Th.e Gonsolidation 
Settlement of a Oircular Eooting, Israel J, 'lech.,, 

Vol. 13, pp. 264-269. 

Ahivin, R,G. and Ulery, H.H. (1962). Tabulated Values for 

Determining the Complete Pattern of Stresses, Strains 
and Homogeneous Half Space, Highway Research Board, 
Bulletiin 342, pp, 1-13. 

Asproudas, S.A. and Desai, O.S. (1975). Analysis and 
Application of a Finite Element Procedure for 
Gonsolidation. Report, Deptt. Oiv. Engrg, , Virginia 
Polrfcech. Inst, and State Univ, , Blackburg, Virginia, 

24061 , May ’1975. 

Babu Shankar, H. (1977) . Stress Distribution Under Embedded 
Circular Footing, Indian Geotech. J. , Vol, 7 (1977), 

Ho , 2 , pp , 178-181 . 

Banerjee, P.K. (1971). Foundation with a Finite Elastic 

Layer. Oiv. Engrg, Publ. Y/ks, Rev., Vol. 66, pp, 1 1 97 - 1 202. 

Baneroee, A. and Jankov, Z.D. (1975). Gircular Mats Under 

arbitrary Loading. J. Struct. Div. Ptoc. ASGE, Vol. 101, Ho. 
ST 10, pp. 21-33-2145. 

Barden, L. (1966). Contact Pressures Under Gircular Slabs. 
Struct. Engrg. Vol. 43, PP. 153-154. 

Bell, R.A. and Iv/akiri, J. (1980), Settlement Goraparison 

Used in Tank Failure Study. J. Geo tech. Engrg. Div., 
Proc. ASGE, Vol. 106, Hb. 412, pp. 155-165. 

Biot M.A. (1941). General Theory of Three-Dimensional _ 

Gonsolidation, J. Appl. Phys. Vol. 12, lo,2, pp, 155-164. 

Bobe, R. and Piotsch, C. (1981). Settlement Calculation by 
a Hew Strength Theory, Proc, 10th Int. Conf, DfitPji, 
Stockholm, V 0 I.I, pp. 53-56. 

BoekMadsen,G. and Lagoni,Pl 1981 ) . Pile Foundation Problems in 
White Chalk. 10th Int. Gonf. SMFE, Stockliolm, V 0 I.I, 

pp. 51 - 34 . 



288 


Booker, J.R. (1973). a lilUmerical Method for the Solution 

of Biot'e Consolidation Iheory. Quart. J. Mech. Appl. 
Maths., Yol. 26, Bo. 4, pp. 457-470. 

Borowicka, H. (1936), Influence of Rigidity of a Circular 

Foundation Slab on the Distribution of Pressure Over 
the Contact Surface. Proc, 1st Int. Conf. SIvJPB, 

Yol. 2, pp. 144-149, 

Borowicka, H, (1930), Druckverteilung Unter Elastischen 

Flatten. Ingenieur Archiv, Yol. 10, Bo, 2, pp. 11 3-125. 

Borowicka, H. (1943). Uber ausmittig belaste starre platten 
au'^ elastich-isotropem Untergrund, Ingenieur -Archiv, 
Berlin, Yol.l, pp. 1-8. 

Boswell, L.F. and Scott, G.R. (1975). A Flexible Circular 

Plate on a Heterogeneous Elxstic Half Space, Influence 
Coefficients for Contact Stress and Settlements, 
Geotechnique , Yol. 25, Bo, 3, pp. 604-610, 

Bozozuk, M, ( 1972 ), Foundation Failure of the Yankleek Hill 
Tower Silo, Performance of Barth and. Earth Supported 
Structures, ASGE Speciality Conf., Yol.l, Pt,2, 
pp. 885-902. 

Boussinesq., M. J. (1885). Application Be s Potent iels, al' 

Etude de 1 'Equilibre et du Mowement Bes Solides 
Elastiques. Gauthier-Yillars, Paris. 

Brinch Hansen, J. (1970). A Revised and Extended Formula for 
Bearing Capacity. ]3ulletin Bo, 28, Banish Geotech, 

Inst., Copenhagen, pp. 5-11. 

Brown, K.F. , Poulos, H.G. and Wiesner, T.J. (1975). Piled 

Ra-f't Foundation Be siffn. Proc. Symp. Raft Foundations, 
CSIRO, Aust. 

Brown, P.T. (1968). The Effect of Local Bearing Failure on 
Behaviour of Rigid Circular R.8.fts. Ii^st. Bngrs., 
Austrrilia, Civil Engrg. 'Trans., Yol, 10, pp. 190-192. 

Brown, P.T. (1969), Numerical Analvses of Uniformly Loaded 

Circular Rafts bn Beep Elastic Foundations. Geotechnique, 
Yol. 19 , pp. 399-404. 



289 


Brown, P.T, (1965). Numerical Analyses of Uniforinly Loaded 
Gircular Rafts on Elastic Layers of Rinite Depth, 
Greotechaique, V0I.I9, pp. 301-306. 

Brown, P.T. (1969). Raft Foundations. P.G. Course on Analysis 
of the Settlement of Foundations. School of Civil 
Engrg., IJniv. of Sydney. 

Brown, P.T, and Gibson, R.S, (1979). Surface Settlement of a 
Finite Elastic Layer Vfeose Modulus Increases Linearly 
With Depth. Univ. Sydney, School of Civ, En^re:. Res. 

Rep. ifo. R 34-5, 1979. 

Buragohain, D. N, and Shah, V.L. (1977). Curved Interface 

Elements for Interaction Problems. Proc, Int. Symp. 

Soil Struct. Interaction, Roorkee, 1977, pp» 197-201. 

Burland, J.B,, Butler, F.G. and Dunican, P. (1966), The Behaviour 
and Design of Large Diameter Bored Piles in Stiff Clay. 
Proc. S37mp. on Large Bored Piles, Inst. Civ. Engrs., 
London, pp. 51-71. 

Butterfield, R. and Banerjee, P.K. (1971). A Rigid Disc 

Embedded in an Elastic Half -Space. Geotech. Engrg., 

V0I.2, Ho. 1, pp. 55-52. 

Carrier, W.D. and Christian, J.T. (1973). Rigid Circular 

Plate .Resting on a Non-homogeneous Elastic Half space. 
Geotebhnique, T0I.23, N0.I, pp.67-84. 

Chakravorty, A.K. and Ghosh, A. (l975). Finite Difference 

Solution for 0ij?cular Plates on Elastic Foundations. 

Int. J. jNumer. Mech. Engrg., Vol.9, PP. 73-84-. 

Ghan,3.C. and Cheung, Y.K. (1974). Contact Pressure of 

Rigid Footings on Elastic Foundations. Civ. Engrg., 

pp. 51-59. 

Chen, W.T. and Engel, P.a. ( 1972). Impact and Contact __^Stress 
Analysis in Multilayer Media. Int, J. Solids Struct., 
V0I.8, Ho. 11, pp. '1257-1281 . 

Cheung, Y.K. and Nag, D.K, (1968). Pl.ates and Beams on EPastic 
Foundations- Linear and Non-linear Behaviour. Geotech- 
nique, Vol. 18, pp. 250-260. 

Ohiarella, G. and Booker, J.R. (1975). The lime -Settlement 

Behaviour of a Rigid Die Resting on a Deep Clay Layer, 
wuart. J. Mech. Appl. Match., Vol. 28, No. 3, 
pp. 317-328. 



290 


Christian, J.'J?, and Boehmer, J.V/, (1970), Plane Strain 

Consolidation by Finite Elements. J. SivIPB DiY. , Proc. 

ASCE, VoI. 96, Mb. SM4, pp. 1435-H57 . 

Clough, R.W. and Woodward, R.J. (1967). Analysis of 

Embankment Stresses and Deformations . J, S3/IFE Div, 

ASCE, VoI.93, Mo. SM.4, pp. 529-549. 

D’Appolonia, D.J. , lambe, P.W. and Poulos, H.G-. (l97l). 

Evaluation of Pore Pressures Beneath an Embankment. J. 
SivIPEDiv., ASCE, Vol.97, pp. 881-897. 

Das, S.C. and G-angopadhyay, C.R. (1973), Undrained Stresses 
and Deformations Under Footings. J. G-eotech. Engrg. 

Div. Proc. ASCE, Vol. 104, Mo. GE1, pp. 11-25. 

Dave, P.C. (1977). Soil Pressure on Annular Foundation. J. 

Inst, Engrs. India, Vol. 58, Part CL 2+3, pp, 56-61. 

Davis, E.H. (1969). Ihe Analysis of the Settlement of 

Foundations, P.G, Course on Analysis of the Settlement 
of Foundations, School of Civil Engrg., Univ, of 
Sydney. 

Davis, E.H. and Poulos, H.G. (1972). Rate of Settlement 
Under Two- and Three-dimensional Conditions, 

Geotechnique, Vol. 22, Mo.1, pp, 95-114. 

De Jong, G.J. (1957). Application of Stress Functions to 
Consolidation Problems, Proc. 4th Int. Oonf, SMFE, 

London, 1957, Vol.1, pp. 320-323, 

Denvar, H. (1981). A Method of Settlement Calculation. 

Proc. 10th Int. Gonf. SWE, Stockholm, Sweden, June 1981, 
Vol. 2, pp. 99-102. 

Desai, C.S. (1968). Solution of Stress-Deformation Problems 
in Soil and Rock Mechanics Using the Finite Element 
Methods. Ph.D. Dissartation, Univ. Texas, Austin. 

Desai, C.S. (1972). Nonlinear Analysis Using Spline Functions. 

J. Sl'IIFS Div., ASCE, Vol. 98, Mo. SM9,pP.967-971 . Also 1971, 
Vo 1 . 97 , SMFE , pp . 1 461 -1 480 . 

Desai, C.S. (1977) inDesai^ O.S. and Christian, J.T. (Eds). 

Numerical Methods in Geotechnical Engineering. McGraw 
Hill Book Co., Mew York, 1977. 



291 


Besai, G.S, and Abel, J.?, (1972). Introduction to the 
Finite Element Method. Van ifestrand Reinhold Go., 
l\few York. 

Desai, C.S. and Oliristian, J.T. (Eds.) (1977), Uumerical 

.Methods in Geomechanical Engineering, McGraw Hill, 1977. 

Desai, C.S., Johnson, I.D. and Hargett, O.M’, ,(1974), Analysis 
of Pile Supported Gravity Lock. J. Geotech. Bngrg, Div. , 
ASGB, Vol, 100, Ho. GT9/pp. 1009-1029. 

Desai, G.S. and Reese, L.G. (1970), Stress Analysis of 

Gircular Footings on Layered Soils. J. Soil Mech. and 
Pound. Div., ASGE, Vol. 96, Ho. SM4, pp.1 289-1 310, 

Dhaliwal, R.S. (1970). Punch Problem for an Elastic layer 

Overlying an Elastic Foundation. Int . J. Engrg, Sci . , 

Vol. 8, Ho. 4, pp. 275-288. 

Dhaliwal, R.S, and Singh, B.M. (1977)-. Annular Punch on an 
Elastic Layer Overlying an Elastic Foundation. Int, 

J. Engrg. Sci., V0I.I5, Ho. 4, pp. 263-270. 

Duncan, J.M. and Ghang, G.Y. (1970). Honlinear Analysis of 

Stress and Strain in Soils, J. S115FE Div,, ASGE, Vol. 96, 

Ho, SM5, pp. 1629-1 655. 

Dunlop, P« and Duncan, J.M. (1970), Development of Failure Around 
Excavated Slopes. J, SMPE Div. aSGE, Vol. 96, SM2, 

pp. 471-493. 

Egorov, K.E. (1965). Calculation of Bed for Foundation with 
Ring Footing. Proc, 6th Int. Gonf. SivIIFS, Montreal, 

1965, V0I.2, pp. 41-45. 

Egorov, K.E., Barvashov, V.A. and Fedor ovsky, V.G, (1975). 

Some Applications of the Elasticity Theory to Design of 
Foundations. Proc, 8th Int. Gonf, SI'ffPE, Moscow, 1973, 
V0I.I, Ho. 5, pp. 69-74. 

Egorov, K.E., Kitaikina, O.V. and Zinovyev, A.V, (1979). _ 

Deformation of Ring Foundation Base Proc. 6th Asian 
Reg. Gonf. Sj^IPE, Singapore, 1979, 7ol.1, pp. 281-284. 

Egorov, K.E., Konovalov, P.A. , Kitaykina, O.V. , Salnikov, 

L.F. and Zinorev, A.V. (1977). Soil Deformations 
Under Gircular Footings. Proc, 9th Int. Gonf, SJ/IFE, 

Tokyo, Japan, July 1977, V0I.I, pp. 489-492, 



292 


Bgorov , K,E. and ITichiporovich, A. A. ( 196 I). Research 
on the Deflexion of Foundations. Proc, 5 th Int. Conf , 

SICPS, Paris, 1961, V 0 I.I, pp, 861-866. 

Egorov, K.E, and Serebrjanyi, R.7. (1963). Determination 
of Stresses in Rigid Circular Foundation. Proc. 2nd 
Asian Reg. Gonf . SMFS, Japan, 1963, T 0 I.I, pp. 246-250. 

Egorov, K.E,, Vronsky, A.V, , Finseva, f.V. and Illynykh, V.A, 
(1980). Study of the Stress and Strain State of a 
Subgrade Under an Eccentric load. Proc, 8 th Danube- 
European Conf. SMFE, Varna, Bulgaria, Sept. 1980, 

Vol. 2, pp. 67-78. 

England, A.H, (1962). A Punch Problem lor a Pransversly Isotropic 
layer, Proc. Gamb. Phil. Soc., Vol. 58 , No, 3, pp. 539-547. 

Eringen, A.G. (1962), Uonlinear Theory of Continuous Media. 
McGraw Hill, Hew York 1962. 

Gazetas, G, (I 98 I). Variational Estimation of the Settlement 
of a Gircular Raft on Anisotropic Soil, Soils and 
Foundations, Vol. 21, Ho. 4, pp. 109-116, 

Gazetas, G. (1982), Axis 3 niimetric Parabolic loading of 

Anisotropic Half -Space, J, Geotech. Engrg. Div, , 

Proc, ASGE, Vol. 108, Ho. GT4, pp. 654 - 660 . 

Ga*etas, G. (1982). Stresses and Displacements in Cross- 

Anisotropic Soil. J. Geotech, Engrg. Div., Proc, ASGE, 
Vol. 108, Ho. GT4, pp. 532-533. 

German! s, E. and Vallippan, S. (1975). lining Subsidence at 
the Graving Dock Site, Hew Castle. Proc, Symp._Recent 
Dev. in the Analysis of Soil Behaviour and their 
Applications to Geotech, Struct,, Univ. H.S.V/, , 

Australia. 

Gerrard O.M, and Harrison, W.J. (1971). Circular loads Applied 
To a Gross Anisotropic Half -Space. Tech. Paper Ho. 8 , 

Div. Appl. Geomech. , GSIRO, Australia. 

Gerrard, G.M. and Wardle, l.J. (1973). Solutions for Point 
loads and Generalised Circular loads Applied to a 
Cross Anisotropic Half-Space. Tech. Paper Ho,13j Div. 
Appl. Geomech. , GSIRO, Austr^alia. 



293 


Grhaboussi, J. and Yfilson, E.L, (1975). Flow of Incompressible 
Fluid in Porous Elastic Media. Int. J. Humber Meth. 
Engrg., Vol.5, pp. 419-442. 

Gibson, R.E. and Lumb, P. (1953). numerical Solution of Same 
Problems in the Consolidation of Olay, Proc, Instn. 

Civ. Engrs., Vol. 2, Hb.1, pp. 182-198. 

Gibson, R.E. , Schiffman, R.L. and Pu, S.l. (1970), Plane 

Strain and Axially Symmetric Consolidation of a Olay 
Layer on a Smooth Impervious Base. Quart. J. Mech. 

Appl. Math; Vol. 25 , No. 4, pp. 505-520. 

Girijavallabhan, C.V. and Reese, L.O. (I 968 ). Finite Element 
Method Applied to Some Problems in Soil Mechanics. 

J. SmE Qiv. ASOE, Vol. 94, Ro. SM2, pp. 473 - 496 . 

Girond , J.P, (I 972 ), Tables Pour le Oaicul des Foundations. 
Dunod, Paris, V 0 I.I. 

Gladwell, G.M. and Iyer, K.R.P. (1974). Unbonded Contact 

Between a Circular Plate and an Elastic Half -Space. 

J.E. Elasticity, Vol. 4, Ro.2, pp. 115-130. 

Hain, S.J. (1975). A. Study of Structure Eoundation Interaction. 
Ph.I). Thesis, Deptt. Civ. Engrg., UHSYT, Aust. 

Hanuska, A. and Hbvotny, B. (1978). Rigid Circular Plate on^ 
Multilayered Half -Space. Stavebricky Cas, 26 (1978), 

No. 516, pp. 383-392. 

Kara, T., Shibuya, T., Koizumi, T. and Nakahara, I. (1975). _ 

An Asymmetric Mixed B.V. Problem of the Thick Elastic 
Plate: Moment Applied to the Annular Rigid Stamp. 

Proc. 25 th Japan National Conf, for Appl. Mech, 1975, 
Theoretical and Applied Mechanics, pp, 505-512, Univ. 
of Tokyo Press. 

Harr, M.E. (1966). Foundations of Theoretical Soil Mechanics. 
McGraw-Hill Book Go., New York. 

Hooper, J.A. (1973). Observations on the Behaviour of a Piled 
Raft Foundation on London Clay. Proc. Inst, Giv, Engrs., 
pt.2, Vol. 55, p. 855. 

Hooper, J.A. (1974). Analysis of a Gircular Raft in Adhesive 
Contact with a Thick Elastic Layer, Geot echni que , 

Vol, 24 , No. 4, pp. 561-580, 



294 


Hooper, J.A. (1975). Elastic Settlement of a Circular Raft in 
Adhesive Contact with a iPransversly Isotropic Medium. 
Geotechnique, 7ol.25, Ro.4, PP. 691-711. 

Hooper, J.A. (1976). Parabolic Adhesive loading of a 

Flexible Raft Foundation. Geotechnique, Vol.26, No, 5, 

PP. 511-525. 

Huang, Y.H, (1969). Finite Element Analysis of Nonlinear 

Soil Media. Proc. Application of FEM in Civ. Engrg., 
Yanderbit Uhiv. , pp. 663-690. 

Hwang, O.P., Morgenstern, V/.R. and Murray, D.Y/. (1971). On 
Solutions of Plane Strain Consolidation Problems by 
FEM. Ganad. Geotech. J. , Yol.S, No.l, pp, 109-118. 

Hwang, C.T., Morgenstern, N..R. and Murray, D.W. (1972). 

Application of FEM to Consolidation Problems. Proc. 

Symp. Appl. FEM Geotech. Engrg., Yicksburg, ilississippi, 
1972, Y0I.2, pp. 739-765. 

Ishkova, A.G. (1951). Bending of a Circular Plate on the 
Elastic Half -Space Under the Action of a Uniform 
Axisymraetric Load, Uch. Zap. Mosk, Gos. Univ,,3, 
pp. 152-160. 

Keer, L.M, (I967), Mixed Bound ary-Yalue Problems for an 

Elastic Half -Space . Proc. Gamb. Phil. Soc., Yol.63> 
pp. 1379-1586. 

Kodner, R.L. (1963). Hyperbolic Stress-Strain Response : Cohesive 
Soils. J. SMFEDiv., ASGE, Yol.89, No. SMI, pp. 115-143. 

Komornik, a. and Mazurik, A, (1977). Restrained Settlements of 
Masonry Buildings, Sbh Int. Oonf. SMFE, lokyo, Japan, 
Proc., Yol.l, pp. 615-618, 

Lee, I.K. (Ed.) (l974). Soil Mechanics-New Horizons. Newnes- 
Butterworths, London. 

Leonards, G.A. (1962). Foundation Engineering. McGraw-Hill 
Book Co . N. Y. 

MadhaVjM.R. (1980). Settlement and Allowable Pressures for^ 

Rin45 or Annular Footings. Indian Getech. J. , Y0I.IO, 

No. 3, pp. 267-271 . 



295 


Madhav, M,R, and Vliikar, P,P. (1980 • Response of Finife 
Saturated Layer to Construction !I? 3 rpe of Loading, 

Int. Symp, Mech, Behaviour of Structural Media, 

Ottawa, Proc. Yol.2, pp. 355-367. 

McRamee, J, and Gibson, R,B. (1960). Plane Strain and Axially 
Syrmetric Problems of the Oonsolidation of a Semi- 
Infinite Olay Stratum. Quart, J. Mech. Appl. Math., 
Vol,13, 110 . 2 , pp. 210-227. 

Mandel, J. (1957). Oonsolidation des couches d'argiles, 

Proc. 4th Int. Conf. SWE, London, 1957, Vol.l, 
pp. 360-367. 

Mandel, J, (1961). Tassements Products par la Consolidation 
d ’une couche d'argile de grande epaisseur. Proc, 5th 
Int. Conf. SMPE, Paris, 1961, Vol.l, pp. 733-736. 

Mandel, J. and Salencon, J. (1969), Force Portante d ’un sol 

Sur une assise rigide. Proc. 7th Int. Conf. SIvIFE, Mexico, 
Vol.2, pp. 157-164. 

Milovic, D.M. (1973). Stresses and Displacements Produced by 
a Ring Foundation. Prof, 8th Int. Conf, Soil ifech. and 
Found. Engrg. , Moscow, Vol.l 3, pp. 167-179. 

Mindlin, R.D. (1936). Force at a Point in the Interior of a 
Semi-Infinite Solid. J, Appl, Phys., Vol,7, i1o.5, 
pp. 195-202. 

Mossakovski, V.I. (1954). The Fundamental General Problem 
of the Theory of Elasticity for a Half-Space v^ith 
a Circular Curve Determining Boundary Conditions, 

Prikl. Mat. Mekh,, Vol.18, pp. 187-196 (in Russian). 

Murray, R.T. (1978) in Scott, C.R. (Ed). Developments in 
Soil Mechanics-1. Applied Science Publishers Ltd,, 
London, 1978. 

Narasirahan, S.L.(1982). Finite Element Analysis of Piled 
Foundations. Ph.D. Thesis (under submission), 

Civ. Engrg. Deptt,, I.I.T. Kanpur. 



296 


ish-ida, Y» (1966), Verliical S'tress and Yen'tical DefoiMna1/ion.s 
of G-round Under a Deep Circular Uniform Pressure in 
the Semi-infinite. Press, 1st Gong. Int . Soc, Rock Mech,, 
Vol.2, pp. 493-498. 

Uayak, U.C. and Zienkiewicz, 0,0. (1972). Slasto -Plastic 

Stress Analysis. A Generalization for Various Consti- 
tutive Relations Including Strain Softening, Int. J, 
Uumer, Meth, Bngrg., Vol,5, pp. 113-135. 

Odemark, U, (1949). Investigation as to the Elastic Properties 
aiid Soils and Design of Pavements According to the 
Theory of Elasticity, Statens Vacginstitut , Stockholm. 

Oden, J.T, (1972). Finite Element Methods for Uonlinear 
Gontinua. McGraw-Hill Book Go,, New York. 

O'Neill, M.W. and Reese, L.G. (1972). Behaviour of Bored Piles 
in Bearaout Olay. J. SMEE Div. , ASGE, Vol.98, No.SM2, 

Eeb. 1972. 

Paria, G. (1957)» Axisymmetric Consolidation for a Porous 
Elastic Material Containing Fluid. J, Maths, Phys. , 

Vol. 36, pp. 338-346. 

Parry, R.H.G. (ed. ) (1971). Stress-Strain Behaviour of Soils, 
Prod, Roscoe Mem. Symp, Cambridge Univ., G.T.Foulis and 
Go. London, 1971 . 

Patankar, M.v. (1982), Non-linear Analysis of Annular and 

Circular Foundations by FEM, J, Instn. Engrs, (India), 
Vol. 62, Pt. Cl 6, pp. 347-552. 

Peck, R.B., Hanson, W.E. and Thornburn, T.H. (1974). Foundation 
Engineering. Wiley Eastern Ltd., New Delhi, 1980, 

Popova, O.V. (1972). Stress and Displacement Distributions in 
a Homogeneous Half -space Below a Circular Foundation, 
Soil Mech. Fdn. Engrg., Vol. 9, No, 2, pp, 86-89. 

Poulos, H.G. (1968). The Behaviour of a Rigid Circular Plate 
Resting on a Finite Elastic layer. Civ. Engrg. Trans. 
Instn. Engrs. Australia, OE 10(2), pp, 213-219. 

Poulos, H.G. and Davis, E.H.(l974). Elastic Solutions for 
Soil and Rock Mechanics, Wiley, New York, 



297 


Poulos, H.&, and Davis^ B.H. (1980). Pile Poundation Analysis 
and Design. John Y/iley and Sons, lew York, 1980. 

Ramberg, W. and Osgood, W.R. (1943). Description of Stress- 
Strain Curves by Three Parameters, latl. Advis. Comm. 
Aeronant,, Tech., I'fote 902, Washington, D.G. 

Rauhans, D. (l977). Tables for Calculating Circular Slabs on 
Elastic Support with Centrally S 3 nnmetric loading. 
Bauingenieur 52 , lo.lO, pp. 387-392. 

Roscoe, K.H, (1970), The Influence of Strains in Soil Mechanics. 

Tenth Rankine Lecture. G-eotechnique, Tol.2o, pp, 129-170. 

Sandhu, R.S. (1977). Variational Principle for Finite Element 

Analysis of Consolidation. Proc. Bum. Meth, in Ceomech,, 
2nd Int. Conf. .lum. Meth. Virginia, June 1976. 

Sandhu, R.S, and 1/Tilson, B.L. (1969), Finite Element Analysis 
of Seepage in Elastic Media. J, Engrg. Mech. Div. , 

Proc. ASCE, Vol.95, lo. EM3, pp. 641-652. 

Sanr^os, J. and Quera, J. (1981). Soil Foundation Interaction 
With Friction at the Interface, Proc, 10th Int, Conf. 
SMFE, Stockholm, Vol.2, pp. 253-258. 

Schiffman, R.L. and Aggarwala, B.D. (1961), Stresses and 

Displacements Produced in a Semi-infinite ' Elastic Solid 
by a Rigid Elliptical Footing. Proc. 5th Int. Conf. 

SMFE, Paris, 1961, Vol.1, pp. 795-801. 

Schultze, E. (l961). Distribution of Stress Beneath a Rigid 
Foundation. Proc, 5th Int. Conf. SMFE, Vg 1.6, 
pp. 807-813. 

Scott, C.R. (Ed.) (1978), Developments in Soil Mechanics-1. 
Applied Science Publishers ltd., London, 1978. 

Scott, R.F, ( 1963 ). Principles of Soil Mechanics, Addison, 
Yfelsey Publ. Co. Inc, 

Selvadurai, A.P.S. (1979). Elastic Analysis of Soil Foundation 
Interaction. Elsevier Scientific Publising Go,, 
Amsterdam, 1979} PP. 223-406. 

Selvadurai, A.P.S. (I 98 O). Flexure of a Thick Circular Raft 
Foundation Resting on an Isotropic Elastic Medium. 

Proc, Int. Conf . 'Structural Foundation in Rock* , 

Sydney, Australia, May 1980. 



298 


Selvadurai, A,P,S, (1980) • Ths Eccentric loadiiag of a Rigid 
Circular Foundation Embedded in an Isotropic Elastic 
Medium. Int. J. taer. Analyt, Methods, Yol.4, ]S'o,2, 
pp. 121..129. ’ ^ ’ 

Selvadurai, A.P.S, and Kenpthorne, R.H. (I98O), Plane Stradn- 
Contact Stress Distribution Beneath a Rigid Footing 
Resting on a Soft Cohesive Soil. Canad. Geotech. J., 
?ol.17, lfo.1, pp. .114-122. 

Shelest, L. A. (1975). Distribution of Stresses and Displacements 
in a Foundation Bed of Finite Thickness Beneath a 
Circular Rigid Foundation. Soil Mech. Fdn, Bngrg., 

V0I.I2, ¥0.6, pp. 404-407. 

Shibuya, T. Koizumi, T. atid Rakahara, I. (1974). An Elastic 
Contact Problem for a Half -space Indented by a Flat 
Annular Rigid Stamp. Int. J, Bngng. Sci.,Yol,l2, Ko,9, 
pp. 759-771. 

Shibuya, T,, Koizumi, T. and Hakahara, I, (1975). An Elastic 
Contact Problem for a Half -Space Indented by a Flat 
Annular Rigid Stamp in the Presence of Adhesion, Proc, 
25th Japan Rational Gong, for Appl, Mech, 1975. 
Theoretical and Applied Mechanics, pp, 497-504, Univ. 
of Tokyo Press. 

Sloan, S.Yif, and Randolph, M.F, (l982), Ifcmerical Prediction 

of Collapse Loads Using Finite Element, Int, J. Rumer, 
Analyt. Meth, Geomech., V0I.6, ifo.1, pp, 45-76. 

Smith, I.M. (1970). A Finite Element Approach to Ela.stic Soil 
Structure Interaction. Canad. Geotech, J. Vol,7 , pp. 
95-105. 

Smith, I.M. and Hobbs, R. (l976). Biot Analysis of Consolidation 
Beneath Embankments, Geotechnique , Yol,26, Iifo,1, pp. 

149-171, 

Sneddon, I.R. (1946), Boussinesq's Problem for a Flat Ended 
Cylinder. Proc. Camb. Phil. Soc, Y0I.42, Ifo.l, 

pp. 29-39. 

Spence, D.A. (1968). Self Similar Solutions to Adhesive 

Contact Problems with Incremental loading. Proc, Roy, 

Soc, A 305, pp. 55-80. 

Svec, O.J, (1974), The Unbounded Contact Problem of a Plate 
on the Elastic Half Space. Computer Meth. in Appl. 

Mech, Engrg, , Yo^,3, PP. 105-113. 



Szabo, G-. (1 970 / 71 ). Die Drimdlagen einer neuen. Festigfciet- 
■tneone, T 0 I.I and 2^ Bauverlag, Wiesbaden. 

Taylor, D.W, (1962). Fundamentals of Soil Mechanics. John 
Wiley and Sons, New York. 

Terzaghi, K. and Peck, R.B, (1948). Soil Mechanics in 
Engineering Practice. John Y/iley and Sons, 

New York. 

Timoshenko, S.P. and Krieger, S,W, (1959). Theory of 

Plates and Shells, McGrraw-Hill Kogakusha Ltd,, Tokyo. 

Valliappan, S. , Hain, S. and Lee, I.K. (Eds. ) (1975 ). Soil 
Mechanics-Jlecent Developments, Proc. Gen, Session 
Symp. Univ, NSW, Australia, July 14-18,19755. 

Valliappan, S, and lee, I.K, (1975). Gonsolidation of Nan- 
homogeneous Anisotropic Layered Soil Media, Proc, 

2nd Australia-NZ Oonf. Geomech, , Brisbane, pp, 67-71. 

Valov, G.M, (1968), Infinite Elastic Layer and Half -Space 

Under the Action of a Ring-Shaped Die, J, Appl, Math. 
Mech. , V 0 I. 52 , No. 5, pp. 917-930. 

Varadarajan, A. and Arora, K.R. (1982), Interaction of 

Circular Footing on Sand Bed System, 4th Int, Oonf. 
Numer. Meth. Geomech,, Edmonton, Vol.2, pp. 945-954. 

Varma, B.S,, Khadilkar, B.S, and (lhandrasekaran, V.S, (1977). 
Circular Footings on Transversly Isotropic Soils, 

Proc, Int, Symp, Soil Struct, Interaction, Roorkee, 
1977, pp. 265-270. 

Verruijt, A, (1977) in Finite Elements in Geomechanics, 

G. Gudehus (Ed.) John Wiley and Sons, London, 1977. 

Volterra, E. (1953). Circular Beams on an Elastic Foundation. 

J. Appl. Mech. Trans. ASMS, Vol.20, No. 2, pp. 227 -2 32. 

Wardle, L.J, (1977). A Computer Program for the Analysis of 

Multiple Complex Circular loads on Layered Anisotropic 
Media. User’s Manual, GSIRO-Div. -appl. Geomech, 
Program Circly, Geomech. Computer Program No. 2, 1977. 

Whitaker, T. and Cooke, R.W. (1966), An Investigation of 
Shaft £md Base Resistances of Large Bored Piles in 
London Olay. Proc. Symp. on Large Bored Piles, 
pp. 7-49. 


300 


Wilson, E.L, (1965). Structural Analysis of Axi symmetric 
Solids. J. Am. Inst, Aero. Astro., Tol.3, H'o.12, 
pp. 2269-2274* 

Winterkorn, H.P. and Pang, H.Y. (Eds.) (1975). Poundation 
Engineering Handbook. Van Ifc strand Reinhold Co., 

New York, 1975. 

Wong, I.H, (1971). Analysis of Braced Excavations, Sc.D, 
Thesis, Deptt. Civ, Engrg, MIT, Cambridge, 

Wood, L.A. (1977). The Economic Analysis of Raft Poundation. 

Int. J. Numer. Analyt. Meth, Geomech,, Y0I.I (1977)-, 
^0.4, pp. 397-405. 

Yamaguchi, H, (1973). Consolidation and Heaving of a Pinite 
Clay Layer, Proc. 8th Int. Conf, SMPB, Moscow, 1973. 

Yegorov, K.E. and Nichiporovich, A. A, (1961) (See Egorov, 

K.E. and Nichiperovich, A.a, (I96 I)). 

Yokoo, Y., Yamagata , K. and Nagaoka, H, (1971). Pinite 

Element Analysis of Consolidation Pollowing Undrained 
Deformation, Soils Pdns, , Vol.11, No, 4, pp, 37-58. 

Yokoo, Y., Yamagata, K, and Nagaoka, H. (1971). Pinite 

Element Method Applied to Biot's Consolidation Theory, 
Soils and Pdns, Vol.ll, N0.I, pp, 29-46, 

Yokoo, Y, , Yamagata, K, and Nagaoka, H. (1971). Variational 
Principles for Consolidation, Soils Pdns,, Vol.ll, 

No. 4, pp. 25-35. 

Zbirohowski-Koscia, K.P. and Gunasekera, D.A, (1970), 

Poundation Settlement and Ground Reaction Calculations 
Using a Digital Computer, 

Zeevaert, L, (1957). Foundation Design and Behaviour of 

Tower Latino Americana in Mexico City, Geotechnique, 
V0I.7, N0.I, p. 115. 

Zeevaert, L, (1972). Poundation Engineering for Difficult 
Subsoil Conditions. Van llbstrand Reirihold Go .Ltd., 

New York. 



Zienkiewicz, 0,0, (1977). The Finite Element Method . Tata McGraw- 
Hill Publ.Go.Ltd. ,¥ew Delhi, 1977. 

Zienkiewicz, 0,0,, Talliappan, S. and King, I.P. (1968). 

Stress Analysis of Rock as a *Ho-Tension’ Material. 
Geotechnique, Yol.18, pp. 56-66. 

Zinovev, A.H, (1979). Determination of the Deformation of 

a Base of Finite Thickness Under an Annular Foundation. 
Translated from Osnovaniya, Fundamenty i Mekhanika 
Gruntov, Ho. 3, PP. 23-25, May-June, 1979. 



